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PREFACE 


This  report  was  prepared  for  the  Flight  Dynamics  Directorate  of  the  Wright 
Laboratory  and  the  Foreign  Aerospace  Science  and  Technology  Center,  Wright-Patterson 
AFB,  Ohio  under  Contract  F33657-89-D-2198.  It  documents  Missile  Datcom  Revision 
6/93.  The  work  was  performed  by  the  McDonnell  Douglas  Aerospace,  St  Louis, 
Missouri,  a  division  of  the  McDonnell  Douglas  Corporation.  The  period  of  performance 
was  August  1992  to  January  1993.  Joseph  W.  Herrmann  and  William  B.  Blake  served 
as  the  Air  Force  project  engineers.  This  report  supersedes  WL  TR-9 1-3039,  produced 
under  Contract  F33615-87-C-3604,  which  documents  Missile  Datcom  Rev  4/91. 


A  list  of  the  Missile  Datcom  principal  investigators  and  individuals  who  made 
significant  contributions  to  the  development  of  this  program  is  provided  below. 


PRINCIPAL  INVESTIGATORS: 
Steven  R.  Vukelich 
Stan  L.  Stoy 
Andrew  A.  Jenn 
Kurt  D.  Bausch 
John  E.  Williams 
Keith  A.  Burns 

CONTRIBUTORS: 

Jerry  E.  Jenkins 
Frank  M.  Sawyer 
Charles  B.  Heath 
William  B.  Blake 
Richard  D.  Samuels 
Joseph  W.  Herrmann 
Kevin  D.  Bruns 
Marvin  E.  Moore 
Kenneth  J.  Deters 


DTCC  QUALITY  INSPECTED  3 


The  Flight  Dynamics  Directorate  is  committed  to  the  continuing  development  of 
Missile  Datcom.  The  development  is  dependent  to  a  large  extent  on  users’  feedback. 
Questions  about  the  program  or  suggestions  for  future  improvements  to  the  program 
should  be  directed  to  Mr.  William  Blake,  WL/FIGC,  2210  Eighth  St  Ste  21,  Wright- 
Patterson  Air  Force  Base,  OH  45433-7531,  phone  (513)  255-6764. 


BACKGROUND  TO  MISSILE  DATCOM 


Missile  Datcom  Revision  6/93  is  the  sixth  in  a  series  of  releases. 

In  the  late  1970’s,  the  Tri-Service  committee  on  Missile  and  Projectile 
Aerodynamics  defined  the  need  for  a  Missile  Datcom  type  prediction  tool.  The  Air  Force 
was  chosen  as  the  lead  service  for  the  effort.  A  contract  was  subsequently  awarded  to 
the  McDonnell  Douglas  Astronautics  Company  (F33615-80-C-3605)  to  recommend 
specific  methods  for  inclusion  into  a  potential  computer  program  and  identify  areas 
where  further  work  was  needed.  The  final  report  from  this  effort,  'Development 
Feasibility  of  Missile  Datcom,"  (AFWAL-TR-81-3130)  was  published  in  October  1981. 

In  September  1981,  The  Missile  Datcom  Development  Contract,  (F33615-81-C- 
3617)  was  awarded  to  the  McDonnell  Douglas  Astronautics  Company.  It  subdivided  the 
effort  into  four  distinct  phases.  The  initial  release  of  the  program  in  August  1984 
presented  the  "Phase  I"  interim  capability.  Cases  run  using  this  version  were  limited 
to  axisymmetric  bodies  with  no  more  than  eight  fins  total  (two  sets  with  up  to  four  fins 
each). 


The  second  release  of  the  program  (Rev  11/85)  represented  the  "Phase  IV" 
capability.  This  was  the  final  version  generated  under  contract  F33615-81-C-3617.  It 
added  capability  for  elliptic  bodies,  inlets  at  supersonic  speeds,  dynamic  derivatives, 
experimental  data  substitution,  and  configuration  incrementing.  It  also  increased  the 
permissible  number  of  fins  to  32  (4  sets  with  up  to  8  fins  each).  Two  volumes  of 
documentation  (User’s  Manual  and  Program  Implementation  Guide),  dated  November 
1985,  were  printed. 

The  third  release  (Rev/  12/88)  coincided  with  the  publishing  of  AFWAL  TR-86- 
3091  This  version  expanded  the  experimental  data  substitution  option  and  dynamic 
derivative  capability.  Error  from  the  11/85  version  were  also  corrected.  Volume  I  of  TR- 
86-3091  (Final  Report)  discusses  the  methods  selected  for  incorporation  into  Missile 
Datcom.  Volume  II  (User’s  Manual)  is  an  updated  version  of  the  November  1985 
manual.  These  reports  are  available  from  DTIC  as  ADA-211086  and  ADA-210128 
respectively. 

The  fourth  release  (Rev  7/89)  added  little  new  capability;  its  primary  purpose  was 
to  correct  coding  errors,  expand  the  body-alone  dynamic  derivative  capability,  and 
modify  the  equivalent  angle-of-attack  formulation.  No  new  documentation  was 
published. 

The  fifth  release  (Rev  4/91)  followed  contracted  efforts  with  Nielsen  Engineering 
and  Research  (F33615-86-C-3626)  and  the  McDonnell  Douglas  Missile  Systems  Company 
(F33615-87-C-3604).  This  version  expanded  the  inlet  capability  to  subsonic  speeds, 
added  methods  for  inlet  additive  drag,  base  plume  effects,  and  body  protuberances. 


Many  other  methods  were  improved  and  coding  errors  were  corrected.  The  User’s 
Manual  for  the  4/91  revision,  WL-TR-91-3039,  reflected  these  changes  and  is  available 
from  DTIC  as  ADA-237817. 

The  6/93  version  of  the  program  is  the  end  product  an  effort  to  further  extend  the 
capability  of  the  code  and  expand  the  user’s  manual.  This  work  was  performed  by 
McDonnell  Douglas  Aerospace  under  contract  F33657-89-D-2198.  The  major  changes 
to  the  program  and  associated  reference  material  are  listed  below: 

(a)  The  program  was  modified  to  be  compatible  with  Unix  based  workstations. 

(b)  A  method  for  trailing  edge  flap  effectiveness  has  been  added.  The  method  is  valid 
at  subsonic  (ref.  1)  and  supersonic  (ref.  2)  speeds.  A  cubic  fairing  between  these 
results  is  used  in  the  transonic  region  (0.8  <  Mach  <  1.4). 

(c)  The  fin-body  carryover  interference  calculation  has  been  expanded  to  include  the 
effect  of  fin  dihedral  (ref.  3). 

(d)  The  inlet  method  has  been  expanded  to  treat  semi-submerged  shapes. 

(e)  The  base  drag  calculation  for  boattails  and  flares  has  been  modified. 

(f)  Coding  errors  have  been  corrected. 

(g)  The  information  in  the  original  Program  Implementation  Guide  (printed 
November  1985)  has  been  updated  and  included  as  part  of  this  User’s  Manual. 


vi 
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1.0  INTRODUCTION 


In  missile  preliminary  design  it  is  necessary  to  quickly  and 
economically  estimate  the  aerodynamics  of  a  wide  variety  of  missile 
configuration  designs.  Since  the  ultimate  shape  and  aerodynamic 
performance  are  so  dependent  upon  the  subsystems  utilized,  such  as  payload 
size,  propulsion  system  selection  and  launch  mechanism,  the  designer  must 
be  capable  of  predicting  a  wide  variety  of  configurations  accurately.  The 
fundamental  purpose  of  Missile  Datcom  is  to  provide  an  aerodynamic  design 
tool  which  has  the  predictive  accuracy  suitable  for  preliminary  design,  and 
the  capability  for  the  user  to  easily  substitute  methods  to  fit  specific 
applications. 
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2.0  PROGRAM  CAPABILITIES 


The  computer  code  is  capable  of  addressing  a  wide  variety  of 
conventional  missile  designs.  For  the  purposes  of  this  document,  a 
conventional  missile  is  one  which  is  comprised  of  the  following: 

•  An  axisymmetric  or  elliptically-shaped  body 

•  One  to  four  fin  sets  located  along  the  body  between  the  nose 
and  base.  Each  fin  set  can  be  comprised  of  one  to  eight 
identical  panels  attached  around  the  body  at  a  common 
longitudinal  position 

•  An  airbreathing  propulsion  system. 

To  minimize  the  quantity  of  input  data  required,  commonly  used 
values  for  many  inputs  are  assumed  as  defaults.  However,  all  program 
defaults  can  be  overridden  by  the  user  in  order  to  more  accurately  model  the 
configuration  of  interest. 

The  following  paragraphs  detail  the  configurations  that  can  be 
analyzed.  Later  paragraphs  describe  the  range  of  aerodynamic  coefficients  that 
can  be  predicted.  Finally,  the  program  constraints  are  discussed. 

2.1  ADDRESSABLE  CONFIGURATIONS 

The  following  configurations  can  be  analyzed: 

•  Circular  or  elliptically-shaped  cross  section  bodies,  with  or 
without  airbreathing  inlets 

•  Fin  alone  (1  to  8  panels  attached  at  the  root) 

•  Body  and  up  to  4  fin  sets  (1  to  8  panels  in  each  fin  set) 

•  The  body  and  fin  set  configurations  with  deflected  fins 

Certain  restrictions  exist  due  to  method  limitations  and  are 
summarized  in  the  following  paragraphs. 
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2.1.1  Axisymmetric  or  Elliptical  Bodies 


Methodology  is  incorporated  that  permits  analysis  of  the  configuration 
components  summarized  in  Table  1.  Due  to  the  types  of  methods  selected 
restrictions  also  apply  to  the  manner  in  which  these  components  are  joined 
to  form  a  complete  configuration: 

•  Subsonic/transonic  speeds  -  The  aerodynamic  methods 
assume  that  the  body  is,  as  a  minimum,  composed  of  a  nose- 
cylinder  combination.  The  afterbody  (boattail  or  flare)  is 
optional,  but  if  used,  it  must  be  attached  to  a  cylindrical  center 
body  whose  length  is  at  least  four  body  diameters;  this 
restriction  minimizes  nose  flow  field  coupling  over  the 
afterbody.  If  an  afterbody  is  specified  it  must  not  be 
cylindrical,  e.g.,  the  base  diameter  must  be  different  than  the 
centerbody  diameter.  Table  2  summarizes  the  other 
restrictions  on  the  configurations. 

•  Supersonic  speeds  -  The  aerodynamic  methods  used  are  not 
restricted  to  nose-cylinder  combinations.  Any  arbitrary  radii 
distribution  can  be  defined  since  theoretical  techniques  are 
employed  at  Mach  numbers  above  1.2.  Care  should  be  taken 
to  avoid  introducing  unexpected  comers  into  the  contour.  If 
the  contour  has  any  concaved  regions  the  marching  may  fail 
due  to  shock  impingement  on  the  body  as  it  starts  to  curve 
out. 

2.1.2  Panels 

The  program  will  accept  inputs  to  describe  most  airfoil  sections  or 
planforms.  Certain  assumptions  and  limitations  are  made  and  summarized 
in  the  following  paragraphs. 

2. 1.2.1  Airfoil  Section  -  The  program  will  accept  virtually  any  symmetrical 
airfoil  section  or  NASA  subsonic  cambered  section.  The  airfoil  section  can  be 
defined  using  a  NACA  designation  or  by  supplying  the  coordinates  of  the 
section.  Circular  arc,  hexagonal,  or  diamond  shaped  sections  can  also  be 
specified.  A  symmetric  hexagonal  cross-section  is  the  default;  its  shape  is 
computed  using  the  planform  inputs.  Hence,  explicit  definition  of  the  airfoil 
section  is  optional.  Although  cambered  airfoil  sections  can  be  input,  their  use 
in  the  code  is  currently  limited  to  subsonic  applications. 
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2. 1.2.2  Planform  -  Each  set  of  fins  may  be  comprised  of  up  to  eight  separate 
panels.  It  is  assumed  that  each  panel  is  geometrically  identical.  Although 
planforms  may  be  described  by  up  to  10  separate  pieces  or  sections,  an 
equivalent  straight-tapered  panel  is  computed  and  used  at  all  speeds.  There  is 
no  capability  to  specify  a  panel  with  outboard  dihedral. 

2.1.3  Airbreathing  Inlets 

Both  axisymmetric  and  two-dimensional  airbreathing  inlet /diverter 
combinations  can  be  defined.  Up  to  20  identical  inlets  can  be  positioned 
around  the  body  at  arbitrary  angles.  Vehicles  with  inlets  can  be  analyzed  at  all 
speeds. 


2.2  TYPES  OF  DATA  COMPUTED 
2.2.1  Aerodynamics 

The  program  computes  the  following  aerodynamic  parameters  as  a 
function  of  angle  of  attack  for  each  configuration: 

Cn  Normal  Force  Coefficient 

Cl  Lift  Coefficient 

Cm  Pitching  Moment  Coefficient 

XcP  Center  of  Pressure  in  calibers  from  the  moment 

reference  center 

Ca  Axial  Force  Coefficient 

Cq  Drag  Coefficient 

Cy  Side  Force  Coefficient 

Cn  Yawing  Moment  Coefficient 

Q  Rolling  Moment  Coefficient 

Cn«  Normal  force  coefficient  derivative  with  angle  of 

attack 

Cma  Pitching  moment  coefficient  derivative  with  angle 

of  attack 

Cyp  Side  force  coefficient  derivative  with  sideslip  angle 

Cnp  Yawing  moment  coefficient  derivative  with 

sideslip  angle 

Cip  Rolling  moment  coefficient  derivative  with 

sideslip  angle 

The  derivative  output  can  be  in  degrees  or  radians.  Partial  output 
results,  which  detail  the  components  used  in  the  calculations,  are  also 
optionally  available. 
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It  should  be  ncted  that  the  drag  force  (and  drag  coefficient)  is  different 
between  the  wind  and  stability  axes  systems  if  the  missile  body  is  at  a  sideslip 

angle  (f3)  to  the  wind.  However,  wind  axis  drag  and  stability  axis  drag  are  the 
same  at  zero  sideslip.  In  Missile  Datcom,  drag  force  methods  are  assumed  to 
be  in  the  stability  axes  system  and  axial  force  methods  are  assumed  to  be  in 
the  body  axes  system  unless  otherwise  noted. 

The  program  has  the  capability  to  perform  a  static  trim  of  the 
configuration,  using  any  fin  set  for  control  with  fixed  incidence  on  the  other 
sets.  The  two  types  of  aerodynamic  output  available  from  the  trim  option  are 
as  follows: 

•  Untrimmed  data  -  Each  of  the  aerodynamic  force  and 
moment  coefficients  are  printed  in  a  matrix,  which  is  a 
function  of  angle  of  attack  and  panel  deflection  angle.  This 
output  is  optional. 

•  Trimmed  data  -  The  trimmed  aerodynamic  coefficients,  and 
trim  deflection  angle,  are  output  as  a  function  of  angle  of 
attack. 

2.2.2  Geometry 

All  components  of  the  configuration  have  their  physical  properties 
calculated  and  output  for  reference  if  requested.  All  data  is  supplied  in  the 
user  selected  system  of  units. 

2.2.3  Other 

The  reference  area  and  reference  length  are  user  defined.  The  user  may 
optionally  select  to  print  the  calculated  body  or  fin  pressure  coefficient 
distributions  at  supersonic  speeds.  Outputs  of  the  partial  aerodynamic  results 
and  a  summary  of  method  extrapolations  are  also  optionally  available. 

2.3  OPERATIONAL  CONSIDERATIONS 

The  code  has  been  written  to  conform  to  the  coding  standards  for  the 
American  National  Standards  Institute  (ANSI)  Standard  FORTRAN  X3.9- 
1978,  often  referred  to  as  FORTRAN  V  or  FORTRAN  77. 

There  are  only  four  exceptions  to  the  ANSI  standards  used  in  the  code: 

•  Transfer  on  end-of-file  -  The  FORTRAN  IV  statement 
IF(EOF(UNIT).NE.O)  is  used  in  the  standard  CDC  compatible 
code.  The  FORTRAN  V  statement  END=labeI  in  the 
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FORTRAN  READ  statements  is  incorporated  in  the  code,  but 
is  inactive  in  the  CDC  code. 

•  NAMELIST  -  The  use  of  namelist  for  input  and  output  (I/O) 
is  used.  Although  this  appears  to  be  a  violation  of  FORTRAN 
IV,  it  is  really  not  since  a  namelist  emulator  has  been  written 
for  the  Missile  Datcom  code  using  FORTRAN  IV. 

•  Mathematical  functions  -  A  few  mathematical  functions  are 
not  considered  "standard,"  such  as  the  trigonometric  tangent. 
Standard  FORTRAN  equivalents  for  these  functions  are 
available  on  request. 

•  PROGRAM  card  -  The  code  was  developed  on  the  Control 
Data  Corporation  CYBER  computers.  This  system  requires 
that  the  first  card  of  the  main  routine  be  a  PROGRAM 
statement.  An  IBM  or  VAX  compatible  version  of  the  code  is 
also  maintained  which  has  a  different  format  for  the  program 
card. 
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3.0  INPUT  DEFINITION 


Inputs  to  the  program  are  grouped  by  "case".  A  "case"  consists  of  a  set 
of  input  cards  which  define  the  flight  conditions  and  geometry  to  be  run. 
Provisions  are  made  to  allow  multiple  cases  to  be  run.  The  successive  cases 
can  either  incorporate  the  data  of  the  previous  case  (using  the  input  card 
SAVE)  or  be  a  completely  new  configuration  design.  The  SAVE  feature,  for 
example,  permits  the  user  to  define  a  body  and  wing  (or  canard)  configuration 
in  the  first  case  and  vary  the  tail  design  for  subsequent  cases. 

The  scheme  used  to  input  data  to  the  computer  program  is  a  mixture  of 
namelist  and  control  cards.  This  combination  permits  the  following: 

•  Inputs  are  column  independent  and  can  be  input  in  any 
order. 

•  All  numeric  inputs  are  related  to  mnemonic  (variable) 
names. 

•  Program  input  "flags"  are  greatly  reduced.  Required  "flags" 
are  identified  by  a  unique  alphabetic  name  which 
corresponds  to  the  option  selected. 

The  program  includes  an  error  checking  routine  which  scans  all  inputs 
and  identifies  all  errors.  This  process  is  a  single-pass  error  checking  routine; 
all  errors  are  identified  in  a  single  "run”.  In  addition,  the  program  checks  for 
necessary  valid  inputs,  such  as  a  non-zero  Reynolds  number.  In  some  cases, 
the  code  will  take  corrective  action.  The  type  of  corrective  action  taken  is 
summarized  later  in  this  section. 

Flexibility  has  been  maintained  for  all  user  inputs  and  outputs.  The 
following  summarize  the  program  generality  available: 

•  The  units  system  can  be  feet,  inches,  meters  or  centimeters. 

The  default  is  feet. 

•  Derivatives  can  be  expressed  in  degree  or  radian  measure. 
Degree  measure  is  the  default. 

•  The  body  geometry  can  be  defined  either  by  shape  type  or  by 
surface  coordinates. 

•  The  airfoil  can  be  user  defined,  NACA,  or  supersonic  shaped 
sections.  The  NACA  sections  are  defined  using  the  NACA 
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designation.  A  hexagonally  shaped  supersonic  section  is  the 
default. 

•  The  configuration  can  be  run  at  a  fixed  sideslip  angle  and 
varying  body  angle  of  attack,  or  a  fixed  aerodynamic  roll  angle 
and  varying  total  angle  of  attack. 

•  The  flight  conditions  can  be  user  defined,  or  set  using  a 
Standard  Atmosphere  model.  The  capability  to  define  wind 
tunnel  test  conditions  as  the  flight  conditions  is  also 
available.  The  default  flight  condition  is  zero  altitude. 

3.1  NAMELIST  INPUTS 

The  required  program  inputs  use  FORTRAN  namelists.  Missile 
Datcom  is  similar  to  other  codes  which  use  the  namelist  input  technique,  but 
differ  as  follows: 

•  Namelist  inputs  are  column  independent,  and  can  begin  in 
any  column  including  the  first.  If  a  namelist  is  continued  to  a 
second  card,  the  continued  card  must  leave  column  1  blank. 

Also,  the  card  before  the  continued  card  must  end  with  a 
comma.  The  last  usable  column  is  number  79  if  column  1  is 
used,  and  column  80  if  column  1  is  blank. 

•  The  same  namelist  can  be  input  multiple  times  for  the  same 
input  case.  The  total  number  of  namelists  read,  including 
repeat  occurrences  of  the  same  namelist  name,  must  not 
exceed  300. 

The  three  namelist  inputs 

$REFQ  SREF=1.,$ 

$REFQ  LREF=2.,$ 

$REFQ  ROUGH=0.001,$ 

are  equivalent  to 

$REFQ  SREF=1  .,LREF=2.,ROUGH=0.001  ,$ 

•  The  last  occurrence  of  a  namelist  variable  in  a  case  is  the 
value  used  for  the  calculations. 


The  three  namelist  inputs 


$REFQ  SREF=1.,$ 

$FLTCON  NMACH=2.,MACH=1 .0,2.0,$ 

$REFQ  SREF=2.,$ 

are  equivalent  to 

$REFQ  SREF=2.,$ 

$FLTCON  NMACH=2.,  MACH=1.0/  2.0,$ 

•  The  namelists  can  be  input  in  any  order. 

•  Only  those  namelists  required  to  execute  the  case  need  be 
entered. 

•  Certain  hollerith  constants  are  permitted.  They  are 
summarized  in  Table  3.  Note  that  any  variable  can  be 
initialized  by  using  the  constant  UNUSED;  for  example, 
LREF=UNUSED  sets  the  reference  length  to  its  initialized 
value. 

All  Missile  Datcom  namelist  inputs  are  either  real  numbers  or  logical 
constants.  Integer  constants  will  produce  a  nonfatal  error  message  from  the 
error  checking  routine  and  should  be  avoided. 

The  namelist  names  have  been  selected  to  be  mnemonically  related  to 
their  physical  meaning.  The  ten  namelists  available  are  as  follows: 


Namelist 

Inputs 

$FLTCON 

Flight  Conditions  (Angles  of  attack, 
Mach  numbers,  etc.) 

$REFQ 

Reference  quantities  (Reference  area, 
length,  etc.) 

$AXIBOD 

Axisymmetric  body  definition 

$ELLBOD 

Elliptical  body  definition 

$PROTUB 

Protuberance  information  and 
geometry 

SFINSETn 

Fin  descriptions  by  fin  set  (n  is  the  fin 
set  number;  1,  2,  3  or  4) 

$DEFLCT 

Panel  incidence  (deflection)  values 

$TRIM 

Trimming  information 

$INLET 

Inlet  geometry 

$EXPR 

Experimental  data 
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Each  component  of  the  configuration  requires  a  separate  namelist 
input.  Hence,  an  input  case  of  a  body-wing-tail  configuration  requires  at  least 
one  of  each  of  the  following  namelist  inputs,  since  not  all  variables  have 
default  values  assigned: 

$FLTCON  to  define  the  flight  conditions 

$AXIBOD  or  $ELLBOD  to  define  the  body 
$FINSET1  to  define  the  most  forward  fin  set 

$FINSET2  to  define  the  first  following  fin  set 

$FINSET3  to  define  the  second  following  fin  set 

$FINSET4  to  define  the  third  following  fin  set 

The  following  namelists  are  optional  since  defaults  exist  for  all  inputs: 

$REFQ  to  define  the  reference  quantities 

$PROTUB  to  define  protuberance  option  inputs 

$DEFLCT  to  define  the  panel  incidence 

(deflection  angles) 

$TRIM  to  define  a  trim  case 

$INLET  to  define  inlet  geometry 

$EXPR  to  define  experimental  input  data 

Defaults  for  all  namelists  should  be  checked  to  verify  the  configuration 
being  modeled  does  not  include  an  unexpected  characteristic  introduced  by  a 
default. 

The  following  sections  describe  each  of  the  namelist  inputs.  Each 
section  is  accompanied  by  a  figure  which  summarizes  the  input  variables, 
their  definitions,  and  units.  Since  the  system  of  units  can  be  optionally 
selected,  the  column  "Units"  specifies  the  generic  system  of  units  as  follows: 

L  Units  of  length;  feet,  inches,  centimeters  or  meters 
F  Units  of  force;  pounds  or  Newtons 
deg  Units  of  degrees;  if  angular,  in  angular  degrees;  if 
temperature,  either  degrees  Rankine  or  degrees 
Kelvin 

sec  Units  of  time  in  seconds 

Exponents  are  added  to  modify  the  above.  For  example,  L2  means  units 
of  length  squared,  or  area.  Combinations  of  the  above  are  also  used  to  specify 
other  units.  For  example,  F/L2  means  force  divided  by  area,  which  is  a 
pressure. 

Since  it  is  difficult  to  discern  the  difference  between  the  number  zero 
"0"  and  the  alphabetic  letter  "O",  it  should  be  noted  that  none  of  the  namelist 
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or  namelist  variable  names  contain  the  number  zero  in  them.  In  general,  the 
number  zero  and  the  letter  "O"  are  not  interchangeable  unless  so  staicd. 


The  program  ascertains  the  configuration  being  modeled  by  the 
presence  of  each  component  namelist,  even  if  no  data  is  entered.  The 
following  rules  for  namelist  input  apply: 

•  Do  not  include  a  namelist  unless  it  is  required.  Once  read,  the 
presence  of  a  namelist  (and,  hence,  a  configuration 
component)  can  only  be  removed  using  the  DELETE  control 
card  in  a  subsequent  case.  Simply  setting  all  variables  to  their 
initialized  values  will  not  remove  the  configuration 
component. 

•  Do  not  include  a  variable  within  a  namelist  unless  it  is 
required.  Program  actions  are  often  determined  from  the 
number  and  types  of  input  provided. 

•  Do  not  over-specify  the  geometry.  User  inputs  will  take 
precedence  over  program  calculations.  Inputs  that  define  a 
shape  that  is  physically  impossible  will  be  used  as  specified. 

The  program  does  not  "fix-up"  inconsistent  or  contradictory 
inputs. 

3.1.1  Namelist  FLTCON  -  Flight  Conditions 

This  namelist  defines  the  flight  conditions  to  be  run  for  the  case.  The 
program  is  limited  to  no  more  than  20  angles  of  attack  and  20  Mach  numbers 
per  case  at  a  fixed  sideslip  angle,  aerodynamic  roll  angle,  altitude,  and  panel 
deflection  angle.  Therefore,  a  "case”  is  defined  as  a  fixed  geometry  with 
variable  Mach  number  and  angles  of  attack. 

The  inputs  are  given  in  Figure  1.  There  are  two  ways  in  which  the 
aerodynamic  pitch  and  yaw  angles  can  be  defined: 

•  Input  ALPHA  and  BETA.  If  BETA  is  input  and  PHI  is  not,  it  is 
assumed  that  the  body  axis  angles  of  attack  (a)  and  sideslip 
angles  (p)  are  defined. 

•  Input  ALPHA  and  PHI.  If  PHI  is  input  and  non-zero,  it  is 
assumed  that  ALPHA  is  the  total  angle  of  attack  (a)  and  PHI 
is  the  aerodynamic  roll  angle  (q). 

•  Input  ALPHA,  BETA  and  PHI.  The  value  for  BETA  is  ignored 
if  PHI  is  non-zero. 
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As  a  minimum  the  following  variables  must  be  defined: 

NALPHA  number  of  angles  of  attack  to  run 

(N  ALPHA  >  2) 

ALPHA  angle  of  attack 

NALPHA) 

NMACH  number  of  Mach 

(NMACH  >  1) 

MACH  or  VINF  Mach  number  or  speed  schedule  (matching 

NMACH) 

The  REN,  TINF  and  PINF  data  must  correspond  to  the  MACH  or  VINF 
inputs.  The  ALPHA  and  MACH  dependent  data  can  be  input  in  any  order;  the 
code  will  sort  the  data  into  ascending  order. 

Reynolds  number  is  always  required.  Three  types  or  inputs  are 
permitted  to  satisfy  the  Reynolds  number  requirement: 

•  Specify  Reynolds  number  per  unit  length  using  REN 

•  Specify  the  altitude  using  ALT,  and  the  speed  using  MACH  or 
VINF  (Reynolds  number  is  computed  using  the  Standard 
Atmosphere  model) 

•  Specify  pressure  and  temperature  using  PINF  and  TINF,  and 
the  speed  using  MACH  or  VINF  (typical  of  data  available 
from  a  wind  tunnel  test) 

User  supplied  data  will  take  precedence  over  program  calculations. 
Hence,  the  user  can  override  any  default  or  Standard  Atmosphere  calculation. 
The  default  condition  is  sea- level  altitude  (ALT=0.)  if  the  wrong  combination 
of  inputs  are  provided  and  the  Reynolds  number  cannot  be  calculated. 

3.1.2  Namelist  REFQ  -  Reference  Quantities 

Inputs  for  this  namelist  are  optional  and  are  defined  in  Figure  2.  A 
vehicle  scale  factor  (SCALE)  permits  the  user  to  input  a  geometry  that  is 
scaled  to  the  size  desired.  This  scale  factor  is  used  as  a  multiplier  to  the  user 
defined  geometry  inputs;  it  is  not  applied  to  the  user  input  reference 
quantities  (SREF,  LREF,  LATREF).  If  no  reference  quantities  are  input,  they 
are  computed  based  upon  the  scaled  geometry.  XCG  is  input  relative  to  the 
origin  (X=0)  and  is  scaled  using  SCALE. 

In  lieu  of  specifying  the  surface  roughness  height  ROUGH,  the  surface 
Roughness  Height  Rating  (RHR)  can  be  specified.  The  RHR  represents  the 
arithmetic  average  roughness  height  variation  in  millionths  of  an  inch. 
Typical  values  of  ROUGH  and  RHR  are  given  in  Table  4. 


schedule  (matching 
numbers  or  speeds 
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3.1.3  Namelist  AXIBOD  -  Axisymmetric  Body  Geometry 

An  axisymmetric  body  is  defined  using  this  namelist.  The  namelist 
input  variables  are  given  in  Figures  3a  and  3b  and  a  sketch  of  the  geometric 
inputs  are  given  in  Figure  4.  The  body  can  be  specified  in  one  of  two  ways: 

OPTION  1:  The  geometry  is  divided  into  nose/  centerbody,  and 
aft  body  sections.  The  shape,  overall  length,  and  base  diameter 
for  each  section  are  specified.  Note  that  not  all  three  body 
sections  need  to  exist  on  a  configuration;  for  example,  a  nose- 
cylinder  configuration  does  not  require  definition  of  an  aft  body. 

OPTION  2:  The  longitudinal  stations  and  corresponding  body 
radii  are  defined,  from  nose  to  tail.  This  option  should  only  be 
selected  if  the  Mach  number  is  greater  than  1.2. 

If  Option  2  is  selected,  the  program  generates  a  bodv  contour  based  on 
the  user  specified  values  of  X,  R,  and  DISCON.  Many  additional  points  in 
between  the  user  specified  input  coordinates  will  be  generated.  The  resulting 
contour  can  contain  more  than  1000  points.  If  the  PRINT  GEOM  BODY 
control  card  is  used,  this  contour  will  be  written  to  tape  unit  3. 

The  program  uses  the  input  value  for  NX  to  determine  which  option  is 
being  used.  If  NX  is  not  input  then  Option  1  inputs  are  assumed.  If  both 
shapes  and  body  coordinates  (Options  1  and  2)  are  used,  the  body  coordinate 
information  will  take  precedence.  NX  can  be  set  to  its  initialized  value  (to 
simulate  the  variable  as  not  input)  by  specifying  "NX=UNUSED". 

It  is  highly  recommended  that  Option  1  be  used  when  possible.  The 
program  automatically  calculates  the  body  contour  based  upon  the  segment 
shapes  using  geometry  generators.  Hence,  more  accurate  calculations  are 
possible.  Even  when  Option  2  is  used,  appropriate  Option  1  inputs  should  be 
included.  This  identifies  where  the  code  should  insert  break  points  in  the 
contour.  If  these  parameters  are  not  input,  they  are  selected  as  follows: 

LNOSE  Length  of  the  body  segment  to  where  the 
radius  first  reaches  a  maximum 
DNOSE  The  dianu  er  at  the  first  radius  maximum 
LCENTR  Length  of  the  body  segment  where  the  radius 
is  constant 

DCENTR  Diameter  of  the  constant  radius  segment 
LAFT  The  remaining  body  length 

DAFT  Diameter  at  the  base 

DEX1T  Not  defined  (implies  that  base  drag  is  not  to 
be  included  in  the  axial  force  calculations) 
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If  DEXIT  is  not  input,  or  set  to  UNUSED,  the  base  drag  computed  for 
the  body  geometry  will  not  be  included  in  the  final  computed  axial  force 
calculations.  To  include  a  "full"  base  drag  increment,  a  zero  exit  diameter 
must  be  specified  (DEXIT=0.). 

If  body  coordinates  are  input  using  the  variables  NX,  X,  R,  and 
DISCON,  and  the  nose  is  spherically  blunted,  the  geometry  must  be 
additionally  defined  using  the  following: 

•  BNOSE  must  be  specified  (even  if  zero) 

•  TRUNC  must  be  set  to  .FALSE. 

•  The  first  five  (5)  points  in  the  X  and  R  arrays  must  lie  on  the 
spherical  nose  cap  (i.e.,  X(l),  X(2),  X(3),  X(4),  X(5),  R(l),  R(2), 

R(3),  R(4),  and  R(5)  are  spherical  cap  coordinates] 

The  following  summarizes  the  input  generality  available: 

•  X(l)  does  not  have  to  be  0.0;  an  arbitrary  origin  can  be 
selected. 

•  Five  shapes  can  be  specified  by  name: 

CONICAL  (CONE)  -  cone  or  cone  frustrum  (default  for 
boattails  and  flares) 

OGIVE  -  tangent  ogive  (default  for  noses) 

POWER  -  power  law* 

HAACK  -  L-V  Haack  (length-volume  constrained)* 

KARMAN  -  von  Karman  (L-D  Haack;  length-diameter 
constrained)* 

•  If  DAFT<DCENTR  the  afterbody  is  a  boattail.** 

•  If  DAFT>DCENTR  the  afterbody  is  a  flare.** 

•  If  LAFT  is  not  input,  aft  body  (boattail  or  flare)  does  not  exist. 

•  applies  to  noses  only 

**  DAFT  must  not  be  equal  to  IX  r  ‘  TK 

The  inputs  for  base-jet  plum*  :  r.unon  effects  are  defined  using 
Option  1.  Incremental  forces  and  moments  due  to  jet  induced  boattail 
separation  and  separation  locations  on  aft  fins  are  calculated  if  these  inputs 
are  used. 
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This  option  should  only  be  run  for  supersonic  cases 
(i.e.  Moo  >  1.2) 


•  The  calculations  will  be  done  for  three  types  of  aft  bodies 
conical  boattail,  ogival  boattail,  or  cylindrical  (i.e.  no  boattail). 

Error  messages  will  be  printed  to  the  output  file  and  the 
calculations  skipped  if  any  other  aft  body  is  defined. 

•  If  BASE=. FALSE,  or  is  not  input  the  calculations  will  be 
skipped. 

•  DEXIT  must  not  equal  zero  if  this  option  is  used. 

•  The  jet  Mach  number  (JMACH),  jet  to  freestream  static 
pressure  ratio  (PRAT),  and  jet  to  freestream  stagnation 
temperature  ratio  (TRAT)  must  be  specified  for  each 
freestream  Mach  number  or  velocity  input  in  the  namelist 
FLTCON.  For  subsonic  or  transonic  freestream  Mach 
numbers  or  velocities,  dummy  values  must  be  input  for 
JMACH,  PRAT,  and  TRAT.  The  user  must  be  careful  to 
match  these  inputs  with  the  proper  freestream  conditions. 

•  If  a  portion  of  the  fins  in  a  fin  set  are  located  on  the  boattail  or 
base,  the  boattail  separation  locations  will  be  calculated  and 
output  at  each  fin  roll  angle.  However,  if  the  fins  do  not 
extend  to  the  boattail  the  separation  locations  will  be  skipped. 

•  Results  may  be  inaccurate  if  excessive  extrapolation  is 
required.  If  extrapolation  occurs,  a  warning  message  will  be 
printed  to  the  output  file.  To  avoid  extrapolation  and 
minimize  inaccuracy,  the  input  parameters  should  be  kept 
within  the  ranges  shown  in  Figure  5. 

3.1.4  Namelist  ELLBOD  -  Elliptical  Body  Geometry 

Elliptically-shaped  cross  section  bodies  are  defined  using  this  namelist. 
The  inputs  are  similar  to  those  for  the  axisymmetric  body  geometry 
(AXIBOD),  and  are  shown  in  Figures  6a  and  6b.  The  types  of  shapes  available, 
and  the  limitations,  are  the  same  as  those  given  for  axisymmetric  bodies. 
However,  the  base-jet  plume  interaction  input  options  in  namelist  AXIBOD 
are  not  available  in  namelist  ELLBOD.  Please  read  Section  3.1.3  for 
limitations. 

Note  that  the  body  cross  section  ellipticity  can  vary  along  the  body 
longitudinal  axis.  Sections  which  are  taller-than-wide  and  wider-than-tall  can 
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be  mixed  to  produce  "shaped"  designs.  The  shape  of  the  sections  is  controlled 
by  the  variables  ENOSE,  ECENTR,  and  EAFT  or  ELLIP,  H  and  W. 

3.1.5  Namelist  PROTUB  -  Protuberance  Geometry 

Missile  protuberances  can  be  input  using  this  namelist.  Axial  force 
coefficient  is  calculated  for  the  protuberances  and  added  to  the  body  axial  force 
coefficient.  Figure  7  shows  the  inputs  required.  Figure  8  shows  the  different 
protuberance  shapes  available.  The  following  defines  the  inputs  required  for 
protuberance  calculations: 

•  NPROT  is  the  number  of  protuberance  sets.  A  protuberance 
set  is  made  up  of  protuberances  at  the  same  axial  location 
with  the  same  size  and  shape.  Therefore,  it  is  only  necessary 
to  describe  the  geometry  of  one  individual  protuberance  per 
set.  The  maximum  number  of  protuberance  sets  is  20. 

•  NLOC  is  the  number  of  protuberances  in  each  protuberance 
set.  NLOC  accounts  for  the  number  of  identical  protuberances 
located  around  the  missile  body  at  a  given  axial  location. 

•  The  following  equation  helps  to  clarify  the  relationship 
between  NLOC  and  NPROT: 

NLOC(l)+NLOC(2)+NLOC(3)+  . .  .  +NLOC(NPROT)  =  (Total 
number  of  protuberances  on  the  missile) 

•  The  axial  location  of  a  protuberance  (XPROT)  should  be  input 
at  the  protuberance  geometric  centroid.  An  approximation  of 
the  centroid  will  be  adequate  for  the  analysis.  The  location  is 
used  to  calculate  the  average  boundary  layer  thickness  over 
the  protuberance  length. 

•  VCYL,  HCYL,  BLOCK,  and  FAIRING  type  protuberances  have 
1  member.  LUG  types  have  4  members  and  SHOE  types  have 
3  members.  (Refer  to  Figure  8) 

•  All  inputs  for  LPROT,  WPROT,  HPROT,  and  OPROT  are  in 
sequential  order  based  upon  the  members  specified  with  the 
protuberance  type  (PTYPE)  input. 

•  The  FAIRING  type  protuberance  should  always  have  a  zero 
offset.  The  code  will  assume  a  zero  offset  even  if  a  non-zero 
offset  is  input. 


18 


More  complex  protuberance  shapes  can  be  analyzed  by  a  component 
build-up  method.  Each  member  is  treated  as  a  separate  protuberance. 
Combinations  of  vertical  cylinders,  horizontal  cylinders,  and  flat  plates  or 
blocks  can  be  input  at  specified  offsets  from  the  missile  body.  If  a  FAIRING 
type  protuberance  is  used  in  a  component  build-up,  the  offset  should  be  zero. 
The  user  must  manually  add  axial  force  of  the  individual  members  of  the 
component  build-up  if  the  total  protuberance  axial  force  is  desired. 

Figure  9  shows  an  example  input  file  for  a  missile  with  several 
protuberances. 

3.1.6  Namelist  FINSETn  -  Define  Fin  Set  n 

Figure  10a  describes  the  variables  needed  to  be  input  for  fin  set 
planform  geometry  descriptions.  Optional  fin  cross-section  inputs  are 
described  in  Figure  10b.  Special  user  specified  fin  cross-sections  can  be  input 
using  the  variables  in  Figure  10c.  The  user  may  specify  up  to  four  non¬ 
overlapping  fin  sets.  The  variable  "n"  in  the  namelist  specifies  the  fin  set 
number.  Fin  sets  must  be  numbered  sequentially  from  the  front  to  the  back  of 
the  missile  beginning  with  fin  set  one.  An  input  error  will  occur  if  "n"  is  zero 
or  omitted.  The  code  allows  for  between  1  and  8  geometrically  identical  panels 
to  be  input  per  fin  set.  The  panels  may  be  arbitrarily  rolled  about  the  body  and 
can  be  given  dihedral. 

Four  types  of  airfoil  sections  are  permitted-hexagonal  (HEX),  circular 
arc  (ARC),  NACA  airfoils  (NACA),  and  user  defined  (USER).  Only  one  type  of 
airfoil  section  can  be  specified  per  fin  set,  and  this  type  is  used  for  all  chord 
wise  cross  sections  from  root  to  tip.  Diamond-shaped  sections  are  considered  a 
special  case  of  the  HEX  type;  hence,  hexagonal  and  diamond  sections  can 
coexist  on  the  same  panel.  The  airfoil  proportions  can  be  varied  from  span 
station  to  span  station. 

The  user  selects  "break  points”  on  the  panel  (Figure  11a).  A  "break 
point"  specifies  a  change  in  leading  or  trailing  edge  sweep  angle.  Also  a  break 
point  may  specify  a  change  in  airfoil  section,  but  the  section  must  be  of  the 
same  type  (i.e.,  a  change  in  section  type  cannot  go  from  a  NACA  to  an  ARC) 
only  the  proportions  can  change.  The  location  of  each  "break  point"  is  defined 
by  specifying  its  semi-span  station  (SSPAN)  from  the  vehicle  centerline  and 
distance  from  the  first  body  station  to  the  chord  leading  edge  (XLE).  The 
"break  point"  chord  leading  edge  array  (XLE)  can  be  defined  by  simply 
specifying  the  root  chord  leading  edge  (XLE(l)]  and  the  sweep  angles  of  each 
successive  panel  segment  if  the  semi-span  stations  are  input.  Note  that  only 
those  variables  that  uniquely  define  the  fin  need  to  be  entered.  Redundant 
inputs  can  lead  to  numerical  inconsistencies  and  subsequent  computational 
errors. 
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The  panel  sweep  angle  (SWEEP)  can  be  specified  at  any  span  station  for 
each  segment  of  the  panels.  If  STA=0.,  the  sweep  angle  input  is  measured  at 
the  segment  leading  edge;  if  STA=1.,  the  sweep  angle  input  is  measured  at  the 
segment  trailing  edge.  Note  that  some  aerodynamic  methods  are  very 
sensitive  to  panel  sweep  angle.  For  small  span  fins,  small  errors  in  the 
planform  inputs  can  create  large  sweep  angle  calculation  errors.  It  is 
recommended  that  exact  sweep  angles  be  specified  wherever  possible;  for 
example,  if  the  panel  trailing  edge  is  unswept,  specifying  SWEEP=0.  and 
STA=1.  will  minimize  calculation  error.  Then  the  leading  edge  sweep  will  be 
computed  by  the  code  internally  using  the  SSPAN  and  CHORD  inputs. 

Plain  trailing  edge  devices  may  be  modelled  in  Missile  Datcom.  This  is 
accomplished  via  the  CFOC  array  which  is  the  flap  chord  to  fin  chord  ratio, 
Cf/c.  Trailing  edge  devices  can  be  either  full  span  or  partial  span  subject  to 
certain  limits  specified  below.  The  trailing  edge  devices  can  not  have  a  taper 
ratio  greater  than  1.0,  and  the  hinge  line  must  be  strsight  regardless  of  the 
number  of  segments  comprising  the  trailing  edge  device.  A  partial  span 
trailing  edge  device  is  specified  by  setting  CFOC=0  for  those  chord /span 
stations  that  are  not  part  of  the  trailing  edge  device.  Examples  of  acceptable 
and  unacceptable  geometries  are  shown  in  Figure  lib  as  well  as  the 
corresponding  input  values  for  the  variable  arrays  CFOC,  CHORD  and 
SSPAN.  A  special  case  where  the  trailing  edge  device  extends  to  the  tip  of  a 
fin  with  a  taper  ratio  of  zero  is  also  shown  in  Figure  lib.  While  any  value  of 
CFOC  will  result  in  the  correct  flap  chord  at  the  tip  (since  the  tip  chord  is 
zero),  the  user  must  specify  a  CFOC =1.0  since  a  value  of  CFOC=0  would 
indicate  the  trailing  edge  device  does  not  exist  at  this  chord/span  station,  the 
user  should  also  be  aware  of  the  following: 

•  All  trailing  edge  deflection  angles  are  measured  with  respect  to  the 
freestream  and  not  relative  to  the  hinge  line.  This  becomes  an 
important  distinction  as  the  hinge  line  sweep  angle  is  increased. 

•  The  variable  SKEW  does  not  apply  to  trailing  edge  devices. 

•  The  hinge  moments  for  trailing  edge  devices  are  not  calculated. 

•  The  increase  in  profile  drag  due  to  trailing  edge  deflection  is  not 
calculated. 

Since  all  panels  are  assumed  to  be  planar  (i.e.,  no  tip  dihedral),  all 
inputs  must  be  "true  view".  Once  the  planform  of  a  single  panel  is  defined, 
all  fins  of  the  set  are  assumed  to  be  identical.  The  number  of  panels  present  is 
defined  using  the  variable  NPANEL.  Each  panel  may  be  rolled  to  an  arbitrary 
position  around  the  body  using  the  variable  PHIF.  PHIF  is  measured 
clockwise  from  top  vertical  center  (looking  forward  from  behind  the  missile) 
as  shown  in  Figure  12.  Each  panel  may  also  contain  a  constant  dihedral.  A 
panel  has  zero  dihedral  when  it  is  aligned  along  a  radial  ray  from  the 
centerline  (see  Figure  12).  The  variable  used  to  specify  dihedral  is  GAM.  GAM 
is  positive  it  the  panel  tip  chord  is  rotated  clockwise. 
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Different  aerodynamics  will  be  computed  depending  upon  whether  the 
FLTCON  namelist  variable  PHI,  or  the  FINSETn  namelist  variable  PHIF,  is 
used  to  roll  the  geometry.  Figure  13  depicts  the  usage  of  the  roll  options.  The 
variable  "PHI"  means  that  the  body  axes  system  is  to  be  rolled  with  the 
missile  body,  whereas  PHIF  keeps  the  aerodynamics  in  a  non-rolled  body  axis, 
but  rather  locates  the  fin  positions  around  the  body.  PHIF  must  be  input  for 
each  panel,  while  PHI  rolls  the  whole  configuration. 

It  is  the  user's  responsibility  to  assure  that  the  fins  are  (1)  on  the  body 
surface,  and  (2)  do  not  lie  internal  to  the  body  mold  line.  The  program  does 
not  check  for  these  peculiarities.  If  SSPAN(1)=0  is  input,  the  program  will 
assume  that  the  panel  semi-span  data  relative  to  its  root  chord  are  supplied. 
The  code  will  automatically  interpolate  the  body  geometry  to  place  the  panel 
on  the  body  surface  with  the  root  chord  parallel  to  the  body  centerline.  See 
Section  3.4  for  modeling  fins  on  body  segments  of  varying  radii. 

When  defining  more  than  one  fin  set,  the  fin  sets  must  never  have 
their  pianforms  overlap  one  another.  There  must  be  sufficient  space  between 
the  forward  fin  trailing  edge  and  aft  fin  leading  edge  to  avoid  violating  the 
assumptions  made  by  the  aerodynamic  computations.  It  is  assumed  by  the 
aerodynamic  model  that  the  vortices  are  fully  rolled  up  when  they  pass  the 
control  points  of  the  next  downstream  set  of  fins.  In  reality  the  vortex  sheet 
does  not  fully  roll  up  until  it  is  at  least  four  semispans  downstream.  If  two  fin 
sets  are  closer  than  this  the  results  may  be  in  error  since  the  use  of  a  vortex 
filament  model  may  introduce  too  much  vorticity.  The  closer  the  spacing  the 
larger  the  error  may  be.  No  algorithm  error  will  result  from  too  close  a  fin  set 
spacing. 

Panels  with  cut-out  portions  can  be  modeled  by  using  one  of  the  ten 
available  fin  segments  as  a  transition  segment.  This  is  accomplished  by  giving 
the  segment  a  small  span,  such  as  0.0001,  and  specifying  the  segment  root  and 
tip  chords  to  transition  into  the  cut-out  portion  of  the  fin. 

3.1.7  Namelist  DEFLCT  -  Panel  Deflection  Angles 

This  namelist  permits  the  user  to  fix  the  incidence  angle  for  each  panel 
in  each  fin  set.  The  variables  are  given  in  Figure  14.  Note  that  the  panel 
numbering  scheme  is  assumed  to  be  that  shown  in  Figure  12.  The  array 
element  of  each  deflection  array  corresponds  to  the  panel  number. 

The  scheme  for  specifying  deflection  angles  is  unique,  yet  concise.  The 
scheme  used  is  based  upon  the  body  axis  rolling  moment: 

"In  Missile  Datcom  a  positive  panel  deflection  is  one  which  will 

produce  a  negative  (counterclockwise  when  viewed  from  the 
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rear)  roll  moment  increment  at  zero  angle  of  attack  and 

sideslip." 

3.1.8  Namelist  TRIM  -  Trim  Aerodynamics 

This  namelist  instructs  the  program  to  statically  trim  the  vehicle 
longitudinally  (Cm=0).  The  inputs  are  given  in  Figure  15.  Note  that  only  one 
fin  set  can  be  used  for  trimming.  The  user  only  specifies  the  range  of 
deflection  angles  desired  using  DELMIN  and  DELMAX;  the  code  will  try  to 
trim  the  vehicle  for  each  angle  of  attack  specified  using  the  allowable  fin 
deflections.  This  option  will  not  trim  the  vehicle  at  a  specific  angle  of  attack  if 
the  deflection  required  is  outside  the  range  set  by  the  values  of  DELMIN  and 
DELMAX. 

The  deflection  sign  convention  used  is  that  described  in  Section  3.1.7; 
hence,  DELMIN  and  DELMAX  are  input  as  if  deflecting  the  panel  to  the 
maximum  will  produce  a  negative  rolling  moment  from  the  panels  resulting 
normal  force  increment.  DELMIN  must  always  be  less  than  DELMAX. 

A  logical  variable,  ASYM,  has  been  included  to  permit  reverse  panel 
deflections.  For  example,  deflecting  all  panels  in  one  sense  results  in  a  rolling 
moment  and  no  pitching  moment.  The  ASYM  flag  will  permit  analysis  of  an 
elevator  (or  pitch  deflection)  effect,  by  deflecting  panels  on  one  side  of  the 
vehicle  only,  with  opposite  panels  mirroring  those  deflections.  Since  a 
maximum  of  eight  panels  are  allowed  in  a  fin  set,  only  four  panels  of  the  fin 
set  can  be  deflected  in  the  reverse  direction  using  the  ASYM  flag.  Both 
trimmed  and  untrimmed  results  are  available  for  output. 

3.1.9  Namelist  INLET  -  Axisymmetric  and  2-Dimensional  Inlet  Geometry 

This  namelist  is  used  to  model  the  inlet  and  diverter  geometry. 
Axisymmetric,  two-dimensional  side  mounted,  and  two-dimensional  top 
mounted  inlets  can  be  described.  The  inlets  may  be  covered  or  uncovered  and 
oriented  in  any  position  about  the  missile  body.  Inlet  normal  force,  pitching 
moment,  side  force,  yawing  moment,  and  axial  force  are  calculated.  The 
methods  are  valid  for  subsonic,  transonic,  and  supersonic  speeds.  Figure  16 
shows  the  INLET  namelist  inputs,  and  Figures  17a,  17b,  17c  and  17d  show  the 
inlet/ diverter  geometry  for  each  type  of  inlet.  The  inlets  may  have  a 
boundary  layer  diverter,  be  conformal  (diverter  height  HDIV=0),  or  be  semi- 
submerged  (diverter  height  HDIVcO).  The  methods  used  for  the  inlets  are 
the  same  regardless  of  wether  the  inlet  has  a  diverter  or  is  semi-submerged, 
and  they  are  not  applicable  to  chin  inlets.  The  variable  HDIV  is  used  to 
determine  whether  a  diverter  exist.  Figure  17d  shows  examples  of  two- 
dimensional  and  axisymmetric  inlets  that  are  conformal  or  semi-submerged. 


22 


•  Inlet  roll  orientation  uses  the  same  convention  as  the  fin 
panel  roll  orientation. 

•  Inlet  height  and  width  or  inlet  diameter  is  input  at  five  axial 
locations  described  in  Figures  17a,  17b,  and  17c: 

1)  leading  edge  or  tip 

2)  cowl  lip  leading  edge 

3)  midbody  start 

4)  boattail  start 

5)  boattail  end 

•  If  the  inlet  is  covered  (COVER=.TRUE.),  no  flow  is  allowed 
into  the  inlet.  The  inlet  is  plugged  between  stations  1  and  2, 
flush  with  the  inlet  face. 

Inlet  additive  drag  or  spillage  drag  can  be  calculated  for  external 
compression  inlets  operating  at  off-design  conditions  (Moo<Mdesign)  for  Mach 
numbers  greater  than  1.  Whenever  flow  spillage  occurs,  the  mass  flow  ratio  is 
less  than  one,  and  additive  forces  are  generated  on  the  deflected  streamtube 
captured  by  the  inlet.  If  the  inlet  operates  on-design,  the  ramp  shock  lies  on 
the  inlet  face  and  on  the  cowl  lip.  In  these  cases,  the  maximum  mass  flow 
ratio  is  one  (zero  spillage)  and  the  minimum  additive  forces  are  zero. 

•  If  the  inlet  is  covered  (COVER=.TRUE.),  the  additive  drag 
calculations  will  be  skipped. 

•  If  ADD=. FALSE.,  or  is  not  input  the  additive  drag  calculations 
will  be  skipped. 

•  Mass  flow  ratio  (MFR)  must  be  specified  for  each  freestream 
Mach  number  or  velocity  given  in  namelist  FLTCON.  For 
Mach  numbers  less  than  1,  dummy  values  must  be  inpuf  for 
MFR.  The  user  must  be  careful  to  match  these  inputs  with 
the  proper  freestream  conditions. 

•  The  additive  drag  is  calculated  at  zero  angle  of  attack  and 
assumed  to  remain  constant  for  all  angles  of  attack. 

3.1.10  Namelist  EXPR  -  Experimental  Data  Substitution 

This  namelist  is  used  to  substitute  experimental  data  for  the  theoretical 
data  generated  by  the  program.  The  variables  to  be  input  are  shown  in  Figure 
18.  Use  of  namelist  EXPR  does  not  stop  the  program  from  calculating 
theoretical  data,  but  rather  the  experimental  data  is  used  in  configuration 
synthesis,  and  it  is  the  experimental  data  that  is  used  for  the  component 
aerodynamics  for  which  it  is  input. 
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Experimental  data  may  be  substituted  for  any  configuration  component 
or  partial  configuration.  Experimental  data  is  input  at  a  specific  Mach 
number.  When  using  namelist  EXPR,  the  case  must  be  run  at  the  Mach 
number  for  which  you  are  substituting  experimental  data.  However,  the 
experimental  data  being  input  may  have  different  reference  quantities  and  a 
different  center  of  gravity  location  than  the  case  being  rim. 

Experimental  data  input  for  a  fin  alone  is  input  as  panel  data,  not  as 
total  fin  set  data.  The  user  should  note  that  experimental  data  for  fin  alone 
Cma  is  not  used  in  the  configuration  synthesis  process.  Instead  fin  alone  Cn« 
(the  experimental  value  if  input)  is  used  to  determine  the  fin  contribution  to 
Cma  during  configuration  synthesis.  If  body  alone  experimental  data  and  body- 
fin  experimental  data  are  input  for  the  same  case  the  body  data  is  ignored  in 
configuration  synthesis.  If  experimental  Cma  data  is  input  for  a  body  +  1  fin  set 
for  a  multi-fin  set  configuration,  the  calculated  contributions  to  Cma  of  the 
other  fin  sets  are  added  to  the  experimental  data. 

Since  the  experin  tal  namelist  forms  the  basis  for  configuration 
incrementing,  the  lateral  uirectional  coefficients  are  included  to  allow  for 
sideslip  cases.  These  coefficients  are  input  the  same  as  the  longitudinal 
coefficients.  However,  if  the  lateral  directional  coefficients  are  input,  the 
lateral  directional  beta  derivatives  will  not  be  computed  our  output. 

The  following  rules  apply  to  the  use  of  namelist  EXPR. 

•  It  is  assumed  that  the  coefficients  in  EXPR  are  for  the  same 
sideslip  and/or  aerodynamic  roll  as  the  case  being  run. 

•  Separate  namelist  EXPR  must  be  specified  for  each  Mach 
number. 

•  Each  namelist  EXPR  must  end  with  a  SEND  card. 

•  Separate  namelist  EXPR  must  be  specified  for  each  partial 
configuration  for  which  experimental  data  is  to  be  input,  (i.e., 
body,  body  +  1  fin  set,  etc) 

•  Separate  namelist  EXPR  must  be  specified  for  each  reference 
quantity  change. 

Example: 

The  user  has  experimental  data  available  for  a  body  +  2  fin  set 
configurations  and  is  interested  in  the  effects  of  adding  a 
booster  containing  a  third  fin  set.  he  would  then  use  namelist 
EXPR  to  input  the  experimental  data.  When  the 
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configuration  is  synthesized,  it  would  use  the  experimental 
data  for  body  +  2  fin  sets  and  theoretical  data  for  fin  set  three. 


3.2  CONTROL  CARD  INPUTS 

Control  cards  are  one  line  commands  which  select  program  options. 
Although  they  are  not  required  inputs,  they  permit  user  control  over 
program  execution  and  the  types  of  output  desired.  Control  cards  enable  the 
following: 

•  Printing  internal  data  array  results  for  diagnostic  purposes 
(DUMP) 

•  Outputting  intermediate  calculations  (PART,  BUILD, 
PRESSURES,  PRINT  AERO,  PRINT  EXTRAP,  PRINT  GEOM, 

PLOT,  NAMELIST,  WRITE,  FORMAT) 

•  Selecting  the  system  of  units  to  be  used  (DIM,  DERIV) 

•  Defining  multiple  cases,  permitting  the  reuse  of  previously 
input  namelist  data  or  deleting  namelists  of  a  prior  case 
(SAVE,  DELETE,  NEXT  CASE) 

•  Adding  case  titles  or  comments  to  the  input  file  and  output 
pages  (*,  CASEID) 

•  Limits  the  calculations  to  longitudinal  aerodynamics  (NO 
LAT) 

3.2.1  Control  Card  -  General  Remarks 

A  total  of  42  different  control  cards  are  available.  There  is  no  limit  to 
the  number  of  control  cards  that  can  be  present  in  a  case.  If  two  or  more 
control  cards  contradict  each  other,  the  last  control  card  input  will  take 
precedence.  All  control  cards  must  be  input  as  shown,  including  any  blanks. 
Control  cards  can  start  in  any  column  but  they  cannot  be  continued  to  a 
second  card.  Misspelled  cards  are  ignored.  Control  cards  can  be  located 
anywhere  within  a  case. 

Once  input,  the  following  control  cards  remain  in  effect  for  all 
subsequent  cases: 

DIM  FT  DIM  IN  DIM  CM  DIMM  FORMAT 

HYPER  INCREMENT  NOGO  NO  LAT  PLOT 

SOSE  WRITE 
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The  following  control  cards  are  effective  only  for  the  case  in  which  they 
appear: 


BUILD 
DUMP  CASE 
PRESSURES 
PRINT  GEOM 


CASEID  DAMP  DELETE 

DUMP  NAME  NAMELIST  PART 

PRINT  AERO  PRINT  EXTRAP 

SAVE  SPIN  TRIM 


These  control  cards  can  be  changed  from  case  to  case: 


DERIV  DEG  DERIV  RAD  NACA 


The  only  control  card  that  can  be  optionally  saved,  from  case-to-case,  is 
the  NACA  card. 


3.2.2  Control  Card  Definition 

Available  control  cards  are  summarized  as  follows: 
BUILD 


This  control  card  instructs  the  program  to  print  the  results  of  a 
configuration  build-up.  All  configurations  which  can  be  built  from  the 
components  defined  will  be  synthesized  and  output,  including  isolated  data 
(e.g.,  body  alone,  fin  alone,  etc.).  Component  build-up  data  is  not  provided  if 
the  TRIM  option  is  selected. 

CASEID 


A  user  supplied  title  to  be  printed  on  each  output  page  is  specified.  Up 
to  72  characters  can  be  specified  (card  columns  8  to  80). 

DAMP 


When  DAMP  control  card  is  input  longitudinal  dynamic  derivatives 
are  computed  and  the  results  output  for  the  configuration.  Dynamic 
derivatives  for  configuration  components  or  partial  configurations  may  be 
output  using  the  PART  or  BUILD  control  cards  respectively. 

DELETE  namel,name2 

This  control  card  instructs  the  program  to  ignore  a  previous  case 
namelist  input  that  was  retained  using  the  SAVE  control  card.  All  previously 
saved  namelists  with  the  names  specified  will  be  purged  from  the  input  file. 
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Any  new  inputs  of  the  same  namelist  will  be  retained.  At  least  one  name 
(namel)  must  be  specified. 


DIM  IN.  DIM  FT.  DIM  CM.  or  DIM  M 

This  control  card  sets  the  system  of  units  for  the  user  inputs  and 
program  outputs.  The  four  options  are  inches  (DIM  IN),  feet  (DIM  FT), 
centimeters  (DIM  CM),  and  meters  (DIM  M).  The  default  system  of  units  is 
feet.  Once  the  system  of  units  has  been  set,  it  remains  set  for  all  subsequent 
cases  of  the  "run". 

DERIV  DEG  or  DERIV  RAD 


All  output  derivatives  are  set  to  either  degree  (DERIV  DEG)  or  radian 
(DERIV  RAD)  measure.  The  default  setting  is  degree.  The  derivative  units 
can  be  changed  more  than  once  during  the  run  by  inputting  multiple  DERTV 
cards. 

DUMP  CASE 


Internal  data  blocks,  used  in  the  computation  of  the  case,  are  written  on 
Tape  6.  This  control  card  automatically  selects  partial  output  (PART). 

DUMP  namel  ,name2 

This  permits  the  user  to  write  selected  internal  data  blocks  or  common 
blocks  on  Tape  6.  At  least  one  name  (namel)  must  be  specified.  The  arrays 
will  be  dumped  in  units  of  feet,  pounds,  degrees  or  degrees  Rankine.  Table  5 
shows  the  common  block  dump  names  and  Tables  8  through  62  provide  a 
definition  of  each  common  block. 

FORMAT  (format) 

This  control  card  is  used  in  conjunction  with  the  WRITE  control  card. 
It  specifies  the  format  of  the  data  to  be  printed  to  tape  unit  3.  The  format  is 
input  starting  with  a  left  parenthesis,  the  form**  and  a  right  parenthesis.  This 
is  exactly  the  same  as  a  FORTRAN  FORMAT  statement.  Because  of  the  code 
structure,  alphanumeric  data  must  not  be  printed.  For  example: 

FORMAT  ((8(2X,F10.4))  is  legal 

FORMAT  (2HX=,F10.4)  is  illegal 

The  default  format  is  8F10.4,  and  will  be  used  if  the  FORMAT  control 
card  is  not  present.  Multiple  formats  can  be  used.  The  last  FORMAT  read  will 
be  used  for  all  successive  WRITE  statements  until  another  FORMAT  is 
encountered.  Hence,  the  FORMAT  must  precede  the  applicable  WRITE. 
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HYPER 


This  control  card  causes  the  program  to  select  ihe  Newtonian  flow 
method  for  bodies  at  any  Mach  number  above  1.4.  HYPER  should  normally  be 
selected  at  Mach  numbers  greater  than  6. 

INCRMT 

This  card  is  used  to  set  the  configuration  incrementing  flag. 
Configuration  incrementing  uses  the  first  case  of  a  run  to  determine 
correction  factors  for  the  longitudinal  and  lateral  aerodynamic  coefficients. 
These  correction  factors  are  computed  by  comparing  theoretical  and 
experimental  values  for  each  coefficient  for  which  data  is  input.  The 
experimental  values  are  input  using  namelist  EXPR.  During  subsequent  cases 
of  the  run,  the  correction  factors  are  applied  to  coefficients  for  which 
experimental  data  was  input  in  the  first  case.  This  provides  the  user  with  a 
method  to  evaluate  changes  in  a  configuration. 

The  INCRMT  card  must  be  input  in  the  first  case  of  a  nm.  The  first  case 
must  be  run  at  the  same  Mach  number  as  the  experimental  data  which  is 
input.  Once  the  increment  flag  is  set  it  cannot  be  deleted  during  that  run. 

The  following  restrictions  apply: 

•  All  cases  of  a  run  must  have  the  same  number  of  fin  sets. 

•  All  cases  of  a  run  must  have  the  same  sideslip  or 
aerodynamic  roll  angle  as  the  first  case  (BETA  or  PHI  as 
specified  in  namelist  FLTCON). 

•  The  first  case  must  be  run  at  exactly  the  same  angles  of  attack 
as  the  experimental  data  being  input. 

•  All  cases  must  be  run  within  the  same  Mach  regime 
(subsonic,  transonic,  or  supersonic)  as  the  experimental  data. 

•  Experimental  data  can  only  be  input  in  the  first  case  and  only 
for  the  complete  configuration.  No  additional  data  can  be 
substituted. 

•  To  increment  Cyp  and  Cno  experimental  data  must  be  input 
for  Cy  and  Cn 
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Use  of  configuration  incrementing  may  or  may  not  increase  the 
accuracy  of  the  results.  The  following  guidelines  will  produce  better  results 
when  using  configuration  incrementing: 

•  The  user  may  run  different  angles  of  attack  in  each  case. 
However,  no  angle  of  attack  should  exceed  the  upper  or 
lower  limit  of  the  angles  of  attack  for  which  experimental 
data  was  input  in  the  first  case. 

•  Experimental  data  should  be  input  at  as  many  angles  of  attack 
as  possible. 

•  The  user  should  remember  that  the  effect  of  a  change  in 
Mach  number  from  case  to  case  is  not  corrected  by  inputting 
experimental  data  at  one  Mach  number  as  is  required. 

NACA 

This  card  defines  the  NACA  airfoil  section  designation  (or  supersonic 
airfoil  definition).  Note  that  if  airfoil  coordinates  and  the  NACA  card  are 
specified  for  the  same  aerodynamic  surface,  the  airfoil  coordinate  specification 
will  be  used.  Therefore,  if  coordinates  have  been  specified  in  a  previous  case 
and  the  SAVE  option  is  in  effect,  the  saved  namelist  must  be  deleted  or  the 
namelist  variable  SECTVP  must  be  changed  for  the  NACA  card  to  be 
recognized  for  that  aerodynamic  surface.  The  airfoil  designated  with  this  card 
will  be  used  for  all  segments  and  panels  of  the  fin  set. 

The  form  of  this  control  card  and  the  required  parameters  are  as 
follows: 


Card  Column(s)  Input(s) 

Purpose 

1  thru  4 

NACA 

The  unique  letters  NACA 
designate  that  an  airfoil  is  to  be 
defined 

5 

Any  delimiter 

6 

1,2,3,  or  4 

Fin  set  number  for  which  the 
airfoil  designation  applies 

7 

Any  delimiter 

8 

1,4,5,6,S 

Type  of  NACA  airfoil  section; 
1-series  (1),  4-digit  (4),  5-digit  (5), 
6-series  (6),  or  supersonic  (S) 

9 

Any  delimiter 

10  thru  80 

Designation 

Input  designation  (see  Table  6); 
columns  are  free-field  (blanks 
are  ignored) 
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Only  fifteen  (15)  characters  are  accepted  in  the  airfoil  designation.  The 
vocabulary  consists  of  the  following  characters: 

01  23456789A,  = 

Any  characters  input  that  are  not  in  the  vocabulary  list  will  be 
interpreted  as  the  number  zero  (0).  Table  6  details  the  restrictions  on  the 
NACA  designation. 

NAMELIST 

This  control  card  instructs  the  program  to  print  all  namelist  data.  This 
is  useful  when  multiple  inputs  of  the  same  variable  or  namelist  are  used. 

NEXT  CASE 

This  card  indicates  termination  of  the  case  input  data  and  instructs  the 
program  to  begin  case  execution.  It  is  required  for  multiple  case  "runs".  This 
card  must  be  the  last  card  input  for  the  case. 

NOGO 

This  control  card  permits  the  program  to  cycle  through  all  of  the  input 
cases  without  computing  configuration  aerodynamics.  It  can  be  present 
anywhere  in  the  input  stream  and  only  needs  to  appear  once.  This  option  is 
useful  for  performing  error  checking  to  insure  all  cases  have  been  correctly  set 
up. 

NO  LAT 

This  control  card  inhibits  the  calculation  of  the  lateral-directional 
derivatives  due  to  sideslip  angle.  Savings  in  computation  time  can  be  realized 
by  using  this  option.  This  option  is  automatically  selected  when  using  TRIM. 

PART 

This  control  card  permits  printing  of  partial  aerodynamic  output,  such 
as  a  summary  of  the  normal  force  and  axial  force  contributors.  Partial  output 
of  the  configuration  synthesis  methods  is  only  provided  if  the  TRIM  option  is 
not  selected.  Use  of  this  card  is  equivalent  to  inputing  all  PRINT  AERO  and 
PRINT  GEOM  control  cards. 


PLOT 


A  data  file  for  use  with  a  post-processing  plotting  program  is  provided 
when  this  control  card  is  used.  A  formatted  file  is  written  to  unit  3.  Appendix 
B  shows  the  format  of  this  data  file. 

PRESSURES 


This  control  card  instructs  the  program  to  print  the  body  and  fin  alone 
pressure  coefficient  distributions  at  supersonic  speeds.  Only  pressure  data  to 
15  degrees  angle  of  attack  for  bodies  and  at  zero  angie  of  attack  for  fins  are 
printed. 

PRINT  AERO  name 


This  control  card  instructs  the  program  to  print  the  incremental 
aerodynamics  for  "name",  which  can  be  one  of  the  following: 


BODY 

for  body  aerodynamics 

FIN1 

for  FINSET1  aerodynamics 

FIN2 

for  FENSET2  aerodynamics 

FIN3 

for  FINSET3  aerodynamics 

FIN4 

for  FTNSET4  aerodynamics 

SYNTHS 

for  configuration  synthesis  aerodynamics 

TRIM 

for  trim/untrimmed  aerodynamics 

BEND 

for  panel  bending  moments 

HINGE 

for  panel  hinge  moments 

INLET 

for  inlet  aerodynamics 

All  options  are  automatically  selected  when  the  control  card  PART  is 
used.  Details  of  the  output  obtained  with  these  options  are  presented  in 
Section  4.2. 

PRINT  EXTRAP 


This  control  card  enables  the  printing  of  method  extrapolation 
messages  produced  during  execution  of  the  case.  Extrapolation  messages  are 
not  normally  provided. 

PRINT  GEOM  name 


This  control  card  instructs  the  program  to  print  the  geometric 
characteristics  of  the  configuration  component  "name",  which  can  be  one  of 
the  following: 


BODY  for  body  geometry 

FIN  1  for  FCMSET1  geometry 

FIN2  for  FINSET2  geometry 

FIN  3  for  FINSET3  geometry 

FIN 4  for  FINSET4  geometry 

INLET  for  inlet  geometry 

If  PRINT  GEOM  BODY  is  selected  and  the  Mach  number  is  greater  than 
1.2,  the  body  contour  coordinates  (X,R)  used  by  the  program  are  written  to 
tape  unit  3.  This  contour  will  contain  many  additional  points  in  between  the 
user  specified  input  coordinates. 

All  options  are  automatically  selected  when  the  control  card  PART  is 

used. 

SAVE 


The  SAVE  card  saves  namelist  inputs  from  one  case  to  the  following 
case  but  not  for  the  entire  run.  This  permits  the  user  to  build-up  or  change  a 
complex  configuration,  case-to-case,  by  adding  new  namelist  cards  without 
having  to  re-input  namelist  cards  of  the  previous  case.  When  changing  a 
namelist  that  has  been  saved,  the  namelist  must  first  be  deleted  using  the 
delete  control  card. 

The  only  control  card  that  can  be  optionally  saved,  case-to-case,  is  the 
NACA  card. 

SOSE 


The  presence  of  this  control  card  selects  the  Second-Order  Shock 
Expansion  Method  for  axisymmetric  bodies  at  supersonic  speeds.  SOSE 
should  be  selected  if  any  Mach  number  is  higher  than  2.0. 

SPIN 


When  the  SPIN  control  card  is  input,  spin  and  magnus  derivatives  are 
computed  for  body  alone.  If  the  configuration  being  run  is  a  body  +  fin  sets, 
the  spin  derivatives  are  still  computed  for  body  alone.  A  PART  or  BUILD  card 
must  be  input  for  body  alone  derivatives  to  be  printed  out. 

TRIM 


This  control  card  causes  the  program  to  perform  a  trim  calculation. 
Component  buildup  data  cannot  be  dumped  if  TTUM  is  selected.  The  use  of 
this  control  card  is  the  same  as  if  the  namelist  TRIM  was  included  except  that 
the  defaults  for  namelist  TRIM  are  used. 
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WRITE  name,  start,  end 

This  control  card  causes  the  common  block  "name"  to  be  printed  to 
tape  unit  3  using  the  most  recent  FORMAT  control  card.  Locations  from 
"start"  to  "end"  are  dumped  (see  Table  5  for  common  block  write  names).  A 
complete  definition  of  each  common  block  is  provided  in  Tables  8  through 
62.  Multiple  WRITE  statements  may  be  input,  and  there  is  no  limit  to  the 
number  which  may  be  present.  The  presence  of  a  WRITE  will  came  the  block 
"name"  to  be  printed  for  all  cases  of  the  run.  The  output  will  be  in  units  of 
feet,  pounds,  degrees,  or  degrees  Rankine.  If  the  PLOT  option  is  also  selected, 
this  output  will  be  "mixed"  with  the  PLOT  file  output  on  tape  unit  3. 

* 


Any  card  with  an  asterisk  (*)  in  Column  1  will  be  interpreted  as  a 
comment  card.  This  permits  detailed  documentation  of  case  inputs. 

3.3  TYPICAL  CASE  SET-UP 

Figure  19  schematically  shows  how  Missile  Datcom  inputs  are 
structured.  This  example  illustrates  a  multiple  case  job  in  which  case  2  uses 
part  of  the  case  1  inputs.  This  is  accomplished  through  use  of  the  SAVE 
control  card.  Case  1  is  a  body-wing-tail  configuration;  partial  output, 
component  buildup  data,  and  a  plot  file  are  created.  Case  2  uses  the  body  and 
tail  data  of  case  1  (the  wing  is  deleted  using  DELETE),  specifies  panel 
deflection  angles  and  sets  the  data  required  to  trim. 

There  is  no  limit  to  the  number  of  cases  that  can  be  "stacked"  in  a 
single  run,  provided  that  no  more  than  300  namelist  inputs  are  "saved" 
between  cases.  If  a  SAVE  control  card  is  not  present  in  a  case,  all  previous  case 
inputs  are  deleted. 

3.3.1  Configuration  Incrementing  Case  Set-up 

A  "configuration  incrementing"  case  set-up  is  shown  in  Figure  20.  This 
figure  shows  the  inputs  for  a  three  case  set-up  fin  parametric  analysis.  The 
first  case  is  the  calibration  case  with  the  remaining  cases  being  used  for  the 
parametric  analysis.  Therefore,  the  first  case  must  contain  both  the  INCRMT 
control  card  and  EXPR  namelist.  These  should  only  appear  in  the  first  case. 
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3.4  SPECIAL  USAGE  OF  INPUT  PARAMETERS 

It  is  possible  to  manipulate  the  input  geometry,  such  that  unique 
configurations  can  be  modeled.  This  section  defines  those  capabilities. 

3.4.1  Locating  Panels  on  Varying  Body  Radii  Segments 

The  fin  panels  can  be  located  anywhere  on  the  geometry.  If  they  are  to 
be  positioned  on  a  varying  radii  segment,  select  the  root  chord  span  station 
[SSPAN(l)]  such  that  the  center  of  the  exposed  root  chord  is  on  the  surface 
mold  line.  Physically  this  places  part  of  the  panel  within  the  body  and  part 
offset  from  the  body. 

If  SSPAN(l)  is  input  precisely  as  zero,  the  code  will  assume  that  panel 
semi-span  stations  relative  to  the  root  chord  are  defined.  It  will  then 
interpolate  the  body  geometry  at  the  root  chord  center  and  add  the  body 
radius  at  this  point  to  the  user  defined  values  in  the  SSPAN  array. 
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Table  1  Body  Addressable  Configurations 


CONFIGURATION 


1.  Nose  Shape 
Conical 

Sharp 
Blunted 
Truncated 
Tangent  Ogive 
Sharp 
Blunted 
Truncated 
Other 

2.  Centerbody  Shape 
Cylinder 

Elliptically  Variable 

3.  Afterbody  Shape 
Boattails 

Conical 
Tangent  Ogive 
Other 
Flares 
Conical 
Ogive 
Other 
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Table  2  Subsoric/Transonic  Method  Limitations 


METHOD 

RANGE  PERMITTED 

ACTION  TAKEN 

IF  EXCEEDED 

Nose  Bluntness  (Cn,  Cm) 

Sharp  Only 

Uses  Sharp  Method 

Conical  Nose  Slope 

0  to  25  Degrees 

Uses  25  Degrees 

Boattail  Shape 

Cone  or  Ogive 

Uses  Cone 

Conical  Boattail  Slope 

Oto  16  Degrees 

Extrapolates 

Ogive  Boattail  Slope 

0  to  28  Degrees 

Extrapolates 

Flare  Shape 

Cone 

Uses  Cone 

Flare  Slope 

Oto  10  Degrees 

Extrapolates 

Airfoil  Vc 

Oto  12% 

Continues,  If  Possible 
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Table  3  Namelist  Alphanumeric  Constants 


NAMELIST 

PERMITTED 

ALPHANUMERIC 

CONVERTED 

VALUE 

CONSTANTS 

(ALL) 

UNUSED 

1.0E-30  (Initialized  Value) 

REFQ 

TURB 

0. 

NATURAL 

1. 

AXIBOD  or  ELLBOD 

CONICAL 

0. 

CONE 

0. 

OGIVE 

1. 

POWER 

2. 

HAACK 

3. 

KARMAN 

4. 

PROTUB 

VCYL 

1. 

HCYL 

2. 

BLOCK 

5. 

FAIRING 

6. 

LUG 

3. 

SHOE 

4. 

FINSETn 

HEX 

0. 

NACA 

1. 

ARC 

2. 

USER 

3. 

INLET 

AXI 

3. 

2DTOP 

1. 

2DSIDE 

2. 

EXPR 

BODY 

1. 

FI 

2. 

F2 

3. 

F3 

4. 

F4 

5. 

BF1 

6. 

BF12 

7. 

BF123 

8. 

BF1234 

9. 
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Table  4  Equivalent  Sand  Roughness 


TYPE  OF  SURFACE 

EQUIVALENT  SAND 
ROUGHNESS 
k  (INCHES) 

RHR 

Aero  dynamically  Smooth 

0.0 

0.0 

Polished  Metal  or  Wood 

0.02  E-3  to  0.08  E-3 

6  to  26 

Natural  Sheet  Metal 

0.16  E-3 

53 

Smooth  Matte  Paint,  Carefully  Applied 

0.25  E-3 

83 

Standard  Camouflage  Paint, 

Average  Application 

0.40  E-3 

133 

Camouflage  Paint, 
Mass-Production  Spray 

1.20  E-3 

400 

Dip  Galvanized  Metal  Surface 

6.0  E-3 

2000 

Natural  Surface  of  Cast  Iron 

10.0  E-3 

3333 

PREFERRED  RHR  VALUES 


APPLICATION 

RHR 

Steel  Structural  Parts 

250 

Aluminum  and  Titanium  Structural  Parts 

125 

Close  Tolerance  Surfaces 

63 

Seals 

32 
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NAMELIST  FLTCON 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

NALPHA 

- 

NUMBER  OF  ANGLES  OF  ATTACK 
(AT  LEAST  2) 

- 

- 

ALPHA 

20 

ANGLE  OF  ATTACK  OR  TOTAL  ANGLE 

OF  ATTACK 

DEG 

“ 

BETA 

- 

SIDESLIP  ANGLE 

DEG 

0. 

PHI 

- 

AERODYNAMIC  ROLL  ANGLE 

DEG 

0. 

NMACH 

- 

NUMBER  OF  MACH  NUMBERS  (AT  LEAST  1) 

- 

- 

MACH  'I 

20 

MACH  NUMBERS 

- 

- 

REN 

20 

REYNOLDS  NUMBER  PER  UNIT  LENGTH 

IKSf© 

* 

ALT 

►  IT 

GEOMETRIC  ALTITUDE 

lie 

0. 

VINF 

f  20 

FREESTREAM  SPEED 

L/SEC(4)  - 

TINF 

20 

FREESTREAM  STATIC  TEMPERATURE 

DEG 

PINF  J 

20 

FREESTREAM  STATIC  PRESSURE 

®  THE  FOLLOWING  COMBINATIONS  SATISFY  THE  REYNOLDS  NUMBER  AND 
MACH  NUMBER  INPUT  REQUIREMENTS 


USEB  INPUT  EBQBBAftlfiflMEmES 


1 .  MACH,  REN 

2.  MACH,  ALT 

3.  VINF,  ALT 

4.  VINF,  TINF,  PINF 

5.  MACH,  TINF,  PINF 


(NONE) 

PINF,  TINF,  REN 
PINF,  TINF,  MACH,  REN 
MACH,  REN 
VINF,  REN 


®  INPUT  AS  1/FT  FOR  ENGLISH  UNITS  AND  1/M  FOR  METRIC  UNITS 
©  INPUT  AS  FT  FOR  ENGLISH  UNITS  AND  M  FOR  METRIC  UNITS 
®  INPUT  AS  FT/SEC  FOR  ENGLISH  UNITS  AND  M/SEC  FOR  METRIC  UNITS 
®  INPUT  AS  LB/FT2  FOR  ENGLISH  UNITS  AND  N/M2  FOR  METRIC  UNITS 


Figure  1  Flight  Conditions  Inputs 
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NAMELIST  REFQ 


ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

SREF 

• 

REFERENCE  AREA 

L*L 

LREF 

- 

REFERENCE  LENGTH  (LONGITUDINAL) 

L 

WMm 

LATREF 

- 

REFERENCE  LENGTH 
(LATERAL-DIRECTIONAL) 

L 

LREF 

ROUGH! /r: 
RHR  J S 

SURFACE  ROUGHNESS  HEIGHT 

L© 

0. 

ROUGHNESS  HEIGHT  RATING 

- 

0. 

XCG 

- 

LONGITUDINAL  POSITION  OF  C.G.  (+AFT) 

L 

0. 

ZCG 

* 

VERTICAL  POSfTION  OF  C.G.  (+UP) 

L 

0. 

SCALE 

- 

VEHICLE  SCALE  FACTOR  (MULTIPLIER  TO 
GEOMETRY) 

“ 

1. 

BLAYER 

BOUNDARY  LAYER  TYPE:  TURB  FOR  FULLY 
TURBULENT.  NATURAL  FOR  NATURAL 
TRANSITION 

TURB 

®  EITHER  CAN  BE  USED 

®  DEFAULT  IS  BODY  MAXIMUM  CROSS-SECTIONAL  AREA.  IF  NO  BODY  IS  INPUT. 
MAXIMUM  FIN  PANEL  AREA  IS  USED. 

®  DEFAULT  IS  BODY  MAXIMUM  DIAMETER.  IF  NO  BODY  IS  INPUT.  MAXIMUM  FIN  PANEL 
MEAN  GEOMETRIC  CHORD  IS  USED. 

©  INPUT  AS  INCHES  FOR  ENGLISH  UNITS  AND  CENTIMETERS  FOR  METRIC  UNITS. 


Figure  2  Reference  Quantity  Inputs 
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NAMELIST  AXIBOD 

OPTION  1  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XO  OR  XO 

- 

LONGITUDINAL  COORDINATE  AT 

L 

0. 

NOSE  TIP 

TNOSE 

- 

NOSE  SHAPE  NAME:  CONICAL.  CONE, 

- 

OGIVE 

OGIVE,  POWER,  HAACK,  OR  KARMAN 

POWER 

- 

EXPONENT,  n,  FOR  POWER  SERIES 

- 

0. 

SHAPES,  (r/R)  =  (x/L)" 

LNOSE 

- 

NOSE  LENGTH 

L 

- 

DNOSE 

- 

NOSE  SECTION  BASE  DIAMETER 

L 

1.0 

BNOSE 

- 

NOSE  BLUNTNESS  RADIUS  OR 

L 

0. 

TRUNCATED  NOSE  RADIUS 

TRUNC 

- 

TRUNCATION  FLAG  (.TRUE.  IF  NOSE  IS 

- 

.FALSE. 

TRUNCATED) 

LCENTR 

- 

CENTERBODY  LENGTH 

L 

0. 

DCENTR 

- 

CENTERBODY  BASE  DIAMETER 

L 

DNOSE 

TAFT 

- 

AFTERBODY  SHAPE  NAME:  CONICAL. 

- 

CONICAL 

J 

T) 

CONE.  OR  OGIVE 

LAFT  _ 

\ 

1/ 

AFTERBODY  LENGTH 

L 

0. 

daftCDJ 

- 

AFTERBODY  BASE  DIAMETER 

L 

- 

DEXIT 

- 

DIAMETER  OF  NOZZLE  EXIT 

L 

- 

BASE 

- 

BASE-JET  PLUME  INTERACTION  FLAG 

- 

.FALSE. 

(.TRUE.  IF  CALCULATIONS  DESIRED) 

BETAN 

> 

NOZZLE  EXIT  ANGLE 

DEG 

- 

JMACH 

|  20 

JET  MACH  NUMBER  AT  NOZZLE  EXIT 

- 

- 

PRAT 

►! 

20 

JET  TO  FREESTREAM  STATIC  PRESSURE 

- 

- 

n 

RATIO 

TRAT  J 

1 

20 

JET  TO  FREESTREAM  STAGNATION 

- 

- 

_l 

TEMPERATURE  RATIO 

®  AFT  BODY  MUST  NOT  BE  CYLINDRICAL  (i.e.  DAFT  NOT  EQUAL  TO  DCENTR) 

<2>  ONLY  REQUIRED  IF  BASE-JET  PLUME  INTERACTION  CALCULATIONS  DESIRED. 
(DEXIT  MUST  NOT  EQUAL  ZERO) 

<D  JMACH,  PRAT.  ANDTRAT  ARE  SPECIFIED  FOR  EACH  FREESTREAM  MACH 
NUMBER  OR  VELOCITY  INPUT  IN  NAMELIST  $FLTCON 
®  DAFT  MUST  BE  NON-ZERO 


Figure  3a  Axisymmetric  Body  Geometry  Inputs  -  Option  1 
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NAMELIST  AXIBOD 

OPTION  2  INPUTS  (USE  ONLY  IF  MACH  >  1.2) 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XOORXO 

- 

LONGITUDINAL  COORDINATE  AT  NOSE  TIP 

0. 

BNOSE 

• 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED 
NOSE  RADIUS 

■ 

0. 

TRUNC 

- 

TRUNCATION  FLAG  (.TRUE.  IF  NOSE  IS 
TRUNCATED) 

■ 

.FALSE. 

DEXIT 

- 

DIAMETER  OF  NOZZLE  EXIT 

- 

NX 

- 

NUMBER  OF  INPUT  STATIONS 

■ 

- 

x  © 

50 

LONGITUDINAL  COORDINATES 

■  jH 

- 

R 

50 

RADIUS  AT  EACH  X  STATION 

1 

• 

DISCON 

20 

INOICES  OF  X  STATIONS  WHERE  THE 
SURFACE  SLOPE  IS  DISCONTINUOUS 

■ 

®  X(NX)  MUST  BE  END  OF  BODY 


Figure  3b  Axisymmetric  Body  Geometry  Inputs  -  Option  2 
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BLUNTED  NOSE  TRUNCATED  NOSE 

BNOSE 


Figure  4  Body  Geometry  Inputs 
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Input  Parameter 

Symbol 

Min. 

Value 

Max. 

Value 

Boa  nail  shape 

-- 

Cylinder.  O 

ine.  Ogive 

Boattajl  fineness  ratio 

UD 

0 

2 

Boattail  terminal  angle 

Pt 

C' 

12° 

Jet  pressure  ratio 

Pj'P- 

0 

10 

Freestream  Mach 
number 

Moo 

2 

5 

Angle  of  Attach 

a 

0° 

8° 

Jet  Mach  number 

Mj 

Moo  -1 

M„  +  1 

Nozzle  terminal  angle 

Pn 

5° 

25° 

Jet  diameter  ratio 

Qj/Db 

0  00 

095 

Jet  temperature  ratio 

Tt  j  /Tt=- 

4 

10 

Not*;  If  Input  parameter  Is  not  between  minimum  and 
maximum  value  the  code  will  extrapolate 


Figure  5  Base-Jet  Plume  Interaction  Parameters 
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NAMELIST  ELLBOD 

OPTION  1  INPUTS 


1 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XO  OR  XO 

- 

LONGITUDINAL  COORDINATE  AT 

NOSE  TIP 

L 

0. 

TNOSE 

“ 

NOSE  SHAPE  NAME:  CONICAL.  CONE. 
OGIVE,  POWER,  HAACK,  OR  KARMAN 

“ 

OGIVE 

POWER 

* 

EXPONENT,  n.  FOR  POWER  SERIES 
SHAPES,  (r/RMx/L)" 

* 

0. 

1  LNOSE 

- 

NOSE  LENGTH 

L 

- 

WNOSE 

- 

NOSE  SECTION  BASE  WIDTH 

L 

1.0 

BNOSE 

“ 

NOSE  BLUNTNESS  RADIUS  OR 

TRUNCATED  NOSE  RADIUS 

L 

0. 

TRUNC 

“ 

TRUNCATION  FUG  (.TRUE.  IF  NOSE  IS 
TRUNCATED) 

- 

FALSE. 

ENOSE 

ELLIPTICITY  AT  NOSE  BASE  (H/W) 

- 

1.0 

LCENTR 

CENTERBODY  LENGTH 

L 

0. 

WCENTR 

CENTERBODY  BASE  WIDTH 

L 

WNOSE 

ECENTR 

ELLIPTICITY  AT  CENTERBODY  BASE  (H/W) 

- 

1.0 

TAFT 

AFTERBODY  SHAPE  NAME:  CONICAL. 
CONE.  OR  OGIVE 

CONICAL 

LAFT 

- 

AFTERBODY  LENGTH 

L 

0. 

WAFT 

“ 

AFTERBODY  BASE  WIDTH 

L  j 

- 

EAFT  J 

- 

ELLIPTICITY  AT  AFT  BODY  BASE  (H/W) 

- 

1.0 

DEXIT 

_ 

DIAMETER  OF  NOZZLE  EXIT 

L 

- 

®  AFT  BODY  MUST  NOT  BE  CYLINDRICAL  (i.e.  WAFT  NOT  EQUAL  TO  WCENTfl) 


Figure  6a  Elliptical  Body  Geometry  Inputs  -  Option  1 
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NAMELIST  ELLBOD 

OPTION  2  INPUTS  (USE  ONLY  IF  MACH  >  1.2) 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XO  OR  XO 

- 

LONGITUDINAL  COORDINATE  AT  NOSE  TIP 

L 

0. 

BNOSE 

• 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED 
NOSE  RADIUS 

L 

0. 

TRUNC 

TRUNCATION  FLAG  (.TRUE.  IF  NOSE  IS 
TRUNCATED) 

" 

FALSE. 

DEXIT 

- 

DIAMETER  OF  NOZZLE  EXIT 

L 

- 

NX 

- 

NUMBER  OF  INPUT  STATIONS 

- 

- 

X  O 

50 

LONGITUDINAL  COORDINATES 

L 

- 

W  ( 2) 

50 

BODY  HALF-WIDTH  AT  EACH  X  STATION 

L 

- 

DISCON 

20 

INDICES  OF  X  STATIONS  WHERE  THE 
SURFACE  SLOPE  IS  DISCONTINUOUS 

- 

- 

H<I> 

50 

BODY  HALF-HEIGHT  AT  EACH  X  STATION 

L 

. 

ELLIP  (5) 

50 

BODY  HEIGHT  TO  WIDTH  RATIO  AT  EACH 
!  X  STATION 

* 

1.0 

®  X(NX)  MUST  BE  END  OF  BODY 

@  ONE  OF  7>!E  FOLLOWING  COMBINATIONS  IS  REQUIRED: 
W  AND  H.  W  AND  ELLIP,  OR  H  AND  ELLIP 


Figure  6b  Elliptical  Body  Geometry  Inputs  -  Option  2 
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NAMELIST  PROTUB 


. 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

NPROT 

- 

NUMBER  OF  PROTUBERANCE  SETS 
(20  MAXIMUM) 

- 

0. 

PTYPE 

20 

PROTUBERANCE  SET  TYPE:  VCYL,  HCYL, 
BLOCK,  FAIRING.  LUG,  OR  SHOE  © 

~ 

* 

XPROT 

20 

LONGITUDINAL  DISTANCE  FROM  MISSILE 
NOSE  TO  PROTUBERANCE  SET 

L 

- 

NLOC  CD 

20 

NUMBER  OF  PROTUBERANCES  IN  EACH 
PROTUBERANCE  SET 

“ 

0. 

LPROT 

100 

LENGTH  OF  EACH  MEMBER  OR 
PROTUBERANCE 

L 

• 

WPROT 

4 

L  100 

P 

WIDTH  OF  EACH  MEMBER  OR 
PROTUBERANCE 

L 

- 

HPROT 

100 

HEIGHT  OF  EACH  MEMBER  OR 
PROTUBERANCE 

L 

• 

OPROT  J 

100 

VERTICAL  OFFSET  OF  EACH  MEMBER 

OR  PROTUBERANCE 

L 

0. 

®  NLOC  ACCOUNTS  FOR  IDENTICAL  PROTUBERANCES  (SAME  SIZE  AND  SHAPE)  LOCATED 
AROUND  THE  BODY  AT  THE  SAME  AXIAL  LOCATION. 

©  LUG  TYPE  HAS  4  MEMBERS.  SHOE  TYPE  HAS  3  MEMBERS.  LPROT,  WPROT,  HPROT,  AND 
OPROT  MUST  BE  SPECIFIED  FOR  EACH  MEMBER. 

®  INPUT  FOR  EACH  PROTUBERANCE  (VCYL.  HCYL,  BLOCK,  OR  FAIRING  TYPE)  OR  EACH 
PROTUBERANCE  MEMBER  (LUG  AND  SHOE  TYPE) 


Figure  7  Protuberance  Inputs 


V 


TOP  VIEW 


3-D  VIEW 


H 


Vertical  Circular  Cylinder 
Oriented  Perpendicular  To 
Flow  (VCYL): 


Horizontal  Circular 
Cylinder  Oriented  In  Flow 
Direction  (HCYL): 


Vertical  Elliptical  Cylinder 
Oriented  Perpendicular  To 
Flow  With  L>  W  Or  W>L 
(VCYL): 


Block  (BLOCK): 


Streamline  Half  Bodies 
(FAIRING): 


Launch  Shoe  (SHOE): 

1)  Base  Member 

2)  Vertical  Member 

3)  Top  Member 


Launch  Lug  (LUG): 

1)  Base  Member 

2)  &  3)  Vertical  Members 

4)  Top  Member 


L-Length  W-Wldth  H-Helght 


Note;  Length,  width,  height,  and  offset  must  be  input  for  each  member  of  launch  lug  and 
launch  shoe  types 

Note;  Offset  Is  the  perpendicular  distance  from  the  missile  mold  line  to  the  bottom  of  the 
protuberance  or  protuberance  member 

Figure  8  Protuberance  Shapes  Available 
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CASE ID  PROTUBERANCE  EXAMPLE  CASE 
DIM  IN 
NO  LAT 

$REFQ  XCG=39 . 0,  $ 

$FLTCON  NMACH=3 . , MACH=0 . 4, 0 . 8, 2 . 0, 

REN=3 . E06, 3 . E06, 3 . E06, ALT=0 . 0 , 

NALPHA=5 . , ALPHA=-8 .  ,  -4 .  ,  0 . , 4 . , 8 . , $ 

$ AX I BOD  TNOSE=OGIVE , LNOSE=12 . 0 , DNOSE= 12 . 0 , 
LCENTR=54 . 0 , DCENTR=12 . 0, 

TAFT=CONE , LAFT=12 . 0 , DAFT =6 . 0 , DEXIT=5 . 0 , $ 
$ PRO TUB  NPROT=4 . , 

P TYPE— FAIRING , VCYL, SHOE, BLOCK, 

XPROT=14 . , 22 . , 39 . , 56 . , 

NLOC=2 . , 4 . , 2 . , 1 . , 

LPROT=5 . , 1 . , 10 . , 10 . , 10 . ,  0 . 5 , 

WPROT=2., 1., 4., 0.25,1.,!., 

HPROT=2 .,0.5,0.1,0.75,0.25,0.25, 

OPROT=0 ., 0. , 0. , 0.1,0. 85, 0.,$ 

$F INSET 1  SSPAN=0 . 0, 9 . 0, CHORD=14 . 0, 8 . 0, 

XLE=64 . 0, SWEEP=0 . 0, STA=1 . 0,NPANEL=4 . , 
PHIF=45 . , 135 . , 225 . , 315 . , $ 

PRINT  GEOM  BODY 
PRINT  AERO  BODY 
SAVE 

NEXT  CASE 


NOTE;  Length,  Width,  and  Height  is  input  for  each  member  of  the  launch  shoe 

Figure  9  Protuberance  Example  Input  File 
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NAMELIST  FINSETn 

NOMINAL  INPUTS 


VARIABLE 

NAME 

DEFINITION 

UNITS 

DEFAULT 

SECTYP 

. 

TYPE  OF  SECTION  TO  BE  DEFINED:  HEX. 

- 

HEX 

SSPAN© 

NACA,  ARC,  OR  USER 

10 

SEMI-SPAN 

L 

- 

CHORD 

10 

PANEL  CHORD  LENGTH  AT  EACH  SSPAN 

L 

- 

XLE 

10 

DISTANCE  FROM  NOSE  TIP  TO  CHORD 

L 

SWEEP  (5) 

LEADING  EDGE  AT  EACH  SSPAN 

10 

SWEEPBACK  ANGLE  AT  EACH  SSPAN 

DEG 

STA 

10 

CHORD  STATION  USED  IN  MEASURING 

- 

1.0 

LER® 

SWEEP  AT  EACH  SSPAN  (0.0=LEADING 
EDGE.  1 .0*TR AILING  EDGE) 

10 

PANEL  LEADING  EDGE  RADIUS  AT  EACH 

L 

SPAN  STATION 

NPANEL 

- 

NUMBER  OF  PANELS  IN  SET  (1-8) 

- 

4 

PHIF 

8 

ROLL  ANGLE  OF  EACH  FIN  MEASURED 

DEG 

© 

CLOCKWISE  FROM  TOP  VERTICAL 
CENTER  LOOKING  FORWARD 

GAM 

8 

DIHEDRAL  OF  EACH  FIN  (POSITIVE  WHEN 

DEG 

PHIF  IS  INCREASED.  SEE  FIG.  12) 

SKEW 

- 

ANGLE  BETWEEN  THE  Y  AXIS  AND  THE 

DEG 

FIN  HINGE  LINE  (POSITIVE  SWEPT  BACK) 

CFOC 

10 

FLAP  CHORD  TO  FIN  CHORD  RATIO 

* 

1.0 

(D  IF  SSPAN(1)=0.0,  INPUTS  ARE  RELATIVE  TO  ROOT  CHORD  NOT  BODY  CENTERLINE 
@  IF  USING  SWEEP.  SPECIFY  ONLY  XLE(1);  IF  USING  XLE  DO  NOT  SPECIFY  SWEEP 
<2>  NOT  REQUIRED  FOR  NACA  AIRFOILS,  REQUIRED  FOR  USER  AIRFOILS 
0  IF  PHIF  NOT  INPUT  THE  NUMBER  OF  PANELS  ARE  EVENLY  SPACED  ABOUT  THE  BODY. 


Figure  10a  Fin  Geometry  Inputs  -  Nominal 
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I 


NAMELIST  FINSETn 

OPTIONAL  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

ZUPPER 

10 

MAXIMUM  THICKNESS  TO  CHORD  RATIO 

OF  UPPER  SURFACE 

- 

0.025 

ZLOWER 

10 

MAXIMUM  THICKNESS  TO  CHORD  RATIO 

OF  LOWER  SURFACE 

- 

ZUPPER 

LMAXU 

10 

FRACTION  OF  CHORD  FROM  SECTION 
LEADING  EDGE  TO  MAXIMUM 
THICKNESS  OF  UPPER  SURFACE 

0.5 

LMAXL 

10 

FRACTION  OF  CHORD  FROM  SECTION 
LEADING  EDGE  TO  MAXIMUM 
THICKNESS  OF  LOWER  SURFACE 

LMAXU 

LFLATU 

10 

FRACTION  OF  CHORD  OF  CONSTANT 
THICKNESS  SECTION  ON  UPPER 
SURFACE 

0.0 

LFLATL 

10 

FRACTION  OF  CHORD  OF  CONSTANT 
THICKNESS  SECTION  ON  LOWER 
SURFACE 

LFLATU 

NOTE;  THESE  PARAMETERS  MUST  BE  INPUT  AT  EACH  SPAN  STATION 


Figure  10b  Fin  Geometry  Inputs  -  Optional 
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NAMELIST  FINSETn 

INPUTS  FOR  "USER”  SECTIONS 


VARIABLE 

NAME 

array 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XCORD 

50 

CHORD  STATION,  FRACTION  OF  CHORD 

- 

. 

MEAN  0 

FROM  LEADING  EDGE 

50 

DISTANCE  BETWEEN  THE  MEAN  LINE  AND 

- 

- 

CHORD  LINE  AT  EACH  XCORD. 

FRACTION  OF  CHORD 

THICK  © 

50 

THICKNESS  TO  CHORD  RATIO  AT  EACH 

- 

- 

YUPPEF0 

XCORD 

1  _  50 

UPPER  SURFACE  COORDINATES. 

- 

- 

FRACTION  OF  CHORD.  AT  EACH  XCORD 

ylowefC? 

J  50 

LOWER  SURFACE  COORDINATES. 

- 

- 

FRACTION  OF  CHORD.  AT  EACH  XCORD 

NOTE;  ALL  VARIABLES  ARE  EXPRESSED  AS  FRACTIONS  OF  CHORD 
LEADING  EDGE  RADIUS  (VARIABLE  LER)  MUST  BE  DEFINED 

CD  EITHER  MEAN  AND  THICK  OR  YUPPER  AND  YLOWER  ARE  REQUIRED 
@  THE  AIRFOIL  SECTION  MUST  BE  CLOSED.  THIS  MEANS  THE  FINAL  CHORD 
STATION  INPUT  MUST  HAVE  EITHER  THICK=0.0  OR  YUPPER* YLOWER. 


+Y 


Figure  10c  Fin  Geometry  Inputs  -  User  Airfoils 
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NOTE;  XLE  IS  MEASURED  FROM  NOSE  TIP 

IF  SSPAN(1)  IS  INPUT  AS  ZERO,  SSPAN  INPUTS 
ARE  RELATIVE  TO  BODY  SURFACE  MOLD  LINE 

Figure  11a  Selecting  Panel  Break  Points 
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S$?AN 


®span  .VoV£  ,25’  jo, 

:FOC  -OOo1225’30 

\  0  '°  30,  O  30  0.0, 


SSPAN  -oS'°  °' 

FOC  **  0.33,  |  0 

V^e.Theflapc. 

\  ‘,>«^t?0a,me 

\  »»  fin  taper  r^n  *he" 
\  ^Jer  fa«>o  «  ?ef0 


Valld  Trailing 


Edaa  Device 


D*finltic 


wm  l  / 


•  Trailing 

.  3  0  Tailing  Edge  Devirec 


F  INSET  2.  PANEL  3 


r 

TOP  VERTICAL 
CENTER 


PANEL  2 


VIEW  LOOKING  FORWARD 
FROM  BEHIND  MISSILE 


Figure  12  Fin  Numbering  and  Orientation 


VIEW  LOOKING  FORWARD 
FROM  BEHIND  MISSILE 


PHI  IS  THE  BODY  ROLL  ANGLE 
PHIF  IS  THE  FIN  PANEL  ROLL  ANGLE 


Figure  13  Roll  Attitude  vs  Fin  Orientation 


NAMELIST  DEFLCT 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

DELTA1 0 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  1  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

DELTA20 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  2  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

DELTA30 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  3  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

DELTA40 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  4  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

XHINGE 

4 

DISTANCE  FROM  COORDINANT  SYSTEM 
ORIGIN  TO  PANEL  HINGE  LINE  FOR  EACF 
FIN  SET 

L 

XO+XLE 

+CR/2© 

SKEW 

4 

SWEEPBACK  OF  HINGE  LINE  FOR  EACH 

FIN  SET 

DEG 

0. 

®  PANEL  NUMBERING  IS  SHOWN  IN  FIGURE  12 
®  DEFAULT  IS  AT  ONE-HALF  THE  EXPOSED  ROOT  CHORD. 
AS  MEASURED  FROM  THE  COORDINANT  SYSTEM  ORIGIN. 


NOTE;  A  POSITIVE  DEFLECTION  ANGLE  PRODUCES  A  NEGATIVE  BODY  AXIS 
ROLLING  MOMENT  AT  ZERO  ANGLE  OF  ATTACK 


Figure  14  Panel  Deflection  Inputs 
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NAMELIST  TRIM 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

Mi I 

■ 

- 

FIN  SET  TO  BE  USED  FOR  TRIM 

- 

1. 

ESSi 

m 

- 

.TRUE.  IF  PANEL  TO  BE  USED 

FALSE. 

■ 

- 

.TRUE.  IF  PANEL  TO  BE  USED 

FALSE. 

PANL3 

- 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL4 

- 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL5 

* 

.TRUE.  IF  PANEL  TO  BE  USED 

FALSE. 

PANL6 

- 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL7 

- 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL8. 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

delminI^J 

V 

MINIMUM  NEGATIVE  DEFLECTION 

DEG 

-25. 

DELMAXr^ 

) 

MAXIMUM  POSITIVE  DEFLECTION 

DEG 

+20. 

ASYM 

* 

8 

.TRUE.  IF  PANEL  WITH  SUBSCRIPT  IS  TO 

BE  DEFLECTED  OPPOSITE  TO  NORMAL 
SIGN  CONVENTION  (ASYMMETRIC 
DEFLECTIONS) 

.FALSE. 

®  DEFAULTS  APPLY  ONLY  IF  ALL  PANL#  DATA  ARE  NOT  INPUT  OR  FALSE. 
®  BOTH  DELM1N  AND  DELMAX  MUST  BE  SPECIFIED 


Figure  15  Trim  Inputs 
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NAMELIST  INLET 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

Inin 

- 

NUMBER  OF  INLETS  (MAXIMUM=20) 

- 

- 

I  intype 

- 

TYPE  OF  INLET:  2DTOP,  2DSIDE,  OR  AXI© 

- 

- 

XINLT 

LONGITUDINAL  DISTANCE  FROM  NOSE  TIP 
TO  INLET  LEADING  EDGE 

L 

- 

XDIV 

LONGITUDINAL  DISTANCE  FROM  INLET 
LEADING  EDGE  TO  DIVERTER 

LEADING  EDGE 

L 

HD  IV 

- 

HEIGHT  OF  DIVERTER  LEADING  EDGE 

L 

- 

LDIV 

- 

LENGTH  OF  DIVERTER 

L 

- 

PHI© 

20 

INLET  ROLL  ORIENTATIONS 

DEG 

- 

X  ' 

5 

INLET  LONGITUDINAL  POSITIONS  RELATIVE 
TO  INLET  LEADING  EDGE 

L 

- 

H© 

►<3) 

5 

INLET  HEIGHTS  AT  THE  LONGITUDINAL 
POSITIONS 

L 

w(dJ 

5 

INLET  WIDTHS  AT  THE  LONGITUDINAL 
POSITIONS 

L 

COVER 

- 

IF  COVER=.TRUE.  INLETS  ARE  COVERED 

- 

FALSE. 

(ramp 

- 

EXTERNAL  COMPRESSION  RAMP  ANGLE 

DEG 

- 

ADD 

IF  ADD=.TRUE.  INLET  ADDITIVE  DRAG  IS 
CALCULATED 

• 

FALSE. 

MFR  © 

20 

MASS  FLOW  RATIO  FOR  EACH  MACH 
NUMBER  (0<  MFR  <  1.0) 

' 

(D  2DTOP:  TWO  DIMENSIONAL  TOP  MOUNTED,  2DSIDE:  TWO  DIMENSIONAL  SIDE  MOUNTED, 
AXI:  AXISYMMETRIC 

<2>  ROLL  POSITIONS  FROM  TCP  VERTICAL  CENTER.  SAME  CONVENTION  AS  FIN  ROLL 
POSITIONS. 

CD  SPECIFY  X,  H.  AND  W  AT  FIVE  INLET  LOCATIONS:  1)  LEADING  EDGE.  2)  COWL  LIP  LEADING 
EDGE,  3)  MIDBODY  START,  4)  START  OF  BOATTAIL,  5)  END  OF  BOATTAIL 
®  NOT  REQUIRED  IF  INTYPE=AXI 
®  IF  INTYPE=AXI,  W=DIAMETER 

®  SPECIFY  MASS  FLOW  RATIO,  MFR,  FOR  EACH  FREE  STREAM  MACH  NUMBER  OR  VELOCITY 
GIVEN  IN  $FLTCON  (ONLY  REQUIRED  IF  ADD  =  TRUE.) 


Figure  16  Inlet  Geometry  Inputs 


SIDE  VIEW 


INPUT  INLET  LONGITUDINAL  COORDINATE  (X),  HEIGHT  (H). 
AND  WIDTH  (W)  AT  EACH  STATION 

STATION  INLET  POSITION 

(1)  LEADING  EDGE  (X(1)  A  H(1)-0) 

(2)  COWL  LIP  LEADING  EDGE 

{3)  MIDBODY  START 

(4)  BOATTAIL  START 

(5)  BOATTAIL  END 


NOTES. 

•  INLET  ROLL  ORIENTATION  IS  SAME  CONVENTION  AS  FIN  ROLL  ORIENTATION. 

•  RAMP  IS  THE  EXTERNAL  COMPRESSION  RAMP  ANGLE  (SHOWN  SHADED  IN  THE  SIDE  VIEW) 

•  HEIGHT  OF  THE  DIVERTER  IS  SPECIFIED  AT  THE  DIVERTER  LEADING  EDGE 

•  THE  DIVERTER  WIDTH  IS  EQUAL  TO  THE  INLET  WIDTH  AT  LDIV 

•  IF  INLET  IS  COVERED  (COVER- TRUE )  A  PLUG  IS  PLACED  BETWEEN  STATIONS  1  AND  2 
FLUSH  WITH  THE  INLET  FACE 


Figure  17a  Top-Mounted  2-D  Inlet /Diverter  Geometry 
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TOP  VIEW 


H(1)=0 


H(1)-0 


TOP  VERTICAL 
CENTER 


INPUT  INLET  LONGITUDINAL  COORDINATE  (X),  HEIGHT  (H), 
AND  WIDTH  (W)  AT  EACH  STATION 

STATION  INLET  POSITION 

(1)  LEADING  EDGE  (X(t)  &  H(1)  -  0) 

(2)  COWL  LIP  LEADING  EDGE 

(3)  MIDBODY  START 

(4)  BOATTAIL  START 

(5)  BOATTAIL  END 


ibo- 


INLET  ALONE 


NOTES, 

•  INLET  ROLL  ORIENTATION  IS  SAME  CONVENTION  AS  FIN  ROLL  ORIENTATION. 

•  RAMP  IS  THE  EXTERNAL  COMPRESSION  RAMP  ANGLE  (SHOWN  SHADED  IN  THE  SIDE  VIEW) 

•  HEIGHT  OF  THE  DIVERTER  IS  SPECIFIED  AT  THE  DIVERTER  LEADING  EDGE 

•  THE  DIVERTER  WIDTH  IS  EQUAL  TO  THE  INLET  WIDTH  AT  LDIV 

•  IF  INLET  IS  COVERED  (COVER-  TRUE  )  A  PLUG  IS  PLACED  BETWEEN  STATIONS  1  AND  2 
FLUSH  WITH  THE  INLET  FACE 


Figure  17b  Side-Mounted  2-D  Inlet/Diverter  Geometry 
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(4) 


(5) 


X(1)-0 

W(1)=0 


TOP  VERTICAL 
CENTER 

♦ 


180- 


INPUT  INLET  LONGITUDINAL  COORDINATE  (X),  AND 
DIAMETER  (W)  AT  EACH  STATION 

STATION  INLET  POSITION 

(1)  TIP  (X(1)  A  W(1)  -  0) 

(2)  COWL  LIP  LEADING  EDGE 

(3)  MIDBODY  START 

(4)  BOATTAIL  START 

(5)  BOATTAIL  END 


INLET  ALONE 


NOTES; 

•  INLET  ROLL  ORIENTATION  IS  SAME  CONVENTION  AS  FIN  ROLL  ORIENTATION. 

•  RAMP  IS  THE  EXTERNAL  COMPRESSION  CONE  HALF-ANGLE  (SHOWN  SHADED 
IN  THE  SIDE  VIEW) 

•  HEIGHT  OF  THE  DIVERTER  IS  SPECIFIED  AT  THE  DIVERTER  LEADING  EDGE 

•  THE  DIVERTER  WIDTH  IS  EQUAL  TO  THE  INLET  DIAMETER  AT  LDIV 

•  IF  INLET  IS  COVERED  (COVER-TRUE.)  A  PLUG  IS  PLACED  BETWEEN  STATIONS  1  AND  2 
FLUSH  WITH  THE  INLET  FACE 


Figure  17c  Axisymmetric  Inlet/ Diverter  Geometry 
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HDIV=0 


Conformal  Inlets 


Semi-Submerged  Inlets 


Figure  17d  Geometry  Definition  For  Conformal  And  Semi-Submerged  Inlets 
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NAMELIST  EXPR 


DEFINITION 

UNITS 

DEFAULT 

MACH 

- 

MACH  NUMBER 

- 

- 

NALPHA 

- 

NUMBER  OF  ANGLES  OF  ATTACK  (2-20) 

- 

- 

ALPHA 

20 

ANGLES  OF  ATTACK  FOR  DATA 

DEG 

- 

SREF 

- 

REFERENCE  AREA  FOR  DATA 

L*L 

© 

LREF 

• 

LONGITUDINAL  REFERENCE  LENGTH 

FOR  DATA 

L 

© 

LATREF 

- 

LATERAL  REFERENCE  LENGTH  FOR  DATA 

L 

LREF 

XCG 

- 

LONGITUDINAL  C.G.  FOR  DATA 

L 

0. 

ZCG 

- 

VERTICAL  C.G.  FOR  DATA 

L 

0. 

CONF 

CONFIGURATION  FOR  DATA 

SELECT  ONE  OF  THE  FOLLOWING 

BODY  -  BODY 

FI  -  WING 

F2  -  TAIL 

F3  -  THIRD  FIN  SET 

F4  -  FOURTH  FIN  SET 

BF1  -  BODY-WING 

BF12  -  BODY-2  FIN  SETS 

BF123  -  BOOY-3  FIN  SETS 

BF1234  -  BODY-4  FIN  SETS 

CN 

20 

Cn  DATA  VS  ALPHA 

- 

CM 

20 

Cm  DATA  VS  ALPHA 

- 

CA 

20 

Ca  DATA  VS  ALPHA 

- 

CY 

20 

Cy  DATA  VS  ALPHA 

- 

CSN 

20 

Cn  DATA  VS  ALPHA 

- 

- 

CSL 

20 

Ci  DATA  VS  ALPHA 

- 

®  DEFAULT  IS  BODY  MAXIMUM  CROSS-SECTIONAL  AREA.  IF  NO  BODY  IS  INPUT. 
MAXIMUM  FIN  PANEL  AREA  IS  USED. 

®  DEFAULT  IS  BODY  MAXIMUM  DIAMETER.  IF  NO  BODY  IS  INPUT.  MAXIMUM  FIN  PANEL 
MEAN  GEOMETRIC  CHORD  IS  USED. 


Figure  18  Experimental  Data  Inputs 
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Table  5  Common  Block  DUMP  and  WRITE  Names 


Table  6  Airfoil  Designation  Using  the  NACA  Control  Card 


INPUT  NACA 
DESIGNATION 

NACA  SERIES 
AIRFOIL 

RESTRICTIONS 

0012.25 

4-Digit 

None.  Fractional  thickness  may  be 
specified. 

23118.50 

5-Digit 

None.  Fractional  thickness  may  be 
specified. 

2406-32 

4-Digit 

modified 

Sixth  digit  specifies  position  of 
maximum  thickness,  (%chord/10),  and 
must  be  a  2, 3, 4, 5,  or  6. 

43006-65 

5-Digit 

modified 

Seventh  digit  specifies  position  of 
maximum  thickness,  (%chord/10),  and 
must  be  a  2, 3, 4, 5,  or  6. 

16-212.25 

1-Series 

Second  digit  specifies  location  of 
minimum  pressure,  (%chord/10),  and 
must  be  a  6, 8,  or  9.  Fractional 
thickness  may  be  specified. 

64-005 

64-205  A=0.6 
63A005 

652A215  A=0.8 

65,2 A21 5  A=0.8 

6-Series 

Second  digit  specifies  location  of 
minimum  pressure,  (%chord/10),  and 
must  be  a  3, 4, 5,  or  6.  The  mean  line 
parameter  (A=xx)  must  be  a  decimal 
between  0.1  and  1.0  (Default  isl.0). 

Se®  Note  1. 

3-30.0-2.5-40.1 

A  B  C  D 

Supersonic 

See  Note  2. 

A  -  Section  type:  1= Double  Wedge 
2=Circular  Arc 
3=Hexagonal 

B  -  Distance  from  leading  edge  to 
position  of  maximum  thickness,  %  of 
chord. 

C  -  Maximum  thickness,  %  of  chord. 

D  -  For  hexagonal  sections,  length  of 
surface  of  constant  thickness,  %  of 
chord. 

Note  1.  The  program  does  not  distinguish  between  a  64,2-220  and  a  64-220 
specification.  The  difference  in  coordinates  between  the  two  is  negligible. 


Note  2.  All  parameters  can  be  expressed  to  0.1%.  The  delimiter must  be  used. 


66 


Figure  19  Typical  "Stacked"  Case  Set-up 
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Figure  20  "Configuration  Incrementing"  Case  Set-up 
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4.0  OUTPUT  DESCRIPTION 

This  section  describes  the  types  of  output  available  from  the  code.  In 
most  cases  the  available  output  is  user  selectable,  that  is,  it  is  not  normally 
provided  and  must  be  specifically  requested  using  a  specialized  control  card. 
This  feature  permits  the  user  to  tailor  the  code  output  to  fit  his  particular 
application  without  extensive  reprogramming.  This  allows  him  to  find  the 
output  that  he  is  interested  in  without  having  to  wade  through  output  that 
does  not  interest  him. 

The  following  four  types  of  output  are  available  from  the  code: 

•  Nominal  output  -  This  output  is  always  provided  by  the  code 
and  consists  of  output  from  the  input  error  checking  module 
(CONERR),  a  listing  of  the  inputs  for  each  case,  and  the  final 
aerodynamic  results  for  the  configuration. 

•  Partial  output  -  This  output  details  the  configuration 
geometry  and  the  intermediate  aerodynamic  calculations. 

Special  control  cards  are  available  so  that  the  user  can  select 
the  quantity  and  types  of  output  desired 

•  External  data  files  -  This  output  permits  the  user  to  create 
external  data  files  which  can  be  used  in  post-processing 
programs,  such  as  plotting  or  trajectory  programs.  Both  fixed 
and  user  defined  format  data  files  can  be  created  with  the 
addition  of  simple  control  cards. 

•  Array  dumps  and  extrapolation  messages  -  This  output 
permits  the  user  to  print  internal  data  arrays  (DUMP)  or  to 
investigate  design  chart  extrapolations  during  the  course  of 
the  case  execution  (PRINT  EXTRAP). 

The  remainder  of  the  section  describes  each  of  these  output  data. 
Examples  of  each  output  page  are  also  included  and  were  created  from  the 
example  problems,  described  in  Appendix  A,  which  can  be  used  as  a  model 
for  setting  up  another,  similar  configuration  or  be  used  as  a  means  to  check 
the  proper  operation  of  the  code. 

4.1  NOMINAL  OUTPUT 

Without  the  use  of  any  program  options  the  code  will  provide  three 
types  of  output.  First,  an  analysis  by  the  input  error  checking  routine  is 
provided.  It  lists  all  input  cards  provided  by  the  user  and  identifies  any  input 
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errors  detected.  Second,  a  listing  of  all  input  cards,  grouped  by  case,  are 
provided;  included  in  this  output  is  an  error  analysis  from  the  major  input 
error  routine  MAJERR.  Finally,  the  total  configuration  aerodynamics  are 
provided  in  summary  form;  one  page  of  aerodynamic  output  is  supplied  for 
each  Mach  number  specified.  The  MAJERR  results  and  the  total  configuration 
aerodynamics  results  are  listed  in  succession  for  each  case. 

4.1.1  Input  Error  Checking 

The  purpose  of  the  input  error  checking  module  is  to  provide  single 
pass  error  checking  of  all  inputs.  If  an  error  is  detected,  it  is  identified  and  an 
appropriate  error  message  provided.  The  error  messages  are  designed  to  be 
self-explanatory.  In  some  cases,  errors  are  automatically  corrected  by  the 
routine,  although  the  routine  was  not  designed  to  be  a  comprehensive  error 
correction  utility. 

The  following  errors  are  automatically  corrected  by  the  code: 

•  No  terminating  comma  on  a  namelist  input  card 

•  No  terminating  "$"  or  "$END"  on  a  namelist  input  ("&"  on 
IBM  systems) 

•  No  terminating  NEXT  CASE  for  the  case  inputs  for  single 
case  or  last  case  inputs. 

Errors  detected  by  the  error  checking  routine  are  considered  either 
"FATAL"  or  "NON-FATAL".  A  "FATAL"  error  is  one  which  will  cause  the 
code  to  terminate  execution  abnormally;  examples  of  "FATAL"  errors  include 
incorrect  spelling  of  any  namelist  name,  incorrect  spelling  of  any  variable 
name,  and  any  drastic  input  error  in  a  namelist  input,  such  as  leaving  out  an 
equals  sign  in  a  constant  definition.  All  "FATAL”  errors  are  clearly  identified 
on  the  output.  A  "NON-FATAL”  error  is  one  which  will  not  cause  the 
program  to  terminate  execution;  an  example  of  a  "NON-FATAL"  error  is 
leaving  off  the  decimal  point  on  numeric  constants  all  Missile  Datcom  inputs 
are  either  REAL  or  LOGICAL  regardless  of  the  variable  name  assigned. 
"NON-FATAL"  errors  will  not  cause  the  code  to  stop  execution,  whereas, 
"FATAL"  errors  will  cause  the  code  to  stop  execution  after  input  error 
checking  has  been  completed. 

An  example  output  from  CONERR  is  shown  in  Figure  21.  This  figure 
illustrates  the  array  of  input  errors  checked  by  CONERR.  Several  additional 
features  of  the  output  are  as  follows: 

•  All  user  defined  input  cards  are  assigned  a  sequential  "line 
number".  This  serves  to  identify  user  inputs  from  the  code 
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generated  inputs  (all  code-created  input  cards  are  not 
identified  with  a  "line  number”).  This  scheme  also  permits 
the  user  to  quickly  identify  those  input  cards  in  error  so  that 
efficient  correction  of  input  errors  can  be  performed. 

•  All  input  cards  are  listed  as  input  by  the  user.  To  the  right  of 
each  input  card  is  a  listing  of  any  errors  encountered  in 
processing  that  card.  If  no  such  error  message  appears  then 
the  input  was  interpreted  as  being  correct. 

•  In  many  cases  alphanumeric  constants  are  available  (see 
Table  3).  Hence  the  user  does  not  need  to  memorize  a 
numeric  scheme  of  "flags”.  Since  some  computers  do  not 
recognize  alphanumeric  constants  as  namelist  constants,  they 
are  automatically  converted  by  the  code  to  their  numeric 
equivalent.  A  message  is  printed  to  identify  the  substitutions 
performed.  The  example  input  ir  Figure  21  shows 
replacements  for  CONE  and  OGIVE. 

In  order  to  permit  column  independent  inputs  the  code  will 
automatically  adjust  some  of  the  input  cards  to  begin  in  columns  1  or  2.  All 
control  cards  will  be  automatically  shifted  to  start  in  column  1;  all  namelists 
which  begin  in  column  1  will  be  shifted  to  column  2.  If  any  input  card  cannot 
be  shifted  to  conform  to  this  scheme,  an  error  message  will  be  produced.  As  a 
general  rule,  column  80  of  namelist  inputs  should  be  left  blank  so  that  the 
code  can  shift  the  card  image,  if  necessary. 

4.1.2  Listing  of  Case  Input  Data 

Figure  22  shows  the  first  page  of  outputs  for  a  case  without  CONERR 
detected  errors.  Then  Figure  23  shows  the  next  page  of  output  which  lists  all 
input  cards  for  the  case  (down  to  the  NEXT  CASE  control  card).  If  the  input 
for  a  case  is  from  a  previous  case  (through  use  of  the  SAVE  control  card)  only 
the  new  case  inputs  are  listed.  All  saved  inputs  are  not  repeated  in  subsequent 
case  input  summaries. 

After  the  case  data  have  been  read,  the  data  set-up  for  the  case  is 
analyzed  by  the  case  major  error  checking  module  (MAJERR).  The  purpose  of 
this  second  error  checking  is  to  insure  that  the  data  input,  although  syntax 
error  free,  properly  defines  a  case  to  be  run.  Examples  of  errors  detected  in 
MAJERR  include  valid  flight  condition  inputs,  valid  reference  condition 
inputs,  and  that  geometry  has  been  defined.  In  most  cases  errors  detected  by 
MAJERR  are  corrected  with  assumed  defaults.  If  any  MAJERR  error  message 
is  produced,  the  user  should  verify  the  "fix-up"  taken  by  the  cede.  In  some 
cases  a  "fix-up"  is  not  possible;  an  appropriate  error  message  and  a  suggestion 
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for  correcting  the  error  is  provided.  If  a  "fix-up”  is  not  possible  the  case  will 
not  run. 

4.1.3  Case  Total  Configuration  Aerodynamic  Output  Summary 

As  shown  in  Figure  24,  the  total  configuration  aerodynamics  are 
provided  in  compact  form  for  easy  review.  The  aerodynamics  are 
summarized  as  a  function  of  angle  of  attack  (ALPHA)  in  the  user  specified 
system  of  units,  the  nomenclature  is  as  follows: 

CN  -  Normal  force  coefficient 

CM  -  Pitching  moment  coefficient 

CA  -  Axial  force  coefficient 

CY  -  Side  force  coefficient 

CLN  -  Yawing  moment  coefficient 

CLL  -  Rolling  moment  coefficient 

CNA  -  Normal  force  coefficient  derivative  with  ALPHA 

CMA  -  Pitching  moment  coefficient  derivative  with 

ALPHA 

CYB  -  Side  force  coefficient  derivative  with  BETA 

CLNB  -  Yawing  moment  coefficient  derivative  with  BETA 

CLLB  -  Rolling  moment  coefficient  derivative  with  BETA 

CL  -  Lift  coefficient 

CD  -  Drag  coefficient 

CL/CD  -  Lift  to  drag  ratio 

XCP  -  Center  of  pressure  from  the  moment  reference 

center  divided  by  reference  length 

All  coefficients  are  based  upon  the  reference  areas  and  lengths  specified 
at  the  top  of  the  output  page.  The  derivatives  CNA  and  CMA  are  computed 
by  numeric  differentiation  of  the  CN  and  CM  curves,  respectively;  precise 
derivatives  are  only  obtained  when  the  angle  of  attack  range  specified  is 
narrow.  The  derivatives  CYB,  CLNB  and  CLLB  are  determined  by  perturbing 
the  sideslip  angle  by  one  degree,  recalculating  the  configuration  forces  and 
moments,  and  then  differencing  with  the  user  specified  orientation.  Hence, 
the  longitudinal  and  lateral  derivatives  will  probably  not  be  numerically 
identical  for  those  conditions  which  should  produce  identical  results  if  they 
were  both  calculated  by  the  same  method. 

A  significant  decrease  in  computational  time  is  realized  when  the 
calculation  of  lateral-directional  derivatives  are  suppressed  using  the  control 
card  NO  LAT.  For  these  cases,  the  CYB,  CLNB,  and  CLLB  data  fields  are  filled 
with  blanks. 

When  selecting  TRIM,  the  output  is  provided  in  a  form  similar  to 
Figure  25  When  running  a  trim  case  the  derivatives  due  to  ALPHA  and 
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BETA  are  not  available.  The  panels  which  were  deflected  to  trim  the 
configuration  are  indicated  by  the  "VARIED"  citation  next  to  them. 

The  format  for  the  values  of  the  numbers  in  the  printed  output  has 
been  assumed  based  on  typical  magnitudes  for  missile  aerodynamic 
coefficients.  In  some  cases,  a  user  specified  reference  area  and/or  length  will 
cause  the  results  to  underflow  or  overflow  the  format  selected.  For  these  cases 
the  user  should  adjust  his  reference  quantities  by  powers  of  ten  to  get  the  data 
to  fit  the  format  specified. 

4.2  PARTIAL  OUTPUT 

Partial  output  consists  of  geometry  calculation  details,  intermediate 
aerodynamic  results,  or  auxiliary  data,  such  as  pressure  distributions.  Each  of 
these  output  types  are  printed  through  the  addition  of  control  cards  input  for 
each  case.  In  all  cases,  partial  output  requested  for  one  case  is  not 
automatically  selected  for  subsequent  cases,  and  the  control  cards  must  be  re¬ 
input.  This  permits  the  user  to  be  selective  on  the  amount  and  types  of 
output  desired. 

A  special  control  card  PART  permits  the  user  to  request  all  geometric 
and  aerodynamic  partial  output.  Due  to  the  amount  of  output  produced,  this 
option  should  be  used  sparingly  or  when  details  of  the  calculations  are 
desired. 

The  following  paragraphs  describe  the  output  received  when  partial 
output  is  requested. 

4.2.1  Geometric  Partial  Output 

Details  of  the  geometry  are  provided  when  the  PART  or  PRINT  GEOM 
control  cards  are  included  in  the  case  inputs.  Figure  26  shows  the  output 
created  when  the  PRINT  GEOM  BODY  control  card  is  used.  Detailed  are  the 
results  of  the  geometric  calculations  for  the  body.  Included  are  such  items  as 
planform  area,  surface  (wetted)  area,  and  the  mold  line  contour. 

If  fins  are  present  on  the  configuration,  two  types  of  fin  geometry  data 
are  produced  when  PRINT  GEOM  FIN1  or  PART  is  requester  As  shown  in 
Figure  27,  the  description  of  the  panel  airfoil  section  is  provided.  Following 
that,  shown  in  Figure  28,  is  a  summary  of  the  major  geometric  characteristics 
of  such  planform;  note  that  fin  planform  geometry  data  is  given  for  one  panel 
of  each  fin  set,  since  it  is  assumed  that  each  fin  of  a  fin  set  is  identical.  If  a 
panel  is  made  up  of  multiple  segments,  the  geometric  data  is  provided  by 
panel  segment  (each  segment  is  assigned  a  number  starting  at  the  root).  Total 
panel  set  of  characteristics  is  also  provided.  This  total  panel  data  represents  an 


73 


equivalent  straight- tapered  panel,  which  is  used  for  most  of  the  aerodynamic 
calculations.  The  thickness-to-chord  ratio  shown  for  each  segment  is  that 
value  at  the  segment  root;  for  the  total  panel,  it  is  an  "effective"  value. 

If  an  airbreathing  inlet  is  specified  the  output  is  similar  to  that  in 
Figure  29.  This  output  reflects  the  user  input  definition  for  the  inlet  design 
specified.  It  is  provided  if  the  PRINT  GEOM  INLET  or  PART  control  cards  are 
included  in  the  input  case. 

4.2.2  Aerodynamic  Partial  Output 

The  output  on  the  configuration  aerodynamics  is  most  extensive  when 
PRINT  AERO  or  PART  is  specified.  Output  is  created  for  the  body  and  each  fin 
set  on  the  configuration.  In  addition,  for  any  subsonic /transonic  Mach 
number  (less  than  1.4)  an  analysis  by  the  Airfoil  Section  Module  is  made, 
which  involves  a  potential  low  analysis  of  the  airfoil  section  using  conformal 
mapping.  If  a  configuration  has  inlets  additional  partial  output  is  included  to 
summarize  the  inlet  external  aerodynamics. 

If  base-jet  plume  interaction  calculations  are  specified  (BASE=.TRUE. 
in  namelist  AXIBOD),  then  there  will  be  one  or  two  separate  pages  of  output. 
Figure  30  shows  an  example  of  the  first  page  of  output.  This  page  will  always 
be  printed  if  BASE=.TRUE.  The  base  pressure  coefficient,  axial  force 
coefficient,  and  freestream  pressure  and  temperature  ratios  are  shown  versus 
angle  of  attack.  Also,  the  incremental  forces  and  moments  due  to  separation 
are  shown  versus  angle  of  attack.  If  extrapolation  of  the  base  pressures  and 
separation  conditions  database  occurs,  a  warning  message  is  printed  explaning 
what  input  variable  required  extrapolation.  A  second  page  of  output 
containing  the  boattail  separation  parameters  will  be  printed  if  there  are  any 
fins  on  the  missile  boattail.  The  separation  location  aft  of  the  nose  and  the 
Mach  cone  angle  are  shown  versus  angle  of  attack  for  each  panel  on  the  fin 
set.  Figure  31  shows  an  example  of  this  page.  This  output  is  provided  if  the 
PRINT  AERO  BODY  or  PART  control  card  is  input. 

The  protuberance  partial  output  is  printed  if  PRINT  AERO  BODY  or 
PART  is  used.  This  output  will  only  be  shown  if  the  namelist  PROTUB  is 
present  in  the  input  file.  Figure  32  is  an  example  of  the  protuberance  output. 
Protuberance  type,  location,  number,  and  axial  force  coefficient  are  listed  for 
each  protuberance  set.  The  total  axial  force  coefficient  or  zero  lift  drag 
coefficient  is  printed  at  the  bottom  of  the  page. 

As  shown  in  Figure  33,  the  body  alone  partial  aerodynamic  output  for 
normal  force  lists  the  axial  force  contributors,  potential  normal  force  (CN- 
POTENTIAL),  viscous  normal  forces  (CN-VISCOUS),  potential  pitching 
moment  (CM-POTENT1AL),  viscous  pitching  moment  (CM-V1SCOUS),  and 
the  crossflow  drag  coefficient  (CDC).  The  cross-flow  drag  proportionality  factor 
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at  subsonic  and  transonic  speeds  is  also  given  for  reference.  These  data  are 
similar  to  that  obtained  for  elliptical  bodies. 

Figure  34  details  the  fin  normal  force  calculations  by  fin  set.  Each 
panel's  contribution  to  the  configuration  normal  force  is  described.  The 
column  titled  CN-POTENTIAL  is  the  potential  contribution  and  the  column 
titled  CN-VISCOUS  is  the  viscous  contribution.  Their  sum  is  given  in  the 
column  titled  CN-TOTAL.  CNAA  is  the  nonlinear  variation  of  normal  force 
due  to  angle  of  attack  and  ALPHA  EQUIV  is  the  panel  angle  of  attack  due  to 
its  roll  position  on  the  body.  Figure  35  illustrates  the  fin  axial  force 
contributors  and  Figure  36  presents  an  example  of  the  fin  pitching  moment 
contributors. 

The  analysis  by  the  Airfoil  Section  Module  is  provided  in  a  format 
similar  to  Figure  37.  If  any  Mach  number  specified  produces  supersonic  flow 
on  the  airfoil  surface,  the  message  "CREST  CRITICAL  MACH  NUMBER 
EXCEEDED"  will  be  printed;  approximation  of  the  airfoil  section  data  is  then 
assumed.  These  fin  aerodynamic  increments  are  repeated  for  each  fin  set  on 
the  configuration.  Note  that  the  Airfoil  Section  Module  assumes  that  the 
panels  have  sharp  trailing  edges.  Any  panel  input  with  a  non-sharp  trailing 
edge  will  have  its  aerodynamic  characteristics  set  as  though  the  airfoil  was 
"ideal".  This  assumption  is  approximate  for  preliminary  design. 

Figure  38  shows  the  aerodynamic  output  available  when  inlets  are 
specified  on  the  configuration.  It  is  provided  when  PRINT  AERO  INLET  or 
PART  is  specified  in  the  case  inputs.  The  aerodynamics  summarized  for  inlets 
can  include  additive  drag  results  if  the  user  input  the  additive  drag 
calculation  flag.  The  maximum  mass  flow  ratio  is  printed  at  the  bottom  of  the 
page  if  the  additive  drag  is  calculated.  If  additive  drag  cannot  be  calculated,  a 
warning  message  is  printed. 

After  the  aerodynamic  details  for  each  component  of  the  configuration 
are  output,  the  aerodynamic  calculations  for  the  synthesis  of  the  complete 
configuration  follows.  For  the  example  case,  fin  set  1  results  would  be 
followed  by  fin  set  2  results  for  each  of  the  following  outputs: 

•  "FIN  SET  PRESENCE  OF  THE  BODY"  -  This  summarizes  the 
aerodynamic  incrementals  of  the  most  forward  set  of  fins 
with  the  influence  of  the  body.  Figure  39  presents  the 
example  of  this  output.  The  left-most  six  columns  include 
the  effect  of  body-on-fin  component  interference.  The  right¬ 
most  columns  represent  the  contribution  to  each  panel  to 
configuration  aerodynamics,  and  include  the  effect  of  body- 
on-fin  interference,  these  values  are,  in  effect,  individual 
panel  loads.  The  panel  characteristic  values  included  are 
"AEQn"  (the  panel  equivalent  (local)  angle  of  attack)  and 


75 


"CNn"  (the  panel  normal  force  coefficient).  The  sign 
convention  is  as  follows:  a  positive  panel  normal  force, 
hence,  equivalent  angle  of  attack,  produces  a  negative  roll 
moment.  Therefore,  panels  on  the  right  side  of  the 
configuration  will  produce  loads  and  angles  of  attack  opposite 
in  sign  to  those  on  the  left  side  of  the  configuration  even 
though  they  produce  the  same  physical  force  loading. 

•  "BODY-FIN  SET"  -  Aerodynamics  for  the  body  plus  most 
forward  set  of  fins  configuration.  It  is  produced  through 
addition  of  the  body  alone  and  wing  in  presence  of  the  body 
incrementals,  described  above.  The  results  shown  in  Figure 
40,  include  the  component  carryover  factors  K-W(B)  (wing  in 
presence  of  the  body  carryover  due  to  angle  of  attack),  K-B(W) 

(body  in  presence  of  the  wing  carryover  due  to  angle  of 
attack),  KK-B(W)  (body  in  presence  of  the  wing  carryover  due 
to  panel  deflection),  XCP-W(B)  (wing  in  presence  of  the  body 
carryover  center  of  pressure),  and  XCP-B(W)  (body  in 
presence  of  the  wing  carryover  center  of  pressure).  This 
output  is  repeated  for  the  body  plus  each  additional  aft  fin  set, 
if  one  exists  on  the  configuration.  Th:s  example  includes  two 
fin  sets  so  the  next  page  of  partial  output  would  look  like 
Figure  41.  If  additional  fin  sets  are  present  on  the 
configuration  additional  pages  are  output  with  each  one 
successively  included. 

•  "CARRYOVER  INTERFERENCE  FACTORS"  -  This  page  of 
partial  output  summarizes  the  carryover  factors  listed  in  the 
paragraph  above.  These  were  included  in  the  body  plus  fln  set 
calculations.  An  example  of  this  output  is  presented  in  Figure 
42. 

•  "COMPLETE  CONFIGURATION"  -  Complete  configuration 
aerodynamics.  This  output  was  illustrated  in  Figure  24.  The 
values  are  obtained  by  summing  the  body-wing  and  tail  in 
the  presence  of  the  wing  flow  field  data. 

In  addition  to  the  output  described  above,  more  data  is  presented  when 
the  BUILD  control  card  is  used.  Static  aerodynamics  are  output  for  each 
configuration  component.  Body  alone  aerodynamics  are  shown  in  Figure  43. 
Fin  alone  aerodynamics  are  shown  for  each  fin  set  present.  Figure  44  shows 
the  output  for  the  first  fin  set.  Static  aerodynamics  for  a  configuration  with 
body  plus  most  forward  set  of  fins  is  given  next.  Figure  45  shows  an  example 
of  this  output.  This  output  is  repeated  for  configurations  including  the  body 
plus  each  additional  fin  set  present. 
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If  the  PRINT  AERO  BEND  or  PART  control  card  is  used,  the  code  will 
compute  and  print  panel  bending  moment  coefficients  for  each  fin  set  on  a 
separate  page.  One  page  is  shown  in  Figure  46.  The  sign  convention  is  that 
assumed  for  the  individual  panel  loads  and  equivalent  angles  of  attack,  noted 
above.  The  bending  moment  coefficients  are  based  upon  the  reference  area 
and  longitudinal  length  given  at  the  top  of  the  page.  The  moments  are 
referenced  about  the  fin-body  structure  specified  by  the  root  chord  span 
station. 

Figure  47  illustrates  the  panel  hinge  moments  coefficients  computed 
when  the  control  cards  PRINT  AFRO  HINGE  or  PART  are  used.  The 
reference  area  and  longitudinal  reference  length  given  at  the  top  of  the  page 
are  used.  All  moments  are  computed  about  the  hinge  line,  which  is  defined 
using  namelist  DEFLCT. 

If  TRIM  is  specified,  the  user  can  selectively  print  the  six  untrimmed 
static  aerodynamic  tables  used  in  the  trim  process.  An  example  is  shown  in 
Figure  48.  The  code  computes  the  six-component  aerodvnamics  at  ten 
deflection  angles  for  each  specified  angle  of  attack,  then  interpolates  for  Cm=0. 
Note  that  this  trim  process  can  be  used  to  create  control  authority  data, 
effectively  giving  the  user  10  deflection  angles,  20  angles  of  attack,  and  20 
Mach  numbers  per  input  case. 

4.2.3  Pressure  Distribution  Data 

If  the  Mach  number  is  supersonic  (M  >  1.2),  the  user  has  the  option  to 
print  the  surface  pressure  distributions  over  the  body  and  fins.  This  option  is 
selected  only  through  the  addition  of  the  control  card  PRESSURES.  Since 
three  body  alone  supersonic  methods  are  available  (Van  Dyke  Hybrid, 
Second-Order  Shock  Expansion  (SOSE),  and  Newtonian  flow)  the  capability 
exists  to  output  the  pressure  distribution  data  from  any  one  of  these  methods. 
The  method  to  be  used  in  the  calculation  of  the  pressure  data  is  controlled 
with  the  control  cards  SOSE  and  HYPER;  if  neither  control  card  is  input,  the 
Van  Dyke  Hybrid  method  is  selected.  Because  of  the  nature  of  the 
calculations,  body  alone  pressures  are  printed  for  angles  of  attack  less  than  or 
equal  to  15  degrees  when  using  the  Hybrid  or  SOSE  techniques. 

The  capability  also  exists  for  the  user  to  output  the  pressure 
distribution  data  over  fins  at  any  Mach  number  greater  than  1.05.  This  option 
is  also  controlled  by  the  PRESSURES  control  card.  Due  to  the  nature  of  the 
method,  only  pressure  distribution  data  at  zero  angle  of  attack  is  presently 
output. 

Figures  49,  50,  and  51  illustrate  typical  output  produced  when 
PRESSURES  is  specified.  The  format  of  Figure  49  is  only  available  when  SOSE 
is  specified;  all  other  body  alone  pressure  methods  produce  output  similar  to 
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Figure  50  for  bodies.  Figure  51  is  representative  of  fin  pressure  distribution 
output.  Note  that  calculation  of  pressures  is  a  time-consuming  process;  much 
higher  computational  times  will  be  required. 

All  body  pressure  distribution  data  is  based  on  a  configuration  that  has 
body  diameter  of  unity;  that  is,  the  configuration  is  expressed  in  calibers  (or 
body  diameters).  The  longitudinal  stations  at  which  pressure  coefficient  data 
is  desired  cannot  be  user  specified;  however,  sufficient  data  is  provided  to 
permit  accurate  interpolation  for  most  applications. 

4.3  DYNAMIC  DERIVATIVES 

As  shown  in  Figure  52,  the  total  configuration  dynamic  derivatives  are 
provided  in  compact  form  for  easy  interpretation.  The  dynamic  derivatives 
are  summarized  as  a  function  of  angle  of  attack  in  the  user  specified  units. 
The  coefficients  provided  are  as  follows: 

CNQ  Normal  force  coefficient  due  to  pitch  rate 

CNAD  Normal  force  coefficient  due  to  rate  of  change  of 
angle  of  attack  (a) 

CMQ  Pitching  moment  coefficient  due  to  pitch  rate 

CMAD  Pitching  moment  coefficient  due  to  rate  of  change 
of  angle  of  attack  (a) 

Note:  For  body  alone  and  body  +  fin  set  data  CMQ  and  CMAD  are 
presented  as  the  sum  CMQ+CMAD. 

The  dynamic  derivatives  are  printed  after  all  static  coefficients  and 
partial  static  aerodynamics  are  printed.  If  a  BUILD  or  PART  card  is  input, 
additional  dynamic  derivatives  for  partial  configurations  and/or 
configuration  components  are  printed. 

4.4  EXTERNAL  DATA  FILES 

The  code  has  the  capability  to  be  used  in  conjunction  with  other 
missile  design  tools,  such  as  post-processing  plotting  programs  or  trajectory 
programs.  Fixed  format  aerodynamic  data  is  output  as  an  external  data  file 
with  the  addition  of  the  PLOT  control  card.  Included  in  this  data  file  are  the 
six  component  forces  and  moments  based  upon  the  user  specified  reference 
quantities.  In  order  to  print  component  buildup  data  to  the  plot  file  the 
DUILD  and  PLOT  control  cards  must  be  present  in  the  case. 
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An  option  to  create  a  user  specified  format  data  file  is  also  available. 
The  control  cards  WRITE  and  FORMAT  have  been  designed  for  easy  access  to 
this  capability. 

Both  the  PLOT  output  and  output  generated  via  the  WRITE  control  are 
written  to  unit  3.  Thus,  if  both  PLOT  and  WRITE  are  used  in  the  same  run, 
the  external  data  file  will  have  both  output  formats  in  the  same  file. 

4.5  EXTRAPOLATION  MESSAGES  AND  ARRAY  DUMPS 

As  shown  in  Figure  53,  the  extrapolation  messages  are  summarized  for 
all  design  charts  which  have  been  extrapolated  during  the  execution  of  the 
case.  Since  many  of  the  aerodynamic  methods  do  not  include  design  charts, 
but  are  either  closed-form  equations  or  complete  theoretical  methods,  this 
option  is  most  useful  in  the  subsonic  and  transonic  Mach  regimes. 
Extrapolation  messages  are  only  provided  if  the  control  card  PRINT  EXTRAP 
appears  in  the  case  inputs.  The  data  titled  "ROUTINE  TRACE-BACK"  lists  the 
subroutines  called  when  the  look-up  was  performed;  "X”  represents  the 
independent  variable  and  "Y"  represents  the  dependent  variable  in  the 
extrapolation. 

When  it  is  necessary  to  examine  the  values  stored  in  internal  data 
arrays  the  DUMP  control  card  can  be  used.  This  control  card  causes  the 
contents  of  the  named  data  arrays  to  be  printed  in  a  form  similar  to  Figure  54. 
Array  dumps  are  provided  for  each  Mach  number  of  the  input  case,  and 
represent  the  data  block  contents  at  aerodynamic  calculation  completion. 

Note  that  all  data  arrays  are  initialized  to  a  constant  named 
"UNUSED",  which  is  preset  to  a  value  of  1  x  10'30.  Hence,  any  array  element 
which  contains  this  constant  was  not  changed  during  execution  of  the  case 
(since  it  is  highly  unlikely  that  this  constant  will  result  from  any  calculation). 
This  scheme  permits  rapid  "tracking"  of  program  calculation  sequences  while 
in  "debug"  mode. 
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Figure  23  Example  of  Default  Substitutions 
for  Incomplete  Case  Inputs 
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Figure  24  Total  Configuration  Aerodynamic  Output  Summary 
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Figure  25  Trimmed  Output  Summary 
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Figure  26  Body  Geometry  Output 
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Figure  27  Airfoil  Geometry  Output 
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Figure  28  Fin  Geometry  Output 
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Figure  31  Base-Jet  Plume  Interaction  Output  -  Page 
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Figure  32  Protuberance  Output 
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Figure  33  Body  Alone  Aerodynamic  Partial  Output 
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Figure  34  Fin  Normal  Force  Partial  Output 
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Figure  35  Fin  Axial  Force  Partial  Output 
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Figure  36  Fin  Pitching  Moment  Partial  Output 
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Figure  38  Inlet  Aerodynamic  Partial  Output 
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Figure  39  Fin  Set  in  Presence  of  the  Body  Partial  Output 
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Figure  40  Body  Plus  Fin  Set  Partial  Output 
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Figure  41  Body  Plus  Two  Fin  Sets  Partial  Output 
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Figure  42  Carryover  Interference  Factors  Partial  Output 
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Figure  43  Body  Alone  Static  Aerodynamic  Partial  Output 
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Figure  44  Fin  Alone  Static  Aerodynamic  Partial  Output 
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Figure  45  Body  Plus  Fin  Static  Aerodynamic  Partial  Output 
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Figure  47  Panel  Hinge  Moment  Partial  Output 
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Figure  48  Untrimmed  Partial  Output 
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Figure  50  Body  Pressure  Distribution  at  Angle  of  Attack 
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Figure  51  Fin  Pressure  Distribution  Output 
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Figure  52  Dynamic  Derivative  Output 
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5.0  IMPLEMENTATION  GUIDE 


This  section  details  the  steps  necessary  to  make  the  code  operational  on 
the  user's  computer  system.  In  addition  this  section  contains  cross  reference 
tables  for  all  the  subroutines  and  common  blocks  used  in  the  computer 
program.  It  also  contains  a  description  of  the  variables  in  the  common  blocks 
that  can  be  written  out  using  the  WRITE  or  DUMP  cards  in  the  input  deck. 

5.1  INSTALLATION  ON  COMPUTER  SYSTEMS 

This  section  details  the  steps  necessary  to  make  the  computer  code 
functional  on  the  user's  computer  system.  Although  conversion  of  the 
program  can  be  easily  accomplished  by  someone  with  a  good  understanding 
of  the  FORTRAN  V  language,  it  is  highly  recommended  that  someone 
familiar  with  the  computer  operating  system  be  consulted. 

5.1.1  Requirements 

In  order  for  the  Missile  Datcom  code  to  be  successfully  implemented  on 
the  user's  computer  system,  there  are  three  requirements  which  must  be  met, 
as  follows: 

•  Language  -  As  received,  the  CDC  code  is  compatible  with 
FORTRAN  IV,  except  as  noted  below.  The  program  can  be  easily 
converted  to  FORTRAN  V  and  the  changes  required  are  described 
later  in  this  section.  The  VAX  compatible  code  is  maintained  as  a 
FORTRAN  V  version. 

•  Namelist  -  The  code  has  been  designed  with  an  internal 
FORTRAN  NAMELIST  emulator  to  allow  the  input  and  output 
(I/O)  to  be  handled  by  namelist  variables.  This  is  an  exception  to 
Standard  FORTRAN  but  with  the  emulator  as  part  of  the  code  the 
program  will  run  under  Standard  FORTRAN.  The  code  is  not 
easily  converted  to  fixed  field,  rather  than  namelist  input. 

•  I/O  Scratch  Files  -  The  code  uses  the  following  logical  file  units:  1, 
2,  3,  4,  5,  6  and  7.  All  file  units  are  accessed  using  formatted  reads 
and  writes.  File  units  1,  2,  4  and  7  are  used  internally;  file  units  3, 
5  and  6  transfer  data  between  the  user  and  the  code. 
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5.1.2  Coding  Changes 

In  order  for  the  code  to  made  operational  on  other  computer  systems, 
several  minor  coding  changes  may  be  required. 

•  Namelist  Delimiter  -  The  namelist  delimiter  used  on  the 
computer  system  is  set  in  the  routine  BLOCK  DATA.  It  is  preset  to 
the  character  for  VAX  computers  in  the  DATA  statement 
defining  the  variable  KAND.  It  must  be  changed  for  other 
computers,  which  employ  the  "&"  character  for  the  namelist 
delimiter. 

•  Defining  UNUSED  -  The  variable  UNUSED  is  defined  in  routine 
BLOCK  DATA.  It  is  preset  to  "l.E-30"  and  should  be  acceptable  for 
most  computers.  If  this  number  is  too  small  for  the  computer 
being  used,  it  may  be  changed  to  a  small,  non-zero,  positive 
number  which  is  no  larger  than  1  x  10'10.  This  constant  is  used  to 
initialize  all  program  data  arrays.  If  a  new  value  is  assigned  it 
must  also  be  incorporated  in  NMLIST,  variable  JUNSD. 

5.1.3  Input/Output 

Seven  file  units  are  used  by  the  program.  They  are  the  logical  units  1, 

2,  3,  4,  5,  6  and  7  which  are  used  as  follows: 

Unit  Usage 

1  All  input  cards  read  from  unit  5  are  written  to  unit  1  by 
CONERR  after  they  have  been  checked  for  errors 

2  Namelists  for  the  input  "case"  are  read  from  unit  1  and  written 
to  unit  2  by  READIN.  The  namelists  for  the  "case"  are  read 
from  unit  2. 

3  Data  dumped  at  user  request  (using  PLOT  or  WRITE  control 
cards)  are  written  to  unit  3  by  PLOT3,  PLTTRM,  and/or  SAVEF. 
The  PLOT  file  format  is  given  in  Appendix  B. 

4  Method  extrapolation  messages  are  printed  in  compact  form  on 
unit  4  by  MESSG  and  read  by  EXTRAP. 

5  User  inputs  are  read  from  unit  5  by  CONERR. 

6  Program  output  is  written  to  unit  6. 

7  The  FORMAT  and  WRITE  control  cards  are  written  to  unit  7  by 
CONTRL  and  read  by  SAVEF. 

The  method  extrapolation  messages  are  written  to  unit  4  as  they  are 
generated  if  a  PRINT  EXTRAP  control  card  is  present.  When  the 
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aerodynamic  calculations  are  completed  for  each  Mach  number,  the  line  END 
MACH  is  written  to  signify  the  end  of  messages  for  that  Mach  number.  For 
each  extrapolation  two  lines  are  written  to  unit  4  as  follows: 


Line 

1 

2 


Contents 

The  names  of  the  routines  called  to  get  to  this  part  of 
the  code.  The  last  name  is  the  routine  currently 
being  executed. 

The  seven  (7)  parameters  defining  the  extrapolation: 
First  Value:  Value  for  the  independent  variable 

Second  Value:  Final  Result 

Lower  independent  value 
Upper  independent  value 
Lower  dependent  value 
Upper  dependent  value 


Third  Value: 
Fourth  Value: 
Fifth  Value: 
Sixth  Value: 


Seventh  Value:  Order  of  extrapolation 


The  data  is  automatically  read  and  formatted  for  output  by  the  program 
in  the  standard  output  file.  Hence,  this  file  need  not  be  retained  or  printed  by 
the  user  at  the  end  of  code  execution. 


5.2  PROGRAM  CROSS  REFERENCES 

This  section  summarizes  the  cross  references  of  subroutines  and 
common  blocks  within  the  code.  The  tables  given  can  be  used  to  verify  the 
proper  location  of  common  blocks  within  the  code,  identify  the  subroutines 
called  by  each  of  the  methods  and,  if  program  revisions  are  made,  facilitates 
subroutine  modification.  The  cross  references  are  given  by  subroutine  name 
and  by  common  block  name. 

5.2.1  Subroutine  Cross  Reference  By  Name 

Figures  55  and  56  summarizes  the  subroutine  cross  references  for  the 
program.  Figure  55  lists  each  subroutine  and  gives  the  subroutines  that  call  it. 
Figure  56  list  each  subroutine  and  gives  the  subroutines  which  it  calls.  For 
convenience,  the  subroutines  are  listed  in  alphabetical  order. 

5.2.2  Common  Block  Cross  Reference  By  Name 

Figures  57  and  58  summarize  the  common  blocks  used  by  the  program. 
Seven  types  of  common  blocks  are  used  by  the  program:  (1)  Internal  data 
management,  (2)  Input  data,  (3)  Aerodynamic  work,  (4)  Geometry,  (5)  Static 
aerodynamic,  (6)  Dynamic  aerodynamic,  and  (7)  Trim  blocks.  These  common 
blocks  are  shown  in  Figure  59  for  each  type  of  data.  For  each  block  type,  the 
data  is  subdivided  by  task  or  configuration,  as  appropriate. 
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Internal  Data  Management  :  These  eight  common  blocks  contain  the 
data  to  control  program  execution. 

CASEID  -  Contains  the  case  title  (  from  the  CASEID  control  card), 
the  case  number  and  flag  to  print  extrapolation  messages  and  the 
flag  to  suppress  calculation  of  the  lateral-directional  derivatives. 

THERY  -  Contains  the  flags  necessary  to  select  the  correct 
supersonic  body  theory. 

CONST  -  Contains  the  program  constants. 

TRACE  -  Contains  the  position  in  the  code  being  executed. 

LOGIC  -  Contains  the  configuration  and  program  option  flags. 

DUMPF  -  Contains  the  logic  flags  to  delete  namelist  inputs  using 
the  DELETE  control  card. 

INPCON  -  Contains  the  data  required  to  validate  namelist  name 
inputs. 

Input  Common  Blocks:  Except  for  namelist  EXPR,  each  input  namelist 
is  assigned  a  separate  common  block  for  data  storage.  The  values  read 
in  from  EXPR  are  stored  in  the  common  blocks  in  which  the 
coefficients  of  the  partial  configurations  for  which  the  experimental 
data  is  provided  are  stored.  The  theoretical  values  are  over  written  by 
the  experimental  data. 


Common  Block 

Namelist 

fix: 

FLTCON 

REFQN 

REFQ 

ABODIN 

AXIBOD 

ABODIN 

ELLBOD 

FSET1 

FINSET1 

FSET2 

FINSET2 

FSET3 

FINSET3 

FSET4 

FINSET4 

DESIG 

FINSET1  -  FINSET4(NACA  control  card) 

INCID 

DEFLCT 

TRIMIN 

TRIM 
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Aerodynamic  Work  Arrays:  These  common  blocks  contain  most  of  the 
intermediate  aerodynamic  calculations,  particularly  those  needed  for 
other  methods. 


Common  Block 

BDWORK 

F1WORK 

F2WORK 

F3WORK 

F4WORK 


Configuration  Component 

Body 

Fin  Set  1 

Fin  Set  2 

Fin  Set  3 

Fin  Set  4 


Geometry  Common  Blocks:  These  common  blocks  hold  the  results  of 
the  geometric  calculations  for  each  component  of  the  configuration. 


Common  Block 

GEOBOD 

GEOFS1 

GEOFS2 

GEOFS3 

GEOFS4 


Configuration  Component 

Body 

Fin  Set  1 

Fin  Set  2 

Fin  Set  3 

Fin  Set  4 


Static  Aerodynamic  Results:  The  final  static  aerodynamic  results  are 
stored  by  configuration  or  configuration  component.  These  arrays 
form  the  Ideal  Output  Matrix  (I.O.M.),  where  the  data  is  stored  in  a 
fixed  format,  regardless  of  configuration  component.  These  data  are 
exactly  that  which  is  printed  as  the  normal  aerodynamic  output. 
Twenty  array  elements  are  reserved  for  each  aerodynamic  parameter, 
which  corresponds  to  the  20  angles  of  attack.  Eleven  aerodynamic 
coefficients  are  defined,  this  pattern  is  used  for  both  the  static  (SBODY, 
SFIN1,  SFIN2,SFIN3,  SFIN4,  SB1,  SB12,  SB123,  SB1234)  and  dynamic 
(DBODY,  DFIN1,  DFIN2,  DFIN3,  DFIN4,  DB1,  DB12,  DB123,  BD1234) 
aerodynamic  results.  Since  a  maximum  of  four  sets  of  fins  are 
permitted,  SB1  (and  DB1)  refers  to  the  body  plus  the  most  forward  fin 
set,  SB  12  (and  DB12)  refers  to  the  body  plus  the  first  and  second  most 
forward  fin  sets,  SB123  (and  DB123)  refers  to  the  body  plus  the  three 
most  forward  fin  sets.  Finally,  SB1234  (and  DB1234)  refers  to  the  body 
plus  all  four  fin  sets. 

Trim  Results:  The  code  uses  ten  deflection  angles,  as  a  function  of 
angle  of  attack,  in  order  to  interpolate  for  the  longitudinal  trim  points. 
The  untrimmed  aerodynamic  results  obtained  are  stored  in  the 
common  blocks  UTRIMD  as  a  function  of  angle  of  attack  and  deflection 
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angle,  the  aerodynamic  coefficients  retained  are  Cn,  Cm,  Ca,  Cy,  Cn, 
and  Q.  The  trimmed  results  are  stored  in  the  common  block  TRIMD, 
and  are  the  trimmed  values  for  Cm,  Ca,  Cy,  Cn,  and  Q. 

Experimental  Data  Substitution:  Experimental  substitution  allows  the 
user  to  input  experimental  data  for  any  part  of  the  configuration  or  any 
partial  configuration.  For  example  the  user  may  input  body  alone 
experimental  data,  fin  alone  experimental  data  and/or  body  +  one  (1) 
fin  set  experimental  data  for  a  configuration  having  a  body  +  two  (2) 
sets  of  fins.  The  experimental  is  substituted  into  the  appropriate 
common  block  and  replaces  the  computed  theoretical  coefficients.  The 
substituted  coefficients  are  then  used  in  the  configuration  synthesis 
process  in  subroutine  SYNTHS. 

Configuration  Incrementing:  Configuration  incrementing  requires  all 
cases  of  an  incrementing  run  to  have  the  same  configuration 
components  (i.e.  if  the  first  case  is  a  body  +  one  (1)  fin  set  then  all  the 
cases  that  follow  must  be  for  a  body  +  one  (1)  fin  set).  Configuration 
incrementing  runs  must  be  thought  of  as  a  sequence  of  case  runs  rather 
than  individual  cases.  When  the  first  case  is  executed,  the 
experimental  data  is  substituted  for  the  configuration  being  run.  The 
experimental  data  is  then  compared  with  the  theoretical  coefficients. 
The  correction  factors  that  are  computed  are  stored  in  the  INC 
common  block.  When  the  subsequent  cases  are  run,  the  correction 
factors  are  applied  to  the  theoretical  coefficients  in  SYNTHS.  These 
corrections  are  made  prior  to  printing  out  the  aerodynamic  results. 
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Figure  55  Subroutine  Cross  Reference  (Listed  by  Calling  Routine) 
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BODYCA  BOCAB  BDCALP  BDCAPR  BDCAWF  &DCAWN  BDCDRV 
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Figure  55  Subroutine  Cross  Reference  (Listed  by  Calling  Routine)  -  Continued 
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Figure  55  Subroutine  Cross  Reference  (Listed  by  Calling  Routine)  -  Continued 
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Figure  55  Subroutine  Cross  Reference  (Listed  by  Calling  Routine)  -  Continued 
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Figure  58  Common  Block /Subroutine  Cross  Reference  (Listed  by  Subroutine)  -  Continued 
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Internal  Common  Management  Blocks 

-  CASEID 

-  CONST 

-  DUMPF 

-  D FLAGS 

-  LOGIC 

-  PARTF 

-  DMPCON 

-  TRACE 

-  THERY 

Input  Common  Blocks 

-  ABODIN 

-  DESIG 

-  DIVERN 

-  FLC 

-  FSET1 

-  FSET2 

-  FSET3 

-  FSET4 

-  INQD 

-  INLETN 

-  REFQN 

-  TOTALC 

-  TRIM  IN 

Internal  Aerodynamic  Work  Arrays 

-  BDWORK 

-  FI  WORK 

-  F2WORK 

-  F3WORK 

-  F4WORK 

-  INC 

-  PAERO 

Geometry  Arrays 

-  GEOBOD 

-  GEORS1 

-  GEOR52 

-  GEOF53 

-  GEOPS4 

Static  Aerodynamic  Results  (I.OA!) 

-  INLTD 

-  SBODY 

-  SB1 

-  SB12 

-  SB123 

-  SB1234 

-  SFIN1 

-  SFIN2 

-  SFIN3 

-  SFIN4 

Dynamic  Aerodynamic  Results  (I.O.M.) 

-  DBODY 

-  DB1 

-  DB12 

-  DB123 

-  DB1234 

-  DDFIN1 

-  DDFIN2 

-  DDFIN3 

-  DDFIN4 

Aerodynamic  Trim  Results 

-  UTRIMD 

-  TRIMD 


Figure  59  Program  Common  Blocks  By  Data  Type 
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5.3  AERODYNAMIC  METHODOLOGY 

This  section  briefly  summarizes  the  method  routines  incorporated  in 
the  Missile  Datcom  code.  Also  covered  in  this  section  is  the  means  to  update 
or  replace  a  method. 

5.3.1  Methods  Incorporated 

The  methods  incorporated  are  summarized  in  Table  7.  Each  method  is 
coded  into  its  own  subroutine  so  that  revision  or  replacement  is  easily 
accomplished.  Detailed  documentation  within  the  code  using  "comment 
cards"  further  describes  the  methods  as  well  as  their  limitations. 

5.3.2  Changing  a  Method 

Replacing  a  component  buildup  method  is  easily  done.  Since  each 
method  is  coded  in  an  individual  subroutine,  simply  replacing  the  method 
subroutine  will  implement  the  new  technique.  A  few  of  the  methods  are 
complex  an  require  several  subroutines;  these  are  called  method  modules. 
Method  modules  substitution  is  more  complex  but  can  still  be  easily 
accomplished.  The  program  development  philosophy,  described  below,  will 
aid  in  method  revision. 

Code  Structure:  The  code  was  developed  using  top-down  design.  This 
development  scheme  was  implemented  by  coding  at  the  top-most 
control  logic  downward  to  integration  of  the  individual  method 
subroutines.  Hence,  the  upper  levels  of  the  code  structure  contain  the 
basic  logic  to  implement  the  component  buildup  methods.  The  lower 
levels  are  the  implemented  methods.  In  most  cases  the  control  logic 
requires  no  changes. 

Method  Coding  Style:  Most  Methods  are  implemented  in  a  single 
subroutine.  Their  inputs  and  outputs  are  passed  through  the 
subroutine  calling  sequence.  Any  method  routine  can  be  extracted  and 
used  in  another  code  without  modification.  In  some  cases  utility 
routines,  such  as  table  look-ups,  are  used;  they  must  also  be  extracted  if 
the  method  is  to  be  used  in  another  code. 

Each  subroutine  includes  a  brief  description  of  the  inputs  and  outputs, 
the  reference  documentation,  and  any  limitations  or  assumptions. 

When  any  routine  contains  data  tables,  two  subroutines  calls  are 
inserted  so  that  the  program  execution  sequence  can  be  "tracked".  The 
call  to  STRACE  at  the  start  of  the  subroutine  places  the  name  of  the 
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routine  being  executed  (parameter  IROUT)  into  the  common  block 
TRACE.  The  subroutine  name  is  removed  from  the  TRACE  common 
block  using  the  call  to  ETRACE  at  the  end  of  the  subroutine.  These 
routines  do  not  have  any  impact  on  other  calculations  of  the  program 
and  the  subroutine  calls  may  be  removed  if  the  code  is  used  in  another 
program. 

Execution  Sequence:  Subroutines  BODY  and  FINS  control  the 
calculation  sequence  for  body  alone  and  fin  alone,  respectively.  These 
two  routines  are  called  by  the  master  aerodynamic  calculation 
subroutine  AERO;  this  is  where  the  Mach  number  and  the  flight 
conditions  are  defined  for  the  user  input  case,  the  full  configuration 
component  build-up  is  done  in  subroutine  SYNTHS. 

The  aerodynamics  are  calculated  in  the  following  sequence:  1)  Normal 
force,  2)  Axial  force,  3)  Pitching  Moment,  4)  Side  Force,  5)  Yawing 
Moment,  6)  Rolling  Moment,  and  7)  the  derivatives  of  the  above  with 
respect  to  angle  of  attack  and  sideslip  angle.  In  some  cases,  the  calling 
sequence  must  not  be  changed  since  subsequent  results  are  dependant 
upon  other  coefficients.  For  example,  drag-to-lift  is  dependant  upon 
normal  force.  Extreme  caution  must  b.  exercised  when  revising  the 
method  execution  sequence.  It  is  recommended  that  the  same  example 
case  be  run  with  both  the  "old"  and  "new"  versions  of  the  code  and  any 
differences  be  reconciled. 

Special  options  of  the  code  such  as  experimental  data  substitution  and 
configuration  incrementing  depend  on  the  methods  by  which  the 
aerodynamic  coefficients  are  calculated.  Both  of  these  Options  are 
executed  in  the  subroutine  SYNTHS.  An  example  of  these  options 
dependence  on  the  computation  methods  is  the  separation  of  Cn  into 
CNo,  Cnp  and  Cnv-  Incrementing  factors  are  applied  to  each  of  these 
components  separately.  Therefore,  a  change  in  the  decomposition  of 
Cn  would  effect  the  configuration  incrementing  option.  Therefore  any 
change  in  the  method  of  computing  an  aerodynamic  coefficient  should 
be  checked  for  synthesis  ramifications. 

Changing  a  Method  Subroutine:  Revising  a  method  which  is  coded 
into  a  single  subroutine  is  as  simple  as  writing  a  routine  with  the  same 
name  and  substituting  it  into  the  program.  Any  changes  to  the 
variables  passed  through  the  routine  calling  sequence  must  also  be 
changed  in  those  routines  that  call  it.  Data  required  which  are  not 
available  in  the  calling  sequence  may  be  optionally  added  by  inserting 
the  appropriate  common  block  (see  Section  5.4).  Care  must  be  taken 
when  using  data  from  a  common  block  to  make  sure  that  it  has  been 
computed  prior  to  its  attempted  use. 
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Changing  a  Method  Module:  Four  methods  are  too  complex  to  be 
included  as  a  single  subroutine.  They  are  the  Airfoil  Section  Module, 
the  Hybrid  Theory  Module,  the  Second  Order  Shock  Expansion  Module 
and  the  Supersonic  Wing  Potential  Flow  Module.  These  techniques 
are  neither  short  nor  easily  changed.  To  replace  each  module  with 
another  technique,  the  following  is  recommended. 

•  Airfoil  Section  Module  -  this  module  starts  with  subroutine 
THEORY.  To  use  another  set  of  airfoil  section  calculations  requires 
the  revision  of  this  subroutine. 

•  Hybrid  Theory  Module  -  The  second-order  potential  flow  solution 
of  Van  Dyke  (Hybrid  Theory)  begins  with  subroutine  HYBRID. 
Replacement  of  this  method  requires  changes  to  subroutine 
SUPPOT. 

•  Second-Order  Shock  Expansion  Theory  Module  -  The  Second- 
Order  Shock  Expansion  method  is  implemented  beginning  with 
SOSE.  Replacement  of  this  method  requires  changes  to  subroutine 
SUPPOT. 

•  Supersonic  Wing  Potential  Flow  Module  -  The  potential  flow 
method  for  supersonic  wave  drag  is  implemented  in  subroutine 
FCAWPF.  Replacement  of  this  method  is  done  in  FINXCA. 

Format  of  PLOT  File:  When  the  PLOT  control  card  is  used,  a  formatted 
data  file  is  written  to  unit  3  which  can  be  used  in  a  separate  plotting 
program.  The  file  format  can  be  seen  in  Appendix  B. 
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Table  7  Summary  of  Methods  Implemented  in  Missile  Datcom 


Coef  f _ ROUTINE  DESCRIPTION  _  REFERENCES 

a  BaaaaagKsaB  ■  i  ■  ■=a=cs=  .n  saaaaa  aa=a=s==;ir,  .  ■  :  i 


aEQ 

ALPEQ  "" 

COMPUTES  EQUIVALENT  ANGLE  OF  ATTACK  OF  FIN 
PANEL 

J  OF  SPACECRAFT  i  ROCKETS, 
JULY- AUGUST  1983 

s* 

ANGDET 

WEDGE  TURN  ANGLE  FOR  WHICH  SHOCK  WILL  BECOME 
DETACHED 

NACA  1135 

EQNS.  138  l  168 

CAb 

BDCAB 

COMPUTES  BASE  DRAG  FOR  BODIES  AT  ALL  SPEEDS 

NASA  TR  R-100 

NSWC  TR-81-156  (NSWC  AERO 
CODE) 

CA,a, 

BDCALP 

COMPUTES  DRAG  DUE  TO  LIFT  FOR  AXISYMMETRIC 

BOD  I ES 

,VMCP  706-280 

CAp 

BDCAPR 

COMPUTES  SUBSONIC  PRESSURE  DRAG  FOR  BODIES 

ALLEN  AND  PERKINS 
(NACA  1048) 

(CAp) 

BT 

BDCAP1 

SUBSONIC-TRANSONIC  PRESSURE  DRAG  COEFFICIENT 
(NO  FRICTION)  OF  AXISYMMETRIC  CONICAL  BODIES 

DATCOM  SECTION  4. 2. 3.1 

<CAp) 

BT 

BDCAP2 

SUBSONIC-TRANSONIC  PRESSURE  DRAG  COEFFICIENT 
(NO  FRICTION)  OF  AXISYMMETRIC  OGIVAL  BODIES 

PAYNE-DTNSRDC/ASED-80/10 

MAY,  1980.  PP.  33-38 

CaP,w 

BDCAWC 

CALCULATES  PRESSURE/WAVE  DRAG  INCREMENT  OF 
CONICAL  NOSE  CYLINDER  AT  TRANSONIC  SPEEDS 

MOORE,  NSWC-TR-80-346, 

P.  28 

CAw 

BDCAWF 

SUBSONIC-TRANSONIC  WAVE  DRAG  COEFFICIENT  OF 
AXISYMMETRIC  CONICAL  BODIES  WITH  FLARE 

AMCP  706-280,  FIGURE  8-33A 

TO  8-33D, P.8-49 

CA 

BDCAWN 

COMPUTES  NOSE  WAVE  DRAG  AT  TRANSONIC  SPEEDS 
USING  CHAUSSEE,  UNSTEADY  EULER  SOLUTION 

NSWC  TR-80-346,  P.28 

c 

^mp 

BDCMP 

COMPUTES  POTENTIAL  PITCHING  MOMENT  FOR  BODIES 

ALLEN  AND  PERKINS 
(NACA  1048) 

^mv 

BDCMV 

COMPUTES  VISCOUS  PITCHING  MOMENT  FOR  BODIES 
ABOUT  USER  SPECIFIED  CENTER  OF  GRAVITY 

ALLEN  AND  PERKINS 
(NACA  1048) 

(CNa) 

BT 

BDCNAB 

COMPUTES  INCREMENTAL  NORMAL  FORCE  SLOPE  DUE 
TO  CONICAL  BOATTAIL  ON  AXISYMMETRIC  BODIES 

NSWC -TR-81-156,  P.110 

(cNa) 

fl 

BDCNAF 

COMPUTE  INCREMENTAL  NORAML  FORCE  SLOPE  DUE  TO 
FLARE  AT  SUBSONIC/TRANSONIC  SPEEDS 

AMCP  706-280,  JULY  1968 

CN 

Na 

BDCNAN 

INTERPOLATE  NORMAL  FORCE  COEFFICIENT  SLOPE 

FOR  CONE-CYLINDER  AND  OGIVE-CYLINDER  AT 

TRANSONIC  SPEEDS 

MBB  TN,  WE2-97 / 69 

MBB  TN,  WE12-88/70 

cn 

Np 

BDCNP 

COMPUTES  POTENTAIL  NORMAL  FORCE  FOR  BODIES 

ALLEN  AND  PERKINS 
(NACA  1048) 

Cm 

wv 

BDCNV 

COMPUTES  VISCOUS  NORMAL  FORCE  FOR  BODIES 

ALLEN  AND  PERKINS 

(NACA  1048) 

BT 

BDXCPB 

LOCATE  AFT  BODY  CENTER  OF  PRESSURE  (SLENDER 

BODY  THEORY) 

NSWC-TR-80-3  16, P  .  4 1 

FL 

BDXCPF 

COMPUTES  SUBSONIC/TRANSONIC  FLARE  CENTER  OF 
PRESSURE 

AMCP  706-280,  JULY  1968 

*cp 

BDXCPN 

INTERPOLATE  LONGITUDINAL  CENTER  OF  PRESSURE 

FOR  CONE-CYLINDER  AND  OGIVE-CYLINDER  AT 
TRANSONIC  SPEEDS 

MBB  TN,  WE2-97/69 

MBB  TN,  WEI 2-88 / 7 0 

Pu 

BFTAU 

DETERMINES  THE  REGION  OF  INFLUENCE  OF  PANEL 

IN  PROXIMITY  TO  ANOTHER 

ONR-CR215-226-4F, 

APPENDIX  D 

cA 

BODYCA 

LINEARLY  REDUCE  PRESSURE  DRAG  IN  MACH  RANGE 

1.0  -  1.2 

DATCOM  SECTION  4.  1.5.1 
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Table  7  Summary  of  Methods  Implemented  in  Missile  Datcom  (Continued) 


Coeff _ ROUTINE  DESCRIPTION _ REFERENCES 


CAf 

CAFRIC 

COMPUTE  FRICTION  COMPONENT  OF  AXIAL  FORCE  AT 
ANGLE  OF  ATTACK 

MDAC-WEST  AERODYNAMIC 
HANDBOOK  (M  8.08C-CD) 

DATCOM  SECTION  4.  1.5.1 

FLUID  DYNAMIC  DRAG (HOERNER) 

CARRYO 

COMPUTE  FIN-BODY  AND  BODY-FIN  CARRY-OVER 
FACTORS  DUE  TO  ANGLE  OF  ATTACK  AND  INCIDENCE 

NACA  1307 

J.  AIRCRAFT,  OCTOBER  1975 

CARRYS 

SUPERSONIC  FIN-BODY  CARRY  OVER  WITH  FINITE 
AFTERBODIES 

AIAA  JOURNAL: 

VOL. 19, NO. 5, MAY  1981, P.661 
VOL.20.NO. 6, JUN  1982, P.855 
VOL.20.NO. 8, AUG  1982, P.1144 

C*C 

CDCS 

COMPUTE  CROSS  FLOW  DRAG  COEFFICIENT 
(JORGENSEN)  DATA  FAIRING  GUIDE  (BAKER) 

NASA  TN  D-6996,  FIGS. 1,2, 3 
AEDC-TR-75-124 

(Sdegi  v 

CLVR 

COMPUTE  PANEL  LIFT  INCREMENT  DUE  TO  BODY 
VORTICES  AS  A  FRACTION  OF  THE  THEORETICAL 
PANEL  LIFT 

NWC  TP-5761 

Cp 

CONEP 

COMPUTE  CONE  PRESSURE  AT  ZERO  ANGLE  OF  ATTACK 
(RASMUSSEN) 

AIAA  JOURNAL,  AUG  1967, 
P.1495 

Cp<a,0) 

CPDIST 

COMPUTE  PRESSURES  AROUND  BODY  (DEJARNETTE 
EQNS.28) 

AIAA  JOURNAL  OF  SPACECRAFT, 
NOV-DSC  1980  P.529 

“2 _ 

EQN4  4 

COMPUTE  MACH  NUMBER  (FROM  EQUATION  44) 

NACA  1135 

X 

n 

V 

FALCP 

CALCULATE  FIN  ALONE  AERODYNAMIC  CENTER  FOR 

AFT  SWEPT  WINGS 

DATCOM  SECTION  4.  1.4. 2, 
FIGURES  26(A)  -  26(F) 

CaLE 

FCALE 

COMPUTE  DRAG  INCREMENT  DUE  TO  LEADING  LEADING 
EDGE  BLUNTNESS 

DATCOM  SECTION  4. 1.5. I 

Cft  (a) 

FCALP 

COMPUTE  DRAG  DUE  TO  LIFT  FOR  FINS  ALONE 

DATCOM  SECTION  4.  1.5. 2 

CA« 

FCAWPF 

COMPUTE  FIN  ALONE  WAVE  DRAG  AT  SUPERSONIC 
SPEEDS 

NSWC  TR-80-346 

C*v 

FCAWS 

SUPERSONIC  FIN  WAVE  DRAG  INCREMENT  TO  AXIAL 
FORCE 

DATCOM  SECTION  4 . 1 . 5 . 1 

C*p,w 

FCAWT 

COMPUTE  FIN  ALONE  TRANSONIC  WAVE  DRAG 

INCREMENT 

DATCOM  SECTION  4.  1.5.1 

ClWAX 

FCLMSB 

COMPUTE  SUBSONIC  MAXIMUM  LIFT  AND  ANGLE  OF 
ATTACK  FOR  MAXIMUM  LIFT  FOR  LOW  ASPECT  RATIO 
WINGS 

DATCOM  SECTION  4. 1.3. 4 

METHOD  3 

%«* 

Clmax 

FCLMSH 

COMPUTE  SUBSONIC  MAXIMUM  LIFT  AND  ANGLE  OF 

ATTACK  FOR  MAXIMUM  LIFT  FOR  HIGH  ASPECT  RATIO 
WINGS 

DATCOM  SECTION  4. 1.3. 4 

METHOD  2 

CNa 

FCNA 

COMPUTES  FIN  ALONE  LINEAR  NORMAL  FORCE  CURVE 
SLOPE.  PER  DEGREE 

DATCOM  SECTION  4. 1.3. 2 

RAS  DATA  SHEETS 

Cflaa 

FCNAAH 

COMPUTES  SUPERSONIC  NON-LINEAR  NORMAL  FORCE 

PER  SIN2a 

DATCOM  SECTION  4.  1.3. 3 

CNaa 

FCNAAS 

COMPUTES  SUBSONIC  NON-LINEAR  NORMAL  FORCE  PER 

siN2a 

DATCOM  SECTION  4.  1.3.  3 

°Naa 

FCNAAT 

COMPUTES  TRANSONIC  NON-LINEAR  NORMAL  FORCE 

PER  SIN2a 

DATCOM  SECTION  4  .  1 .  3 . j 

Cfla 

FCNASB 

COMPUTES  FIN-ALONE  LIFT  CURVE  SLOPE  FOR  TWO 
(2)  PANELS  CONNECTED  AT  THE  ROOT  CHORD 
(LOWRY-POLHAMUS) 

DATCOM  FIGURE  4  .  1  .  3 . 7 

CNa 

FCNASP 

COMPUTES  SUPERSONIC  LIFTING  SURFACE  CN  PER 
ALPHA 

DATCOM  FIGURE  4.1.3.2-56 
DATCOM  FIGURE  4.1,3.7-60 

DAC  SM  13110 

CNa 

FCNATR 

COMPUTES  FIN  ALONE  TRANSONIC  LIFT  CURVE  SLOPE 
FOR  TWO ( 2 )  PANELS  (EXPOSED  PLANFORM) 

BRITISH  DATA  SHEETS: 

S.  01 .03. 06,  S. 01 . 03. 05 

S. 01. 03. 04,  S. 02. 03. 03 

S. 08. 01. 02 
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Table  7  Summary  of  Methods  Implemented  in  Missile  Datcom  (Continued) 


Coeff  ROUTINE  DESCRIPTION _ REFERENCES _ 

"cN‘aa_"  ==^^G60B  "  '  INTERPOLATE  FOP  FTN  CNAA  I  DATCOM  FIGU»E  DATi-om^^^On"  4  .  1  .  3 

_ 4.1.3. 3-60B) _ 

CAb  FINCAB  COMPUTE  BASE  DRAG  OF  FTN  TRAILING  EDGE  NWC  TR-2796  r  ~ 


EMPIRICAL 

0.0 

<  M 

< 

3.1 

30 

0.0 

<  M 

< 

1.2 

2D 

1.2 

<  M 

< 

3.1 

1/M2 

M  > 

4.0 

o 

> 

V 

FINCAP 

FIN  ALONE  SUBSONIC  PRESSURE  DRAG 

FLUID  DYNAMIC  DRAG 

(HOERNER) 

(AD EG) 5 

FINFIN 

DETERMINE  DELTA-ALP HA-EQUIVALENT  OF  DEFLECTED 
FIN  IN  PRESENCE  OF  ADJACENT  CRUCIFORM  FINS 

NEAR  TR  125 

e* 

FSDETA 

CALCULATE  SHOCK  DETACHMENT  ANGLE-OF-ATTACK  ON 
FINS 

DATCOM  SECTION  4. 1.3. 3, 

PAGE  4.1.3.3-33  STEP  2B 

NACA  1135  EQN  138  t  168 

Xcp 

FWDXAC 

CALCULATE  FIN  ALONE  AERODYNAMIC  CENTER  FOR 
FORWARD  SWEPT  FINS 

AFWAL  TR-84-3084 

GETETA 

COMPUTE  SUBSONIC  CROSS  FLOW  DRAG 
PROPORTIONALITY  FACTOR 

AEDC-TR-75-124 

HYPER 

COMPUTES  PRESSURES  BY  MODIFIED  NEWTONIAN 

THEORY 

NASA  TND-176 

KWBALP 

DETERMINE  RATIO  K-W (B) /K-W (B) SBT  DUE  TO  ANGLE 
OF  ATTACK 

EMPIRICAL  CORRELATION  OF 

TEST  DATA  (TEMPORARILY  SET 

TO  1.0) 

LOADF 

PRESSURE  LOADING  FUNCTIONS  (DEJARNETTE) 

AIAA  JOURNAL  CF  SPACECRAFT, 
NOV-DEC  1980  P.529 

CN' 

PANLCN 

COMPUTES  TOTAL  FIN  NORMAL  FORCE  OF  A 

CRUCIFORM  FIN  SET  USING  THE  CONCEPT  OF 
EQUIVALENT  ANGLE  OF  ATTACK 

AIAA  PAPER  T>-il53 
(NIELSON) 

Er£ 

S-R 

SFWRW 

COMPUTES  SURFACE  VORTEX  LATERAL  POSITION  FOR 
LIFTING  SURFACE/BODY  VORTEX  INTERFERENCE 

NACA  TR-1307 

(PITTS,  NIELSON,  KAATARI ) 

es 

SHOCK 

GET  SHOCK  SHAPE  ANGLE  FROM  EQN.42  OF 

DEJARNETTE 

AIAA  JOURNAL  OF  SPACECRAFT, 
NOV-DEC  1980  P.529 

CF 

SKINF 

COMPUTES  SKIN  FRICTION  DRAG  USING  VAN  DRIEST 
METHOD  II  FOR  TURBULENT  FLOW  AND  BLAIUS  FOR 
LAMINAR  FLOW.  EMPIRICAL  CORRELATIONS  ARE 

USED  FOR  TRANSITION  AND  SURFACE  ROUGHNESS 

CALCULATIONS 

MDAC-WEST  AERODYNAMIC 

HANDBOOK  (M  8.080-CD) 

DATCOM  SECTION  4. 1.5.1 

FLUID  DYNAMIC  DRAG ( HOERNER) 

SLOPE 

COMPUTES  SUBSONIC  AIRFOIL  SECTION  LIFT  CURVE 
SLOPE,  AERODYNAMIC  CENTER,  AND  CRITICAL  MACH 
NUMBER 

AFFDL-TR-7 1-87 

SOSE 

PRESSURES  BY  MODIFIED  NEWTONIAN  AND  IMPROVED 

2ND  ORDER  SHOCK  EXPANSION  OF  DEJARNETTE 

AIAA  JOURNAL  OF  SPACECRAFT, 
NOV-DF.C  1980  P.529  (NSWC 

AERO  CODE) 

Iv 

SURINT 

COMPUTES  INTERFERENCE  FACTOR  FOR  VORTEX 
INTERACTION  WITH  LIFTING  SURFACE 

NACA  TR-1307,  APPENDIX  5 
(PITTS,  NIELSON,  KAATARI) 

SVTRAK 

COMPUTES  WING  VORTEX  HORIZONTAL  AND  VERTICAL 
LOCATION  AT  CENTER  OF  PRESSURE 

NACA  TR-1307 

(PITTS,  NIELSON,  KAATARI' 

VANDYK 

COMPUTES  SECOND  ORDER  AXIAL  AND  FIRST  ORDER 
CROSS  FLOW  PERTURBATION  VELOCITY  COMPONENTS 

NSWC  AERO  CODE 

Iv 

VRINTS 

COMPUTES  INTERFERENCE  FACTOR  FOR  VORTEX 
INTERACTION  WITH  LIFTING  SURFACE 

NACA  TR-1307,  APPENDIX  B 
(PITTS,  NIELSON,  KAATARI) 

Yep 

YCP 

SUBROUTINE  TO  COMPUTE  LATERAL  CENTER  OF 
PRESSURE 

AIAA  PAPER  91-0708 
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5.4  COMMON  BLOCK  DEFINITIONS 


This  section  summarizes  the  contents  of  each  major  common  block 
within  the  program.  If  the  block  of  data  is  input  using  namelist,  the  namelist 
name  is  noted.  If  the  array  can  be  dumped  using  the  DUMP  control  card  or 
written  using  the  WRITE  control  card,  the  name  is  shown  in  the  proper 
blank. 


All  data  is  stored  in  the  foot-pound-second-degree-Rankine  system  of 
units  regardless  of  the  dimensional  units  set  by  the  user.  A  data  element  is 
determined  "unused"  when  the  value  of  the  number  is  equal  to  1  x  10*30.  (See 
section  5.1.2  for  a  complete  definition  of  UNUSED) 

The  blocks  of  data  are  listed  alphabetically  in  this  section  according  to 
their  COMMON  BLOCK  names. 
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Table  8  AXIBOD  Namelist  Inputs  (Common  Block  ABODIN) 


DEFINITION  OF  COMMON  BLOCK  ABODIN 
(DUMP  NAME  BDIN  ,  WRITE  NAME  ABODIN) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

NX 

Number  of  body  stations 

- 

2 

XO 

Xo 

Body  apex  station 

ft 

3-52 

X 

X 

Longitudinal  stations 

ft 

53-102 

R 

r 

Radii 

ft 

103 

TNOSE 

Nose  shape  type 

0.  =  Conical  or  Cone 

1.  =  Tangent  ogive 

2.  =  Power  series 

3.  =  Haack 

4.  =  Von  Karam 

104 

LNOSE 

L 

Actual  nose  length 

ft 

105 

DNOSE 

D 

Nose  base  diameter 

ft 

106 

BNOSE 

b 

Nose  bluntness  radius 

ft 

107 

TRUNC 

.TRUE,  if  truncated 

- 

108 

LCENTR 

L 

Centerbody  length 

ft 

109 

DCENTR 

D 

Centerbody  base  diameter 

ft 

110 

TAFT 

Afterbody  shape  type 

0.  =  Conical  or  Cone 

1.  =  Ogive 

111 

LAFT 

L 

Afterbody  length 

ft 

112 

DAFT 

D 

Afterbody  base  diameter 

ft 

113 

POWER 

n 

Exponent  for  power  series  nose 

- 

114-133 

DISCON 

Indices  of  X  stations  where  surface  slope  is 

discontinuous 

- 

134-183 

ELLIP 

e 

Ellipticity  of  body  at  each  X  station 

- 

184-233 

H 

h 

Height  of  body  at  each  X  station 

ft 

234 

ENOSE 

Cn 

Nose  base  ellipticity 

- 

235 

ECENTR 

ec 

Centerbody  base  ellipticity 

- 

236 

EAFT 

e. 

Afterbody  base  ellipticity 

- 

237 

DEXIT 

^tah  ! 

Nozzle  exit  diameter  at  base 

ft 

NOTE:  Table  continues  on  next  page. 
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Table  8  A XI BOD  Namelist  Inputs  (Common  Block  ABODIN)  -  Continued 


DEFINITION  OF  COMMON  BLOCK  ABODIN 
(DUMP  NAME  BDIN  ,  WRITE  NAME  ABODIN) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

238 

BASE 

.TRUE,  if  jet  interaction  calculated 

- 

239 

BETAN 

Nozzle  exit  angle 

deg 

240-259 

JMACH 

Jet  Mach  number  at  nozzle  exit 

- 

260-279 

PRAT 

Jet  to  freestream  static  pressure  ratio 

- 

280-299 

TRAT 

Jet  to  freestream  stagnation  temperature  ratio 

- 

300 

PROTUB 

.TRUE,  if  protuberance  drag  is  calculated 

- 

301 

NPROT 

Number  of  protuberance  sets 

- 

302-321 

PTYPE 

Protuberance  set  type: 

1.  =  Vertical  cylinder 

2.  =  Horizontal  cylinder 

3.  =  Launch  lug 

4.  =  Launch  shoe 

5.  =  Block 

6.  =  Fairing 

7.  =  Component  build  (Not  Used) 

322-341 

XPROT 

Longitudinal  distance  from  missile  nose  to 
protuberance  set 

ft 

342-361 

NLOC 

Number  of  protuberances  in  each  protuberance 

- 

set 

362-381 

BLDMEM 

Number  of  protuberance  types  in  component 

build-up 

382-481 

BLDTYP 

Types  of  protuberances  in  component  build-up 

- 

482-581 

LPROT 

Length  of  each  member  or  protuberance 

ft 

582-681 

WPROT 

Width  of  each  member  or  protuberance 

ft 

682-781 

HPROT 

Height  of  each  member  or  protuberance 

ft 

782-881 

OPROT 

Vertical  offset  of  each  member  of  protuberance 

ft 
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Table  9  ELLBOD  Namelist  Inputs  (Common  Block  ABODIN) 


DEFINITION  OF  COMMON  BLOCK  ABODIN 
(DUMP  NAME  BDIN  ,  WRITE  NAME  ABODIN) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

NX 

Number  of  body  stations 

- 

2 

xo 

Xo 

Body  apex  station 

ft 

3-52 

X 

X 

Longitudinal  stations 

ft 

53-102 

R 

r 

Radii 

ft 

103 

TNOSE 

Nose  shape  type 

0.  =  Conical  or  Cone 

1.  =  Tangent  ogive 

2.  =  Power  series 

3.  =  Haack 

4.  =  Von  Karam 

104 

LNOSE 

L 

Actual  nose  length 

ft 

105 

DNOSE 

D 

Nose  base  diameter 

ft 

106 

BNOSE 

b 

Nose  bluntness  radius 

ft 

107 

TRUNC 

.TRUE,  if  truncated 

- 

108 

LCENTR 

L 

Centerbody  length 

ft 

109 

DCENTR 

D 

Centerbody  base  diameter 

ft 

110 

TAFT 

Afterbody  shape  type 

0.  =  Conical  or  Cone 

1.  =  Ogive 

111 

LAFT 

L 

Afterbody  length 

ft 

112 

DAFT 

D 

Afterbody  base  diameter 

ft 

113 

POWER 

n 

Exponent  for  power  series  nose 

- 

114-133 

DISCON 

Indices  of  X  stations  where  surface  slope  is 

discontinuous 

" 

134-183 

ELL  IP 

e 

Ellipticity  of  body  at  each  X  station 

- 

184-233 

H 

h 

Height  of  body  at  each  X  station 

ft 

234 

ENOSE 

<5n 

Nose  base  ellipticity 

- 

235 

ECENTR 

Centerbody  base  ellipticity 

- 

236 

EAFT 

e. 

Afterbody  base  ellipticity 

- 

237 

DEXIT 

_ _ 

Nozzle  exit  diameter  at  base 

ft 

NOTE:  Table  continues  on  next  page. 
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Table  9  ELLBOD  Namelist  Inputs  (Common  Block  ABODIN)  -  Continued 


DEFINITION  OF  COMMON  BLOCK  ABODIN 
(DUMP  NAME  BDIN  ,  WRITE  NAME  ABODIN) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

238 

BASE 

.TRUE,  if  jet  interaction  calculated 

- 

239 

BETAN 

Nozzle  exit  angle 

deg 

240-259 

JMACH 

Jet  Mach  number  at  nozzle  exit 

- 

260-279 

PRAT 

Jet  to  freestream  static  pressure  ratio 

- 

280-299 

TRAT 

Jet  to  freestream  stagnation  temperature  ratio 

- 

300 

PROTUB 

.TRUE,  if  protuberance  drag  is  calculated 

- 

301 

NPROT 

Number  of  protuberance  sets 

- 

302-321 

PTYPE 

Protuberance  set  type: 

1.  =  Vertical  cylinder 

2.  =  Horizontal  cylinder 

3.  =  Launch  lug 

4.  =  Launch  shoe 

5.  =  Block 

6.  =  Fairing 

7.  =  Component  build  (Not  Used) 

322-341 

XPROT 

Longitudinal  distance  from  missile  nose  to 

protuberance  set 

ft 

342-361 

NLOC 

Number  of  protuberances  in  each  protuberance 

- 

set 

362-381 

BLDMEM 

Number  of  protuberance  types  in  component 
build-up 

382-481 

BLDTYP 

Types  of  protuberances  in  component  build-up 

- 

482-581 

LPROT 

Length  of  each  member  or  protuberance 

ft 

582-681 

WPROT 

Width  of  each  member  or  protuberance 

ft 

682-781 

HPROT 

Height  of  each  member  or  protuberance 

ft 

782-381 

OPROT 

Vertical  offset  of  each  member  of  protuberance 

ft 
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Table  10  Body  Aerodynamic  Work  Array  (Common  Block  BDWOKK) 


DEFINITION  OF  COMMON  BLOCK  BDWORK 
(DUMP  NAME  BDWK  ,  WRITE  NAME  BDWORK) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

■ 

CNP 

CNp 

Potential  nonnal  force  vs.  a 

- 

CMP 

Cmp 

Potential  pitching  moment  vs.  a 

- 

41-60 

CNVIS 

Cnv 

Viscous  normal  force  vs.  a 

- 

61-80 

CMVIS 

Qny 

Viscous  pitching  moment  vs.  a 

- 

81-100 

CAPR 

CAp,W 

Pressure/wave  axial  force  vs.  a 

- 

101-120 

CAF 

CAf 

Friction  axial  force  vs.  a 

- 

121-140 

CABASE 

CAb 

Base  axial  force  vs.  a 

- 

141 

ETA 

Cross-flow  proportionality  factor 

- 

142-161 

CDC 

Qc 

Cross-flow  drag  coefficient  vs.  a 

- 

162-181 

CAPROT 

Protuberance  axial  force  coefficient 

_ 

182-201 

BOTDCA 

Boattail  incremental  axial  force  due  to  separation 

202-221 

BOTDCN 

Boattail  incremental  nonnal  force  due  to 
separation 

- 

222-241 

BOIDCM 

Boattail  incremental  pitching  moment  due  to 
separation 
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Table  11  Case  Identification  (Common  Block  CASEID) 


DEFINITION  OF  COMMON  BLOCK  CASEID 

(DUMP  NAME 

.WRITE NAME  CASEID  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1-74 

'5 

'6-175 

176 

177 

178 

179 

180 


IDCASE 

KOUNT 

NAMSV 

CASE 

NOEXTR 

NOLAT 

IR 

IPAGE 


ase  I.D.,  one  character  per  element 

Number  of  saved  namelists 

Saved  namelist  order  (packed  3  per  element) 

Case  number 

.TRUE,  if,  no  extrapolation  messages 
.  TRUE  .  if  no  lat-dir  derivatives  to  be  computed 

Run  number  for  plot  file 
Page  number  of  output 


Table  12  Program  Constants  (Common  Block  CONST) 


DEFINITION  OF  COMMON  BLOCK  CONST 

(DUMP  NAME 

.WRITE NAME  CONST  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

Value  for* 

Value  for  180/rc 

Value  for  unused  (1  x  10"^) 

Namelist  delimineter  ($) 
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Table  13  Dynamic  Derivatives  for  Body  and  Finset  1  (Common  Block  DB1) 


DEFINITION  OF  COMMON  BLOCK  DB1 
(DUMP NAME  DB1  .WRITE NAME  DB1  ) 


UNITS 


Normal  force  due  to  pitch  vs.  a 

Normal  force  due  to  time  rate  of  change  of  angle 

of  attack  vs.  a 

Total  damping  in  pitch  vs.  a 
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Table  14  Dynamic  Derivatives  for  Body  and  Finset  1,2  (Common  Block  DB12) 


DEFINITION  OF  COMMON  BLOCK  DB12 
(DUMP NAME  DB12  .WRITE NAME  DB12  ) 


VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1-20 

T2CNQ 

CNq 

Normal  force  due  to  pitch  vs.  a 

1/rad 

21-40 

T2CNAD 

CN- 

Normal  force  due  to  time  rate  of  change  of  angle 

1/rad 

of  attack  vs.  a 

41-60 

T2PTCH 

Cm- 

Total  damping  in  pitch  vs.  a 

1/rad 

Table  15 


amic  Derivatives  for  Body  and  Finset  1,2,3  (Common  Block  DB123) 


DEFINITION  OF  COMMON  BLOCK  DB123 
(DUMP NAME  DB13  .WRITE NAME  DB123  ) 


LOCATION 


e 


41-60 


VARIABLE 

NAME 


T3CNQ 

T3CNAD 


T3PTCH 


ENGINEERING 

SYMBOL 


UNITS 


Normal  force  due  to  pitch  vs.  a 

Normal  force  due  to  time  rate  of  change  of  angle 

of  attack  vs.  a 

Total  damping  in  pitch  vs.  a 


Table  16  Dynamic  Derivatives  for  Body  and  Finset  1 ,2,3,4  (Common  Block  DB1234) 


DEFINITION  OF  COMMON  BLOCK  DB1234 
(DUMP NAME  DBM  .WRITE NAME  DB1234  ) 


LOCATION 


21-40 


VARIABLE 

NAME 


T4CNQ 

T4CNAD 


ENGINEERING 

SYMBOL 


UNITS 


Normal  force  due  to  pitch  vs.  a 

Normal  force  due  to  time  rate  of  change  of  angle 

of  attack  vs.  a 


41-60 


T4PTCH 


Total  damping  in  pitch  vs.  a 


Table  17  Bod 


amic  Derivatives  (Common  Block  DBODY) 


DEFINITION  OF  COMMON  BLOCK  DBODY 
(DUMP NAME  DBOD  .WRITE NAME  DBODY  ) 


LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

1-20 

3CNQ 

CNq 

21-40 

BCNAD 

CNa 

41-60 

BPTCH 

Cm- 

61 

CLIP 

C»P 

62 

CYP 

CYp 

63 

CNP1 

64 

CNP3 

65 

CNP5 

66 

CNPY5 

Normal  force  due  to  pitch  vs.  a 

Normal  force  due  to  time  rate  of  change  of 

angle  of  attack  vs.  a 

Total  damping  in  pitch  vs.  a 

Rolling  moment  due  to  roll  rate 

Side  force  due  to  roll  rate 

1st  order  yawing  moment  (magnus)  coefficient 

3rd  order  yawing  moment  (magnus)  coefficient 

5th  order  yawing  moment  (magnus)  coefficient 

5  degree  secant  slope  of  magnus  moment 
coefficient  (@  5  degrees  yaw) 


UNITS 


1/rad 

1/rad 

1/rad 

l/sin(a) 

l/sin3(a) 

l/sin5(a) 

l/sin(a) 
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Table  18  Finset  1  Dynamic  Derivatives  (Common  Block  DDFIN1) 


DEFINITION  OF  COMMON  BLOCK  DDF  INI 
(DUMP NAME  DF1  .WRITE NAME  DFIN1  ) 


VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1-20 

F1CNQ 

CNn 

Normal  force  due  to  pitch  vs.  a 

1/rad 

21-40 

F1CNAD 

CNa 

Normal  force  due  to  time  rate  of  change  of  angle 

1/rad 

of  attack  vs.  a 

41-60 

F1CMQ 

cmq 

Pitching  moment  due  to  pitch  rate  vs.  a 

1/rad 

61-80 

F1CMAD 

Pitching  moment  due  to  time  rate  of  change  of 

1/rad 

angle  of  attack  vs.  a 

Table  20  Finset  3  Dynamic  Derivatives  (Common  Block  DDFIN3) 


DEFINITION  OF  COMMON  BLOCK  DDFIN3 

(DUMP  NAME  DF3 

.WRITE NAME  DFIN3  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

Table  21  Finset  4  Dynamic  Derivatives  (Common  Block  DDFIN4) 


DEFINITION  OF  COMMON  BLOCK  DDF1N4 
(DUMP NAME  DF4  .WRITE NAME  DFIN4 


DEFINITION 


UNITS 


Normal  force  due  to  pitch  vs.  a 

Normal  force  due  to  time  rate  of  change  of  angle 

of  attack  vs.  a 

Pitching  moment  due  to  pitch  rate  vs.  a 
Pitching  moment  due  to  time  rate  of  change  of 
angle  of  attack  vs.  a 


Table  22  NACA  Designation  (Common  Block  DESIG) 


DEFINITION  OF  COMMON  BLOCK  DESIG 
(DUMP NAME  .WRITE NAME  DESIG  ) 


LOCATION 


1-320 


VARIABLE 

ENGINEERING 

NAME 

SYMBOL 

DEFINITION 

NNACA 


NACA  designation  by  fin  set  (80  each) 


Table  23  Delete  Flags  for  Input  NAMELISTs  (Common  Block  DFLAGS) 


DEFINITION  OF  COMMON  BLOCK  DFLAGS 
(DUMP NAME  .WRITE NAME  DFLAGS  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

DFLT 

Delete  prior  SFLTCON  flag 

- 

2 

DREF 

Delete  prior  $REFQ  flag 

- 

3 

DAXI 

Delete  prior  SAXIBOD  flag 

- 

4 

DFIN1 

Delete  prior  SFTNSETl  flag 

- 

5 

DFTN2 

Delete  prior  SFINSET2  flag 

- 

6 

DFIN3 

Delete  prior  $FTNSET3  flag 

- 

7 

DF1N4 

Delete  prior  $FTNSET4  flag 

- 

8 

DDEFL 

Delete  prior  SDEFLCT  flag 

- 

9 

DTRIM 

Delete  prior  $TRIM  flag 

- 

10 

DELLB 

Delete  prior  $ELLBOD  flag 

- 

11 

DINLET 

Delete  prior  $  INLET  flag 

- 

12 

DARBOD 

Delete  prior  $ARBBOD  flag 
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Table  24  Pump  Array  Flags  (Common  Block  DUMPF) 


DEFINITION  OF  COMMON  BLOCK  DUMPF 
(DUMP NAME  .WRITE NAME  DUMPF  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

—i  it  a»] 

Flag  to  dump  GEOBOE) 

- 

2 

LF1GM 

Flag  to  dump  GEOFS1 

- 

3 

LF2GM 

Flag  to  dump  GEOFS2 

- 

4 

LF3GM 

Flag  to  dump  GEOFS3 

- 

5 

LF4GM 

Flag  to  dump  GEOFS4 

- 

6 

LATMP 

Flag  to  dump  FLC  and  TOTALC 

- 

7 

LBDWK 

Flag  to  dump  BDWORK 

- 

8 

LFLCT 

Flag  to  dump  FLC 

- 

9 

LINLD 

Flag  to  dump  INLTD 

- 

10 

LINPT 

Flag  to  dump  input  data 

- 

11 

LFLTC 

Flag  to  dump  FLCT 

- 

12 

LREFQN 

Flag  to  dump  REFQN 

- 

13 

LB  DIN 

Flag  to  dump  ABODIN  orEBODIN 

- 

14 

LF1IN 

Flag  to  dump  FSET1 

- 

15 

LF2IN 

Flag  to  dump  FSET2 

16 

LF3IN 

Flag  to  dump  FSET3 

17 

LF4LN 

Flag  to  dump  FSET4 

18 

LINLEN 

Flag  to  dump  INLETN 

19 

UOM 

Flag  to  dump  the  I.O.M. 

20 

LSBOD 

Flag  to  dump  SBODY 

21 

LSF1 

Flag  to  dump  SFIN1 

22 

LSF2 

Flag  to  dump  SFIN2 

23 

LSF3 

Flag  to  dump  SFIN3 

24 

LSF4 

Flag  to  dump  SFIN4 

- 

25 

LSB1 

Flag  to  dump  SB  1 

- 

26 

LSB12 

Flag  to  dump  SB  12 

- 

27 

LSB123 

Flag  to  dump  SB  123 

- 

28 

LS1234 

Flag  to  dump  SB1234 

- 

29 

LBOD 

Flag  to  dump  DBODY 

- 

30 

LDF1 

Flag  to  dump  DFIN1 

- 

31 

LDF2 

Flag  to  dump  DFIN2 

- 

32 

LDF3 

Flag  to  dump  DFIN3 

- 

33 

LDF4 

Flag  to  dump  DFIN4 

- 

34 

LDB1 

Flag  to  dump  DB1 

- 

35 

LDB12 

Flag  to  dump  DB12 

- 

36 

LDB123 

Flag  to  dump  DB 123 

- 

37 

LD1234 

Flag  to  dump  DB1234 
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Table  25  Finset  1  Aerodynamic  Work  Array  (Common  Block  FI  WORK) 


DEFINITION  OF  COMMON  BLOCK  F1WORK 
(DUMP  NAME  F1WK  .  WRITE  NAME  F1WORK) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

RHOl 

PLE 

Panel  effective  L.E.  radius 

ft 

2 

TMAX1 

WC>MAX 

Panel  effective  maximum  t/c 

- 

3 

KSHAR1 

k 

Wave  drag  parameter  for  section 

- 

4-23 

CCLA1 

Cl« 

Section  Ci^  vs.  Mach 

1/deg 

24-43 

XAC1 

X* 

Section  Xac  vs.  Mach 

- 

44-63 

CMC041 

cm cJ4 

Section  Cm  about  c/4 

- 

64 

CNALF1 

CNa 

Single  panel  CNa 

1/deg 

65-84 

CNAAF1 

^Naa 

Single  panel  Cn^  vs.  a 

1/nd2 

85-104 

CNLF1 

CNL 

Fin  set  linear  Cn 

- 

105-124 

CNNLF1 

°NnL 

Fin  set  non-linear  Cn 

- 

125-144 

CNF1AT 

Cn 

Fin  set  total  Cn 

- 

145 

XCPL1 

XCPL 

Single  panel  linear  C.P. 

ft 

146 

XCPNL1 

XCPNL 

Single  panel  non-linear  C.P. 

ft 

147-166 

CMFL1 

CmL 

Fin  set  linear  Cm 

- 

167-186 

CMFNL1 

CmNL 

Fin  set  non-linear  Cm 

- 

187-206 

CMF1AT 

Cm 

Fin  set  total  Cm 

- 

207 

CA0F1 

Cao 

Single  panel  Caq 

- 

208-227 

CANLF1 

canl 

Single  panel  (Ca-Caq)  vs  a 

- 

228-247 

ALPTF1 

aj 

Interpolated  a  for  panel  char. 

deg 

248-267 

CNFIT1 

CNj 

Interpolated  Cn  for  panel  char. 

- 

268 

AH 

Oideal 

Ideal  a  for  section 

deg 

269 

AL01 

OOL 

Zero  lift  a  for  section 

deg 

270 

CU1 

CLjdeal 

Ideal  CL  for  section 

- 

271-290 

CLM1 

^Lmax 

Maximum  airfoil  section  Cl 
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Table  26  Finset  2  Aerodynamic  Work  Array  (Common  Block  F2WORK) 


DEFINITION  OF  COMMON  BLOCK  F2WORK 
(DUMP  NAME  F2WK  ,  WRITE  NAME  F2WORK) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

RH02 

PLE 

Panel  effective  L.E.  radius 

ft 

2 

TMAX2 

(^max 

Panel  effective  maximum  t/c 

- 

3 

KSHAR2 

k 

Wave  drag  parameter  for  section 

- 

4-23 

CCLA2 

CJq 

Section  Cfo  vs.  Mach 

1/deg 

24-43 

XAC2 

Xac 

Section  Xac  vs.  Mach 

- 

44-63 

CMC042 

Cm^ 

Section  Cm  about  c/4 

- 

64 

CNALF2 

°Na 

Single  panel  CNa 

1/deg 

65-84 

CNAAF2 

CNaa 

Single  panel  Cn^  vs.  a 

1/rad2 

85-104 

CNLF2 

CNl 

Fin  set  linear  Cn 

- 

105-124 

CNNLF2 

°nnl 

Fin  set  non-linear  Cn 

- 

125-144 

CNF2AT 

Cn 

Fin  set  total  Cn 

- 

145 

XCPL2 

XCPl 

Single  panel  linear  C.P. 

ft 

146 

XCPNL2 

XCPnl 

Single  panel  non-linear  C.P. 

ft 

147-166 

CMFL2 

CmL 

Fin  set  linear  Cm 

- 

167-186 

CMFNL2 

CmNL 

Fin  set  non-linear  Cm 

- 

187-206 

CMF2AT 

Cm 

Fin  set  total  Cm 

- 

207 

CA0F2 

Cao 

Single  panel  Caq 

- 

208-227 

CANLF2 

CANL 

Single  panel  (Ca-Caq)  vs  a 

- 

228-247 

ALPTF2 

°j 

Interpolated  a  for  panel  char. 

deg 

248-267 

CNFIT2 

CNj 

Interpolated  Cn  for  panel  char. 

- 

268 

AI2 

oqdeal 

Ideal  a  for  section 

deg 

269 

AL02  | 

°OL 

Zero  lift  a  for  section 

deg 

270 

CLI2 

CLjdeal 

Ideal  CL  for  section 

- 

271-290 

CLM2 

^■Lmax 

Maximum  airfoil  section  Cl 
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Table  27  Finset  3  Aerodynamic  Work  Array  (Common  Block  F3WORK) 


DEFINITION  OF  COMMON  BLOCK  F3WORK 
(DUMP  NAME  F3WK  ,  WRITE  NAME  F3WORK) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

RH03 

PLE 

Panel  effective  L.E.  radius 

& 

2 

TMAX3 

^c)max 

Panel  effective  maximum  t/c 

- 

3 

KSHAR3 

k 

Wave  drag  parameter  for  section 

- 

4-23 

COLA3 

Cl« 

Section  C^  vs.  Mach 

1/deg 

24-43 

XAC3 

xac 

Section  Xac  vs.  Mach 

- 

44-63 

CMC043 

Cmc/4 

Section  Cm  about  c/4 

- 

64 

CNALF3 

°Na 

Single  panel  CNa 

1/deg 

65-84 

CNAAF3 

^Noux 

Single  panel  Cn^  vs.  a 

1/rad2 

85-104 

CNLF3 

Cnl 

Fin  set  linear  Cn 

- 

105-124 

CNNLF3 

°nnl 

Fin  set  non-linear  Cn 

- 

125-144 

CNF3AT 

Cn 

Fin  set  total  Cn 

- 

145 

XCPL3 

XCPl 

Single  panel  linear  C.P. 

ft 

146 

XCPNL3 

XCPNL 

Single  panel  non-linear  C.P. 

ft 

147-166 

CMFL3 

OnL 

Fin  set  linear  Cffl 

- 

167-186 

CMFNL3 

CmNL 

Fin  set  non-linear  Cm 

- 

187-206 

CMF3AT 

Cm 

Fin  set  total  Cm 

- 

207 

CA0F3 

Cao 

Single  panel  Cao 

- 

208-227 

CANLF3 

CAnl 

Single  panel  (Ca-Caq)  vs  a 

- 

228-247 

ALPTF3 

Interpolated  a  for  panel  char. 

deg 

248-267 

CNFIT3 

CNj 

Interpolated  Cn  for  panel  char. 

- 

268 

AI3 

oqdeal 

Ideal  a  for  section 

deg 

269 

AL03 

aCL 

Zero  lift  a  for  section 

deg 

270 

CU3 

CLjdeal 

Ideal  CL  for  section 

- 

271-290 

CLM3 

^-^Lfnax 

Maximum  airfoil  section  Cl 

180 


Table  28  Finset  4  Aerodynamic  Work  Array  (Common  Block  F4WORK) 


DEFINITION  OF  COMMON  BLOCK  F4WORK 
(DUMP  NAME  F4WK  ,  WRITE  NAME  F4WORK) 

LOCATION 

VARIABLE 

ENGINEERING 

DEFINITION 

UNITS 

NAME 

SYMBOL 

1 

RH04 

PLE 

Panel  effective  L.E.  radius 

ft 

2 

TMAX4 

^C^MAX 

Panel  effective  maximum  t/c 

- 

3 

KSHAR4 

k 

Wave  drag  parameter  for  section 

- 

4-23 

CCLA4 

Ciq 

Section  vs.  Mach 

1/deg 

24-43 

XAC4 

Xgc 

Section  Xac  vs.  Mach 

- 

44-63 

CMC044 

Section  Cm  about  c/4 

- 

64 

CNALF4 

°Na 

Single  panel  Cnq 

1/deg 

65-84 

CNAAF4 

cNaa 

Single  panel  Cn^  vs.  a 

1/rad2 

85-104 

CNLF4 

Cnl 

Fm  set  linear  Cn 

- 

105-124 

CNNLF4 

°nnl 

Fin  set  non-linear  Cn 

- 

125-144 

CNF4AT 

Cn 

Fm  set  total  Cn 

- 

145 

XCPL4 

XCPl 

Single  panel  linear  C.P. 

ft 

146 

XCPNL4 

XCPnl 

Single  panel  non-linear  C.P. 

ft 

147-166 

CMFL4 

cmL 

Fin  set  linear  Cm 

- 

167-186 

CMFNL4 

CfnNL 

Fm  set  non-linear  Cm 

- 

187-206 

CMF4AT 

Cm 

Fin  set  total  Cm 

- 

207 

CA0F4 

cAo 

Single  panel  Caq 

__ 

208-227 

CANLF4 

CAnl 

Single  panel  (C\-Caq\  vs  a 

- 

228-247 

ALPTF4 

Interpolated  a  for  panel  char. 

deg 

248-267 

CNFIT4 

CNj 

Interpolated  Cn  for  panel  char. 

- 

268 

AI3 

“ideal 

Ideal  a  for  section 

deg 

269 

AL04 

«OL 

Zero  lift  a  for  section 

deg 

270 

CLI4 

CLideal 

Ideal  CL  for  section 

- 

271-290 

CLM4 

^Lmax 

Maximum  airfoil  section  Cl 

Table  29  FLTCON  Namelist  Inputs  (Common  Block  FLC) 


DEFINITION  OF  COMMON  BLOCK  FLC 
(DUMP NAME  FLT  .WRITE NAME  FLC 


LOCATION 


25-44 


46-65 

66-85 

86-105 

106-125 


VARIABLE  ENGINEERING 
NAME  SYMBOL 


NALPHA 


ALPHA 

BETA 


NMACH 

MACH 


VINF 

TINF 

PINF 


Number  of  angles  of  attack 
Angles  of  attack 
Sideslip  angle 
Roll  angle 

Number  of  Mach  numbers 
Mach  numbers 
Geometric  altitude 
Reynolds  number 
Free-stream  velocity 
Free-stream  temperature 
Free-stream  pressure 


UNITS 


182 


Table  30  FINSET1  Namelist  Inputs  (Common  Block  FSET1) 


DEFINTHON  OF  COMMON  BLOCK  FSET1 

(DUMP  NAME  FUN 

.WRITE NAME  FSET1  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

SECTYP 

Type  of  airfoil  section 

1.  =  NACA 

- 

0.  =  Hexagonal/diamond 

2.  -  Circular  arc 

3.  =  User  defined 

2-11 

SSPAN 

W2 

Semi-span  stations 

ft 

12-21 

LMAXU 

(t/c)maxu 

L.E.  to  max  t/c,  upper  surface 

- 

22-31 

LFLATU 

lu 

Length  of  constant  t/c,  upper  surface 

- 

3141 

LMAXL 

(t/c)maxi 

L.E.  to  max,  t/c,  lower  surface 

- 

42-51 

LFLATL 

ll 

Length  of  constant  t/c,  lower  surface 

- 

52-61 

CHORD 

c 

Chord  length 

ft 

62-71 

THICKU 

0/c)u 

t/c  of  upper  section 

- 

72-81 

THICKL 

(t/c)i 

t/c  of  lower  section 

- 

82-91 

SWEEP 

A 

Sweep-back  angle 

deg 

92-101 

STA 

T1 

Station  for  measuring  sweep 

- 

102-111 

XLE 

Xle 

Station  for  L.E.  of  chord 

ft 

112 

NPANEL 

Number  of  panels 

- 

113-162 

XCORD 

X/c 

X/c  of  section 

- 

163-212 

MEAN 

Ymfc 

Y/c  of  section  (mean  line) 

- 

213-262 

THICK 

Vc 

t/c  of  section  (thickness  distribution) 

- 

263-312 

YUPPER 

Yu 

Y/c  of  upper  surface 

- 

313-362 

YLOWER 

Yl 

Y/c  or  lower  surface 

- 

363 

FINPHI 

4> 

Roll  angle  of  fins 

deg 

364-373 

LER 

rLE 

Fin  L.E.  radius 

ft 

374-381 

GAM 

r 

Fin  dihedral  angle 

deg 

382-389 

PHIF 

Fin  <j)  from  top  vertical  center 

deg 

390-399 

CFOC 

Cjfc 

Flap  chord  to  total  chord  ratio 

183 


Table  31  FINSET2  Namelist  Inputs  (Common  Block  FSET2) 


DEFINITION  OF  COMMON  BLOCK  FSET2 
(DUMP NAME  F2IN  .WRITE NAME  FSET2  ) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

l 

SECTYP 

Type  of  airfoil  section 

1.  =  NACA 

- 

2-11 

SSPAN 

Wl 

0.  =  Hexagonal/diamond 

2.  =  Circular  arc 

3.  =  User  defined 

Semi-span  stations 

ft 

12-21 

LMAXU 

(t/c)maxu 

L.E.  to  max  t/c,  upper  surface 

- 

22-31 

LFLATU 

lu 

Length  of  constant  t/c,  upper  surface 

- 

31-41 

LMAXL 

(t/c)maxi 

L.E.  to  max,  t/c,  lower  surface 

- 

42-51 

LFLATL 

II 

Length  of  constant  t/c,  lower  surface 

- 

52-61 

CHORD 

c 

Chord  length 

ft 

62-71 

THICKU 

(t/c)u 

t/c  of  upper  section 

- 

72-81 

THICKL 

(t/C)l 

t/c  of  Iowa*  section 

- 

82-91 

SWEEP 

A 

Sweep-back  angle 

deg 

92-101 

STA 

n 

Station  for  measuring  sweep 

- 

102-111 

XLE 

XLE 

Station  for  L.E.  of  chord 

ft 

112 

NPANEL 

Number  of  panels 

- 

113-162 

XCORD 

X/c 

X/c  of  section 

163-212 

MEAN 

Ym fc 

Y/c  of  section  (mean  line) 

- 

213-262 

THICK 

tfc 

t/c  of  section  (thickness  distribution) 

- 

263-312 

YUPPER 

Yu 

Y/c  of  upper  surface 

- 

313-362 

YLOWER 

Y| 

Y/c  or  lower  surface 

- 

363 

FINPHI 

Roll  angle  of  fins 

deg 

364-373 

LER 

rLE 

Fin  L.E.  radius 

ft 

374-381 

GAM 

r 

Fin  dihedral  angle 

deg 

382-389 

PHIF 

Fin  $  from  top  vertical  center 

deg 

390-399 

CFOC 

C (fc 

Flap  chord  to  total  chord  ratio 
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Table  32  FINSET3  Namelist  Inputs  (Common  Block  FSET3) 


DEFINITION  OF  COMMON  BLOCK  FSET3 

(DUMP  NAME  F3IN 

.WRITE NAME  FSET3  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

SECTYP 

Type  of  airfoil  section 

1.  =  NACA 

- 

0.  =  Hexagonal/diamond 

2.  =  Circular  arc 

3.  =  User  defined 

2-11 

SSPAN 

W2 

Semi-span  stations 

ft 

12-21 

LMAXU 

(t/c)maxu 

L.E.  to  max  t/c,  upper  surface 

- 

22-31 

LFLATU 

lu 

Length  of  constant  t/c,  upper  surface 

- 

31-41 

LMAXL 

(t/cJmaxi 

L.E.  to  max,  t/c,  lower  surface 

42-51 

LFLATL 

U 

Length  of  constant  t/c,  lower  surface 

- 

52-61 

CHORD 

c 

Chord  length 

ft 

62-71 

THICKU 

(l/c)u 

t/c  of ,  pper  section 

- 

72-81 

THICKL 

(t/c)l 

t/c  of  lower  section 

- 

82-91 

SWEEP 

A 

Sweep-back  angle 

deg 

92-101 

STA 

Station  for  measuring  sweep 

- 

102-111 

XLE 

Xle 

Station  for  L.E.  of  chord 

ft 

112 

NPANEL 

Number  of  panels 

- 

113-162 

XCORD 

X/c 

X/c  of  section 

- 

163-212 

MEAN 

Ymfc 

Y/c  of  section  (mean  line) 

• 

213-262 

THICK 

tfc 

t/c  of  section  (thickness  distribution) 

- 

263-312 

YUPPER 

YU 

Y/c  of  upper  surface 

- 

313-362 

YLOWER 

Yl 

Y/c  or  lower  surface 

- 

363 

FINPHI 

0 

Roll  angle  of  fins 

deg 

364-373 

LER 

rLE 

Fin  L.E.  radius 

ft 

374-381 

GAM 

r 

Fin  dihedral  angle 

deg 

382-389 

PHIF 

4>f 

Fin  4>  from  top  vertical  center 

deg 

390-399 

CFOC 

cpfc 

Hap  chord  to  total  chord  ratio 

185 


Table  33  FINSET4  Namelist  Inputs  (Common  Block  FSET4) 


DEFINITION  OF  COMMON  BLOCK  FSET4 

(DUMP  NAME  F41N 

.WRITE NAME  FSET4  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

SECTYP 

Type  of  airfoil  section 

1.  =  NACA 

■■ 

0.  =  Hexagonal/diamond 

2.  =  Circular  arc 

3.  =  User  defined 

2-11 

SSPAN 

m 

Semi-span  stations 

ft 

12-21 

LMAXU 

(t/c)maxu 

L.E.  to  max  t/c,  upper  surface 

- 

22-31 

LFLATU 

lu 

Length  of  constant  t/c,  upper  surface 

- 

31-41 

LMAXL 

(t/c)maxi 

L.E.  to  max,  t/c,  lower  surface 

. 

42-51 

LFLATL 

ll 

Length  of  constant  t/c,  lower  surface 

- 

52-61 

CHORD 

c 

Chord  length 

ft 

62-71 

THICKU 

(t/cXl 

t/c  of  upper  section 

- 

72-81 

THICKL 

mi 

t/c  of  lower  section 

- 

82-91 

SWEEP 

A 

Sweep-back  angle 

deg 

92-101 

STA 

Station  for  measuring  sweep 

- 

102-111 

XLE 

Xle 

Station  for  LJE.  of  chord 

ft 

112 

NPANEL 

Number  of  panels 

- 

113-162 

XCORD 

X/c 

X/c  of  section 

- 

163-212 

MEAN 

Ymfc 

Y/c  of  section  (mean  line) 

- 

213-262 

THICK 

t/c 

t/c  of  section  (thickness  distribution) 

- 

263-312 

YUPPER 

Yu 

Y/c  of  upper  surface 

- 

313-362 

YLOWER 

Yi 

Y/c  or  lower  surface 

- 

363 

FINPHI 

$ 

Roll  angle  of  fins 

deg 

364-373 

LER 

rLE 

Fin  L.E.  radius 

ft 

374-381 

GAM 

r 

Fin  dihedral  angle 

deg 

382-389 

PHIF 

Fin  4>  from  top  vertical  center 

deg 

390-399 

CFOC 

Cffc 

Flap  chord  to  total  chord  ratio 

186 


Table  34  Body  Geometry  Data  (Common  Block  GEOBOD) 


DEFINITION  OF  COMMON  BLOCK  GEOBOD 
(DUMP  NAME  GEOB  .WRITE  NAME  GEOBOD  ) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

2 

SpNose 

SWNOSE 

(Sp>N 

(Sw)jsi 

Nose  section  planform  area 

Nose  section  wetted  area 

ft2 

ft2 

3 

VOLNOS 

Vn 

Nose  section  volume 

ft3 

4 

XCVNOS 

(Xc.vJn 

Nose  section  volume  centroid 

ft 

5 

XCPNOS 

(Xc.P.)n 

Nose  section  planform  centroid 

ft 

6 

FRNOSE 

(Fr)n 

Nose  Section  fineness  ratio 

- 

7 

THEOLN 

1N* 

Nose  section  theoretical  length 

ft 

8 

THEOFN 

Frn* 

Nose  section  theoretical  fineness  ratio 

- 

9-16 

17-24 

25 

SPLAN 

Sp 

Repeat  of  1-  8  for  centerbody 

Repeat  of  1  -8  for  aft  body 

Total  planform  area 

ft2 

26 

SWET 

Sw 

Total  wetted  area 

ft2 

27 

VOL 

V 

Total  volume 

ft3 

28 

XCENTV 

Xc.V. 

Total  volume  centroid 

ft 

29 

XCENTP 

Xc.P. 

Total  planform  centroid 

ft 

30 

FR 

Fr 

Total  fineness  ratio 

- 

31 

DBASE 

db 

Base  diameter 

ft 

32 

SBASE 

Sb 

Base  area 

ft2 

33 

DMAX 

dmax 

Maximum  diameter 

ft 

34 

SMAX 

Smax 

Maximum  cross-section  area 

ft2 

35 

LTOTAL 

1T 

Total  body  length 

ft 

36 

BPNOSE 

(Spv>N 

Vertical  planform  area  of  nose 

ft2 

37 

BCPNOS 

(XC.P.v)n 

Vertical  planform  centroid  of  nose 

ft 

38 

BPCENT 

(SPv>C 

Vci  deal  planform  area  of  center  body 

ft2 

39 

BCPCEN 

(XC.P.y)c 

Vertical  planform  centroid  of  center  body 

ft 

40 

BPAFT 

(Spv)A 

Vertical  planform  area  of  aft  body 

ft2 

41 

BCPAFT 

(XcP.v)A 

Vertical  planform  centroid  of  aft  body 

ft 

42 

BPLAN 

Spy 

Tot?1  vertical  planform  area 

ft2 

43 

BCENTP 

Xc.P.y 

Total  vertical  planform  centroid 

ft 

44 

ECSPN 

^Peqv^N 

Equivalent  circular  nose  planform  area 

ft2 

45 

ECSPC 

^Peqv^C 

Equivalent  circular  centerbody  planform  area 

ft2 

46 

ECSPA 

(^Peqv^A 

Equivalent  circular  aftbody  planform  area 

ft2 

47 

ECSPT 

^Peqv 

Total  equivalent  circular  planform  area 

*  r 

ft2 

187 


I 


61-70 

71-80 

81-90 

91-100 

101-110 

111-120 


121-130 

131-140 

141-150 

151-160 

161-170 

171 

172 


Table  35  Finset  1  Geometry  Data  (Common  Block  GEOFS1) 


DEFINITION  OF  COMMON  BLOCK  GEOFS1 
(DUMP NAME  F1GM  .WRITE NAME  GEOFS1  ) 


VARIABLE 

NAME 


CBAR 

SWEPLE 

SWEP25 

SWEP50 

SWEP75 

SWEPTE 

XMGCLE 

XMGC25 

XMGC50 

XMGC75 

XMGCTE 

YMEANC 

ASPCT 

TAPER 

TO\rERC 

XCENT 

FINAR 

FINTPR 

FINSPN 

FINSP 

FINWET 

TCEFF 

FINXCG 

FINYCG 

SWAVLE 

SWAV25 

SWAV50 

SWAV75 

SWAVTE 

FINCBR 

DELTAY 

DELTAD 

XCENF 

XOVC 


ENGINEERING 

SYMBOL 


Q 

AL.E.j 

A.25ci 

A.50q 

A.75cj 

Aj.E.| 

XlEj 

X.25cj 

XsOq 

X;75ci 

XTEi 


X 

W1 

Sp 

Sw 

(t/c)e£f 

XMGC 

YmGC 

Al.E. 

A.25c 

A.50c 

A.75c 

At.E. 

c 

Ay 

5le 

Xc 

X/c 


DEFINITION 


orm  area  or  segment 
Mean  geometric  chord  of  segment 
Sweep  angle  at  L .  E .  of  segment 

Sweep  angle  at  c/4  of  segment 

Sweep  angle  at  c/2  of  segment 

Sweep  angle  at  3c/4  of  segment 

Sweep  angle  at  T.  E.  of  segment 

Distance  from  Cj  L.E.  to  M.G.C.  of  segment 

Distance  from  Cr  L.E.  to  .25  M.G.C.  of  segment 

Distance  from  Cj  L.E.  to  .50  M.G.C.  of  segment 

Distance  from  <*  L.E.  to  .75  M.G.C.  of  segment 

Distance  from  Cr  L.E.  to  T.E.  M.G.C.  of 
segment 

Semi-span  location  of  M.G.C.  of  segment 
Aspect  ratio  of  segment 
Taper  ratio  of  segment 

Thickness  to  chord  ratio  of  segment 
Planform  area  centroid  station  of  segment 

Panel  aspect  ratio  (overall) 

Panel  taper  ratio  (overall) 

Panel  exposed  semi-span  (overall) 

Panel  planform  area  (overall) 

Panel  wetted  area  (overall) 

Panel  effective  t/c  (overall) 

Panel  M.G.C.  L.E.  from  root  chord  L.E. 

Panel  lateral  M.G.C.  position  from  C.L. 

Panel  leading  edge  sweep,  effective 

Panel  c/4  sweep,  effective 

Panel  c/2  sweep,  effective 

Panel  3c/4  sweep,  effective 

Panel  T.E.  sweep,  effective 

Mean  geometric  chord 
Datcom  parameter  DY 

Panel  leading  edge  wedge  angle 

Panel  area  centroid  from  root  chord  L.E. 

Axial  position  of  section  max  thickness 


UNITS 


188 


Table  36  Finset  2  Geometry  Data  (Common  Block  GEOFS2) 


DEFINITION  OF  COMMON  BLOCK  GEOFS2 
(DUMP NAME  F2GM  .WRITE NAME  GEOFS2  ) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

iinnami 

Spj 

Planform  area  of1  segment 

ft2 

mmm 

CBAR 

Q 

Mean  geometric  chord  of  segment 

ft 

SWEPLE 

AL.E.j 

Sweep  angle  at  L .  E .  of  segment 

deg 

31-40 

SWEP25 

A,25ci 

Sweep  angle  at  c/4  of  segment 

deg 

41-50 

SWEP50 

A.50c; 

Sweep  angle  at  c/2  of  segment 

deg 

51-60 

SWEP75 

A.75cj 

Sweep  angle  at  3c/4  of  segment 

deg 

61-70 

SWEPTE 

AT.E.i 

Sweep  angle  at  T.  E.  of  segment 

deg 

71-80 

XMGCLE 

XlEj 

Distance  from  cr  L.E.  to  M.G.C.  of  segment 

ft 

81-90 

XMGC25 

X.25ci 

Distance  from  cr  L.E.  to  .25  M.G.C.  of  segment 

ft 

91-100 

XMGC50 

Distance  from  cr  L.E.  to  .50  M.G.C.  of  segment 

ft 

101-110 

XMGC75 

X.75CJ 

Xtei 

Distance  from  cr  L.E.  to  .75  M.G.C.  of  segment 

ft 

111-120 

XMGCTE 

Distance  from  cr  L.E.  to  T.E.  M.G.C.  of 
segment 

ft 

121-130 

YMEANC 

Yj 

Semi-span  location  of  M.G.C.  of  segment 

ft 

131-140 

ASPCT 

ARj 

Aspect  ratio  of  segment 

- 

141-150 

TAPER 

X.i 

t/Ci 

Taper  ratio  of  segment 

- 

151-160 

TOVERC 

Thickness  to  chord  ratio  of  segment 

- 

161-170 

XCENT 

Xci 

Planfoim  area  centroid  station  of  segment 

- 

171 

FINAR 

AR 

Panel  aspect  ratio  (overall) 

- 

172 

FINTPR 

X 

Panel  taper  ratio  (overall) 

- 

173 

FINSPN 

m 

Panel  exposed  semi-span  (overall) 

ft 

174 

FINSP 

Sp 

Panel  planform  area  (overall) 

ft2 

175 

FINWET 

sw 

Panel  wetted  area  (overall) 

ft2 

176 

TCEFF 

(t/c)eff 

Panel  effective  t/c  (overall) 

177 

FINXCG 

XMGC 

Panel  M.G.C.  L.E.  from  root  chord  L.E. 

ft 

178 

FINYCG 

YmGC 

Panel  lateral  M.G.C.  position  from  C.L. 

ft 

179 

SWAVLE 

Al.E. 

Panel  leading  edge  sweep,  effective 

deg 

180 

SWAV25 

A.  25c 

Panel  c/4  sweep,  effective 

deg 

181 

SWAV50 

A. 50c 

Panel  c/2  sweep,  effective 

deg 

182 

SWAV75 

A.75c 

Panel  3c/4  sweep,  effective 

deg 

183 

SWAVTE 

At.e. 

Panel  T.E.  sweep,  effective 

deg 

184 

FINCBR 

C 

Mean  geometric  chond 

ft 

185 

DELTAY 

Ay 

Datcom  parameter  DY 

- 

186 

DELTAD 

5lE 

Panel  leading  edge  wedge  angle 

deg 

187 

XCENF 

Xc 

Panel  area  centroid  from  root  chord  L.E. 

ft 

188 

XOVC 

X/c 

Axial  position  of  section  max  thickness 

189 


Table  37  Finset  3  Geometry  Data  (Common  Block  GEOFS3) 


DEFINITION  OF  COMMON  BLOCK  GEOFS3 
(DUMP NAME  F3GM  .WRITE NAME  GEOFS3  ) 


LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

Spi 

Planform  area  of  segment 

ft2 

CBAR 

Q 

Mean  geometric  chord  of  segment 

ft 

SWEPLE 

Al.E.j 

Sweep  angle  at  L .  E .  of  segment 

deg 

\m,Wm 

SWEP25 

A.25cj 

Sweep  angle  at  c/4  of  segment 

deg 

41-50 

SWEP50 

A.50c; 

Sweep  angle  at  c/2  of  segment 

deg 

51  60 

SWEP75 

A.75c; 

Sweep  angle  at  3c/4  of  segment 

deg 

61-70 

SWEPTE 

AT.E.j 

Sweep  angle  at  T.  E.  of  segment 

deg 

71-80 

XMGCLE 

XlEi 

Distance  from  cr  L.E.  to  M.G.C.  of  segment 

ft 

81-90 

XMGC25 

X.25ci 

Distance  from  cr  L.E.  to  .25  M.G.C.  of  segment 

ft 

91-100 

XMGC50 

X.50q 

Distance  from  Cf  L.E.  to  .50  M.G.C.  of  segment 

ft 

101-110 

XMGC75 

X.75q 

*TEj 

Distance  from  Cf  L.E.  to  .75  M.G.C.  of  segment 

ft 

111-120 

XMGCTE 

Distance  from  cr  L.E.  to  T.E.  M.G.C.  of 
segment 

ft 

121-130 

YMEANC 

Yj 

Semi-span  location  of  M.G.C.  of  segment 

ft 

131-140 

ASPCT 

AR, 

Aspect  ratio  of  segment 

- 

141-150 

TAPER 

Xi 

Taper  ratio  of  segment 

- 

151-160 

TOVERC 

t/Ci 

11110101685  to  chord  ratio  of  segment 

- 

161-170 

XCENT 

Xq 

AR 

Planform  area  centroid  station  of  segment 

- 

171 

FINAR 

Panel  aspect  ratio  ^overall) 

- 

172 

FINTPR 

A. 

Panel  taper  ratio  (overall) 

- 

173 

FINSPN 

IV2 

Panel  exposed  semi-span  (overall) 

ft 

174 

FINSP 

sp 

Panel  planform  area  (overall) 

fl2 

175 

FINWET 

Sw 

Panel  wetted  area  (overall) 

ft2 

176 

TCEFF 

(t/c)eff 

Panel  effective  t/c  (overall) 

- 

177 

FINXCG 

XmGC 

Panel  M.G.C.  L.E.  from  root  chord  L.E. 

ft 

178 

FINYCG 

YmGC 

Panel  lateral  M.G.C.  position  from  CJL. 

ft 

179 

SWAVLE 

Al.E. 

Panel  leading  edge  sweep,  effective 

deg 

180 

SWAV25 

A.25c 

Panel  c/4  sweep,  effective 

deg 

181 

SWAV50 

A.  50c 

Panel  c/2  sweep,  effective 

deg 

182 

SWAV75 

A.75c 

Panel  3c/4  sweep,  effective 

deg 

183 

SWAVTE 

At.E. 

Panel  T.E.  sweep,  effective 

deg 

184 

FINCBR 

c 

Mean  geometric  chord 

ft 

185 

DELTAY 

Ay 

Datcom  parameter  DY 

- 

186 

DELTAD 

8le 

Panel  leading  edge  wedge  angle 

deg 

187 

XCENF 

Xc 

Panel  area  centroid  from  root  chord  L.E. 

ft 

188 

XOVC 

X/c 

Axial  position  of  section  max  thickness 
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Table  38  F*nset  4  Geometry  Data  (Common  Block  GEOFS4) 


DEFINITION  OF  COMMON  BLOCK  GEOFS4 
(DUMP NAME  F4GM  .WRITE NAME  GEOFS4  ) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

imm^i 

Spj 

Planform  area  of  segment 

ft2 

11-20 

CBAR 

Cj 

Mean  geometric  chord  of  segment 

ft 

121-30 

SWEPLE 

AL.E.j 

Sweep  angle  at  L .  E  .  of  segment 

deg 

31-40 

SWEP25 

A.25ci 

Sweep  angle  at  c/4  of  segment 

deg 

41-50 

SWEP50 

A.50ci 

Sweep  angle  at  c/2  of  segment 

deg 

51-60 

SWEP75 

A.75ci 

Sweep  angle  at  3c/4  of  segment 

deg 

61-70 

SWEPTE 

AT.E.i 

Sweep  angle  at  T.  E.  of  segment 

deg 

71-80 

XMGCLE 

XlEj 

Distance  from  cr  L.E.  to  M.G.C.  of  segment 

ft 

81-90 

XMGC25 

X.25ct 

Distance  from  cr  L.E.  to  .25  M.G.C.  of  segment 

ft 

91-100 

XMGC50 

X^Oci 

Distance  from  L.E.  to  .50  M.G.C.  of  segment 

ft 

101-110 

XMGC75 

X.75cj 

Distance  from  cr  L.E.  to  .75  M.G.C.  of  segment 

ft 

111-120 

XMGCTE 

XTEi 

Distance  from  cr  L.E.  to  T.E.  M.G.C.  of 
segment 

ft 

121-130 

YMEANC 

Yj 

Semi-span  location  of  M.G.C.  of  segment 

ft 

131-140 

ASPCT 

ARj 

Xj 

Aspect  ratio  of  segment 

- 

141-150 

TAPER 

Taper  ratio  of  segment 

- 

151-160 

TOVERC 

t/q 

Thickness  to  chord  ratio  of  segment 

- 

161-170 

XCENT 

Xci 

AR 

Planform  area  centroid  station  of  segment 

- 

171 

FINAR 

Panel  aspect  ratio  (overall) 

- 

172 

F1NTPR 

X 

Panel  taper  ratio  (overall) 

- 

173 

FINSPN 

m 

Panel  exposed  semi-span  (overall) 

ft 

174 

FINSP 

sp 

Panel  planform  area  (overall) 

ft2 

175 

FIN  WET 

sw 

Panel  wetted  area  (overall) 

ft2 

176 

TCEFF 

(t/c  )eff 

Panel  effective  t/c  (overall) 

- 

177 

FINXCG 

Xmgc 

Panel  M.G.C.  L.E.  from  root  chord  L.E. 

ft 

178  i 

FINYCG 

YmGC 

Panel  lateral  M.G.C.  position  from  C.L. 

ft 

179 

SWAVLE 

Al.E. 

Panel  leading  edge  sweep,  effective 

deg 

I8u 

swa  vr 

A.25c 

Panel  c/*  sweep,  effective 

deg 

181 

SWAV50 

A. 50c 

Panel  c/2  sweep,  effective 

deg 

182 

SWAV75 

A.75c 

Panel  3c/4  sweep,  effective 

deg 

183 

SWAVTE 

at.e. 

Panel  T.E.  sweep,  effective 

deg 

184 

FINCBR 

C 

Mean  geometric  chord 

ft 

185 

DELTAY 

Ay 

Datcom  parameter  DY 

- 

186 

DELTAD 

5le 

Panel  leading  edge  wedge  angle 

deg 

187 

XCENF 

Xc 

Panel  area  centroid  from  root  chord  L.E. 

ft 

188 

XOVC 

x/c 

Axial  position  of  section  max  thickness 
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Table  39  Panel  Incidence  Data  (Common  Block  INCID) 


DEFINITION  OF  COMMON  BLOCK  INCID 

(DUMP  NAME 

.WRITE  NAME  INCID  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1-8  DELTA1  Si  Panel  deflections  for  fin  set  1  deg 

9-16  DELTA2  52  Panel  deflections  for  fin  set  2  deg 

17-24  DELTA3  S3  Panel  deflections  for  fin  set  3  deg 

25-32  DELTA4  5 4  Panel  deflections  for  fin  set  4  deg 

33-36  XHINGE  Xhl  Station  of  panel  hinge  line  for  fin  sets  1-4  ft 

37-40  SKEW  Ahl  Hinge  line  skew  for  each  panel,  +aft  deg 


Table  40  INLET  Namelist  Inputs  (Common  Block  INLETN) 


DEFINITION  OF  COMMON  BLOCK  INLETN 

(DUMP  NAME  INU 

.WRITE  NAME  INLETN  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

WIN 

Number  of  inlets 

2 

INTYPE 

Type  of  inlet 

0  =  2  dimensional  top  mounted 

- 

1  =  2  dimensional  side  mounted 

2  =  axisymmetric 

j 

XINLT 

Longitudinal  distance  from  nose  tip  to  inlet 
leading  edge 

ft 

A 

XDIV 

Longitudinal  distance  from  inlet  leading  edge  to 

ft 

diverter  leading  edge 

5 

HDIV 

Height  of  diverter  leading  edge 

ft 

6 

LD(V 

Length  of  diverter 

ft 

7-26 

PHI 

4> 

Inlet  roll  orientations 

deg 

27-31 

XI 

Inlet  longitudinal  positions  relative  to  inlet  ieading 
edge 

ft 

32-36 

HI 

Inlet  height  at  the  longitudinal  positions 

ft 

37-41 

WI 

Inlet  widths  at  the  longitudinal  positions 

ft 

42 

ICOVER 

if  .TRUE.  Inlets  are  covered 

- 

43 

RAMANG 

External  compression  ramp  angle 

deg 

44 

IAD 

if  .TRUE.  Inlet  additive  drag  is  calculated 

- 

45-64 

MFR 

Mass  flow  ratio  for  each  Mach  number 
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Table  41  Inlet  Incremental  Aerodynamics  (Common  Block  INLTD) 


DEFINITION  OF  COMMON  BLOCK  INLTD 
(DUMP NAME  INLD  .WRITE NAME  INLTD  ) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1-20 

CNINLT 

CNi 

Inlet  increment  for  normal  force  coefficient 

21-40 

CMINLT 

Cmi 

Inlet  increment  for  pitching  moment  coefficient 

- 

41-60 

CAINLT 

cm 

Inlet  increment  for  axial  force  coefficient 

- 

61-80 

CYINLT 

Cyi 

Inlet  increment  for  side  force  coefficient 

- 

81-100 

CLNILT 

Cnj 

Inlet  increment  for  yawing  moment  coefficient 
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Table  42  Namelist  Inputs  Names  (Common  Block  INPCON) 


DEFINITION  OF  COMMON  BLOCK  INPCON 
(DUMP NAME  .WRITE NAME  INPCON  ) 


LOCATION 


VARIABLE  ENGINEERING 
NAME  SYMBOL 


UNITS 


Position  of  first  letter  for  each  name 
Number  of  characters  in  each  name 
Namelist  names  valid  as  input 
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Table  43  Program  Execution  Logic  Flags  (Coir.mon  Block  LOGIC) 


DEFINITION  OF  COMMON  BLOCK  LOGIC 
(DUMP NAME  .WRITE NAME  LOGIC  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

Dump  case  dag 

2 

LDAMP 

Dynamic  derivatives  flag 

3 

LBUILD 

Configuration  buildup  flag 

4 

LNACA 

NACA  airfoil  designation  flag 

5 

LDERDG 

Degree  measure  flag 

6 

LDERRD 

Radian  measure  flag 

7 

LPART 

Partial  output  flag 

g 

LNAME 

Namelist  print  flag 

. 

9 

LPLOT 

Plot  flag 

10 

LFLT 

$FLTCON  input  flag 

. 

11 

LREFQ 

$REFQ  input  flag 

12 

LAXIS 

$AXIBOD  input  flag 

. 

13 

LF1N1 

$HNSET1  input  flag 

. 

14 

LFIN2 

SFINSET2  input  flag 

. 

15 

LFIN3 

$FINSET3  input  flag 

. 

16 

LFIN4 

$FINSET4  input  flag 

. 

17 

LDEFL 

SDEFLCT  input  flag 

18 

LTRIM 

$TRIM  input  flag 

19 

LDIMIN 

Inches  units  flag 

* 

20 

LDIMFT 

Feet  units  flag 

* 

21 

LDIMCM 

Centimeters  units  flag 

22 

LDIMM 

Meters  units  flag 

23 

LELLB 

SELLBOD  input  flag 

24 

LINLET 

$INLET  input  flag 

25 

LEXPR 

Experimental  data  input  flag 

26 

LICRMT 

Configuration  incrementing  flag 

. 

27 

LSPIN 

Body  magnus  coefficient  flag 

28 

LARBOD 

$  ARB  BOD  input  flag 

- 

NOTE:  All  flags  are  initially  false.  The  flag  is  set  if  the  array  element  is  true. 
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Table  44  Partial  Aerodynamic  Coefficients  (Common  Block  PAERO) 


DEFINITION  OF  COMMON  BLOCK  PAERO 

(DUMP  NAME 

,  WRITE  NAME  PAERO  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

21-24 

25-664 

665-1304 

1305-1944 

1945-2584 

2585-2664 

2665-2744 

2745-2824 

2825-2904 

2905-2984 

2985-3064 

3065-3144 

3145-3224 

3225-3304 

3305-3384 

3385-3464 

3465-3544 

3545-3624 

3625-3704 

3705-3784 

3785-3864 

3865-3944 

3945-4024 


AKWB 

AKKBW 

AKKWB 

XCPBW 

XCPWB 

DAQ 

WAQ 

CNW 

QLW 

CNWB 

CMWB 

CAWB 

CYWB 

CSNWB 

CSLWB 

CNBW 

CMBW 

CABW 

CYBW 

CSNBW 

CSLBW 

CN 

CM 

CA 

CY 

CSN 

CSL 


mg  carryover  on  body 
Body  carryover  on  wing 
Wing  incidence  carryover  cm  body 
Body  incidence  carryover  on  wing 
Body  c.p.  in  presence  of  wing 
Wing  c.p.  in  presence  of  body 
Change  in  equivalent  angle  of  attack 
Equivalent  angle  of  attack 
Wing  along  normal  for  coefficient 
Wing  alone  rolling  moment 
Wing  (body)  normal  force  coefficient 
Wing  (body)  pitching  moment  coefficient 
Wing  (body)  axial  force  coefficient 
Wing  (body)  side  force  coefficient 
Wing  (body)  yawing  moment  coefficient 
Wing  (body)  rolling  moment  coefficient 
Body  (wing)  normal  force  coefficient 
Body  (wing)  pitching  moment  coefficient 
Body  (wing)  axial  force  coefficient 
Body  (wing)  side  force  coefficient 
Body  (wing)  yawing  moment  coefficient 
Body  (wing)  rolling  moment  coefficient 
Total  normal  force  coefficient 
Total  pitching  moment  coefficient 
Total  axial  force  coefficient 
Total  side  force  coefficient 
Total  yawing  moment  coefficient 
Total  rolling  moment  coefficient 
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Table  45  REFQ  Namelist  Inputs  (Common  Block  REFQN) 


DEFINITION  OF  COMMON  BLOCK  REFQN 

(DUMP  NAME  REFQ 

,  WRITE  NAME  REFQN  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1 

SREF 

Sref 

Reference  area 

ft2 

2 

LREF 

(Lref^on 

Longitudinal  reference  length 

ft 

3 

LATREF 

(Lref)lat 

Lateral  reference  length 

ft 

4 

ROUGH 

h* 

Surface  roughness  height 

in 

5 

XCG 

XC.G. 

Longitudinal  center  of  gravity 

ft 

6 

ZCG 

^C.G. 

Vertical  center  of  gravity 

ft 

7 

SCALE 

Vehicle  scale  factor 

- 

8 

BLAYER 

Boundary  layer  type 

0.  =  turbulent 

- 

1.  =  natural  transition 

9 

RHR 

Surface  roughness  height  rating 
(ROUGH  =3x  KT6 (RHR)] 
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Table  46  Body  and  Finset  1  Ideal  Output  Matrix  (I.O.M.)  (Common  Block  SB1) 


DEFINITION  OF  COMMON  BLOCK  SB1 
(DUMP NAME  SB1  .WRITE NAME  SB1  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

CN1  ~~~ 

Normal  force  vs.  a 

CM1 

Cm 

Pitching  moment  vs.  a 

41-60 

CA1 

CA 

Axial  force. vs.  a 

61-80 

CY1 

Cy 

Side  force  vs.  a 

81-100 

CSN1 

Cn 

Yawing  moment  vs.  a 

. 

101-120 

CSL1 

Cl 

Rolling  moment  vs.  a 

121-140 

CNA1 

Cn« 

CNa  vs-  a 

1/deg 

141-160 

CMA1 

Cmot 

Cma  vs.  a 

1/deg 

161-180 

CYB1 

CYp 

Cyp  vs.  a 

1/deg 

181-200 

CLNB1 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

CLLB1 

Clp 

Cjp  vs.  a 

1/deg 

Table  47  Body  and  Finset  1,2 I.O.M.  (Common  Block  SB12) 


DEFINITION  OF  COMMON  BLOCK  SB12 

(DUMP  NAME  SB  12 

.WRITE NAME  SB12  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

- 1 

i 


-40 
-60 
61-80 
81-100 
101-120 
121-140 
141-160 
161-180 
181-200 
201-220 


1 

CM12 

CA12 

CY12 


CNA12 

CMA12 

CYB12 

CLNB12 

CLLB12 


Cm 

CA 

Cy 

Cn 

Cl 

cNa 

Cm-, 


Normal  force  vs.  a 
Pitching  moment  vs.  a 
Axial  force,  vs.  a 
Side  force  vs.  a 
Yawing  moment  vs.  a 
Rolling  moment  vs.  a 
CNa  vs-  a 
Cma ® 

Cyp  vs.  a 
Cnp  vs.  a 
Cip  vs.  a 
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Table  48  Body  and  Finset  1,2,3  I.O.M.  (Common  Block  SB123) 


DEFINITION  OF  COMMON  BLOCK  SB  123 
(DUMP NAME  SB  13  .WRITE NAME  SB  123  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINmON 

UNITS 

Cn 

Normal  force  vs.  a 

B 

CM13 

Cm 

Pitching  moment  vs.  a 

41-60 

CA13 

CA 

Axial  force,  vs.  a 

61-80 

CY13 

Cy 

Side  force  vs.  a 

81-100 

CSN13 

Cn 

Yawing  moment  vs.  a 

101-120 

CSL13 

Ci 

Rolling  moment  vs.  a 

121-140 

CNA13 

CNa 

Cnq  vs.  a 

1/deg 

141-160 

CMA13 

Cma 

Cma  vs.  a 

1/deg 

161-180 

CYB13 

Cyp 

Cyp  vs.  a 

1/deg 

181-200 

CLNB13 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

CLLB13 

Qp 

Qp  vs.  a 

1/deg 

Table  49  Body  and  Finset  1 ,2,3/4  I.O.M.  (Common  Block  SB1234) 


DEFINITION  OF  COMMON  BLOCK  SB1234 
(DUMP NAME  SB  14  .WRITE NAME  SB1234  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

ON  14 

Normal  force  vs.  a 

21-40 

CM  14 

Cm 

Pitching  moment  vs.  a 

41-60 

CA14 

CA 

Axial  force. vs.  a 

61-80 

CY14 

Cy 

Side  force  vs.  a 

81-100 

1  ..  H 

c„ 

Yawing  moment  vs.  a 

101-120 

Cl 

Rolling  moment  vs.  a 

121-140 

CNA14 

CNa 

CNa  vs.  a 

1/deg 

141-160 

CMA14 

Cma 

Cma  vs*  a 

1/deg 

161-180 

CYB14 

cYp 

Cyp  vs.  a 

1/deg 

181-200 

CLNB14 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

CLLB14 

Cip 

Cip  vs.  a 

1/deg 
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Table  50  Body  I.O.M.  (Common  Block  SBODY) 


DEFINITION  OF  COMMON  BLOCK  SBODY 
(DUMP NAME  SBOD  .WRITE NAME  SBODY  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

■■SOHH 

Normal  force  vs.  a 

_ 

CM 

Cm 

Pitching  moment  vs.  a 

41-60 

CA 

CA 

Axial  force,  vs.  a 

_ 

61-80 

CY 

Cy 

Side  force  vs.  a 

. 

81-100 

CLN 

Cn 

Yawing  moment  vs.  a 

_ 

101-120 

CLL 

Q 

Rolling  moment  vs.  a 

_ 

121-140 

CNA 

cNa 

CNa  vs-  a 

1/deg 

141-160 

CMA 

Cma 

Cma ® 

1/deg 

161-180 

CYB 

Cyp 

Cyp  vs.  a 

1/deg 

181-200 

CLNB 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

CLUB 

Qp 

Qp  vs.  a 

1/deg 
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Table  51  Finset  1 1.O.M.  (Common  Block  SFIN1) 


DEFINITION  OF  COMMON  BLOCK  SFIN1 
(DUMP NAME  SF1  .WRITE NAME  SFIN1  ) 

VARIABLE 

ENGINEERING 

LOCATION 

SYMBOL 

DEFINITION 

UNITS 

■ 

mmsoHH 

Nonnal  force  vs.  a 

FI  CM 

Cm 

Pitching  moment  vs.  a 

_ 

41-60 

F1CA 

Ca 

Axial  force.vs.  a 

. 

61-80 

F1CY 

Cy 

Side  force  vs.  a 

. 

81-100 

F1CLN 

Cn 

Yawing  moment  vs.  a 

_ 

101-120 

F1CLL 

Q 

Rolling  moment  vs.  a 

_ 

121-140 

F1CNA 

CNa 

Cnq  vs*  ® 

1/deg 

141-160 

F1CMA 

Cma 

Cma  vs.  a 

1/deg 

161-180 

F1CYB 

CYp 

Cyp  vs.  a 

1/deg 

181-200 

F1CLNB 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

F1CLLB 

Clp 

Cip  vs.  a 

1/deg 

204 


Table  52  Finset  2  I.O.M.  (Common  Block  SFIN2) 


DEFINITION  OF  COMMON  BLOCK  SFIN2 
(DUMP NAME  SF2  .WRITE NAME  SFIN2  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

Cn 

Nonnal  force  vs.  a 

- 

F2CM 

Cm 

Pitching  moment  vs.  a 

- 

41-60 

F2CA 

Ca 

Axial  force.vs.  a 

- 

61-80 

F2CY 

Cy 

Side  force  vs.  a 

- 

81-100 

F2CLN 

Cn 

Yawing  moment  vs.  a 

- 

101-120 

F2CLL 

Ci 

Rolling  moment  vs.  a 

- 

121-140 

F2CNA 

CNa 

Cnu  vs.  a 

1/deg 

141-160 

F2CMA 

Cma 

Cma  ® 

1/deg 

161-180 

F2CYB 

Cyp 

Cyp  vs.  a 

1/deg 

181-200 

F2CLNB 

Clip 

Cnp  vs.  a 

1/deg 

201-220 

F2CLLB 

Qp 

Cjp  vs.  a 

1/deg 
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Table  53  Finset  3  I.O.M.  (Common  Block  SFIN3) 


DEFINITION  OF  COMMON  BLOCK  SF1N3 
(DUMP NAME  SF3  .WRITE NAME  SFIN3  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

Cn 

Normal  force  vs.  a 

F3CM 

Cm 

Pitching  moment  vs.  a 

F3CA 

CA 

Axial  force.vs.  a 

61-80 

F3CY 

Cy 

Side  force  vs.  a 

81-100 

F3CLN 

Cn 

Yawing  moment  vs.  a 

_ 

101-120 

F3CLL 

Cl 

Rolling  moment  vs.  a 

121-140 

F3CNA 

CNa 

CNa  vs.  a 

l/deg 

141-160 

F3CMA 

Cma 

Cma  vs.  a 

l/deg 

161-180 

F3CYB 

CYp 

Cyp  vs.  a 

l/deg 

181-200 

F3CLNB 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

F3CLLB 

Cip 

Cip  vs.  a 

l/deg 
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Table  54  Finset  4  I.O.M.  (Common  Block  SFIN4) 


DEFINITION  OF  COMMON  BLOCK  SFIN4 
(DUMP NAME  SF4  .WRITE NAME  SFIN4  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

Cn 

Normal  force  vs.  a 

iSlSlSw 

F4CM 

Cm 

Pitching  moment  vs.  a 

- 

F4CA 

Ca 

Axial  force,  vs.  a 

_ 

6-80 

F4CY 

Cy 

Side  force  vs.  a 

_ 

81-100 

F4CLN 

Cn 

Yawing  moment  vs.  a 

- 

101-120 

F4CLL 

Cl 

Rolling  moment  vs.  a 

- 

121-140 

F4CNA 

CNa 

Cnq  vs.  a 

1/deg 

141-160 

F4CMA 

Cma 

Cma  vs.  (X 

1/deg 

161-180 

F4CYB 

Cyp 

Cyp  vs.  a 

1/deg 

181-200 

F4CLNB 

Cnp 

Cnp  vs.  a 

1/deg 

201-220 

F4CLLB 

Clp 

Cip  vs.  a 

1/deg 
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Table  55  Methodology  Flags  (Common  Block  THERY) 


DEFINITION  OF  COMMON  BLOCK  THERY 
(DUMP  NAME  .WRITE  NAME  THERY  ) 


UNITS 


.TRUE,  if  SOSE  to  be  used 
.TRUE,  if  pressure  data  to  be  output 
.TRUE,  if  Hybrid  method  to  be  used 
.TRUE,  if  Newtonian  theory  to  be  used 
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Table  56  Configuration  Total  Attitude  (Common  Block  TOTALC) 


DEFINITION  OF  COMMON  BLOCK  TOTALC 
(DUMP NAME  FLCT  .WRITE NAME  TOTALC  ) 


LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1-20 

BALPHA 

a 

Body  axis  angle  of  attack 

deg 

21-40 

BBETA 

P 

Body  axis  sideslip  angle 

deg 

41-60 

BPHI 

4> 

Roll  angle 

deg 

61-80 

ALPT0T 

a' 

Total  angle  of  attack 

deg 
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Table  57  TRACE  Subroutine  Variables  (Common  Block  TRACE) 


DEFINITION  OF  COMMON  BLOCK  TRACE 

(DUMP  NAME 

.WRITE  NAME  TRACE  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1  LEVEL  Program  level 

. .  Routine  names  by  level  (2  locations  per  name,  10 

2-21  IROUTN  ]evels)  J 
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Table  58  Trimmed  Aerodynamic  Data  Arrays  (Common  Block  TRIMD) 


DEFINITION  OF  COMMON  BLOCK  TRIMD 
(DUMP  NAME  ,  WRITE  NAME  TRIMD  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1-20 

DELTRM 

^TRTM 

^TRIM  vs.  a 

deg 

21-40 

CNTRM 

CNTRIM 

CNTRIM  vs.  a 

- 

41-60 

CATRM 

CAtRIM 

CATRIM  vs.  a 

- 

61-80 

CYTRM 

CYTRIM 

CYTRIM  vs.  a 

- 

81-100 

CLNTRM 

CnTRIM 

CnTRIM  vs.  a 

- 

101-120 

CLLTRM 

ClTRIM 

ClTRIM  vs.  a 
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Table  59  Trim  Inputs  (Common  Block  TRIMIN) 


DEFINITION  OF  COMMON  BLOCK  TRIMIN 
(DUMP  NAME  .WRITE  NAME  TRIMIN  ) 

LOCATION 

VARIABLE 

NAME 

ENGINEERING 

SYMBOL 

DEFINITION 

UNITS 

1 

SET 

Fin  set  for  trimming 

- 

2 

PANL1 

.TRUE,  if  used  for  trim 

- 

3 

PANL2 

.TRUE,  if  used  for  trim 

- 

4 

PANL3 

.TRUE,  if  used  for  trim 

- 

5 

PANL4 

.TRUE,  if  used  for  trim 

- 

6 

PANL5 

.TRUE,  if  used  for  trim 

- 

7 

PANL6 

.TRUE,  if  used  for  trim 

- 

8 

PANL7 

.TRUE,  if  used  for  trim 

- 

9 

PANL8 

.TRUE,  if  used  for  trim 

- 

10 

DELMIN 

Smin 

Minimum  deflection  angle 

deg 

11 

DELMAX 

Smax 

Maximum  deflection  angle 

deg 

12-19 

ASYM 

.TRUE,  if  panel  is  deflected  asymmetric 
(reverses  sign  convention) 
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Table  60  Untrimmed  Aerodynamic  Data  Arrays  (Common  Block  UTRIMD) 


DEFINITION  OF  COMMON  BLOCK  UTRIMD 
(DUMP  NAME  ,  WRITE  NAME  UTRIMD  ) 

VARIABLE 

ENGINEERING 

LOCATION 

NAME 

SYMBOL 

DEFINITION 

UNITS 

1-20 

CN 

Cn 

Cn  vs.  a  for  8i 

- 

21-40 

Cn  vs.  a  for  62 

- 

41-60 

Cn  vs.  a  for  83 

- 

61-80 

Cn  vs.  a  for  84 

- 

81-100 

Cn  vs.  a  for  85 

- 

101-120 

Cn  vs.  a  for  86 

- 

121-140 

Cn  vs.  a  for  87 

- 

141-160 

Cn  vs.  a  for  8s 

- 

161-180 

Cn  vs.  a  for  89 

- 

181-200 

Cn  vs.  a  for  810 

- 

201-400 

CM 

Cm 

Cm  vs.  a,8  (see  Cn  pattern) 

- 

401-600 

CA 

cA 

Ca  vs.  0,8  (see  Cn  pattern) 

- 

601-800 

CY 

Cy 

Cy  vs.  0,5  (see  Cn  pattern) 

- 

801-1000 

CLN 

Cn 

Cn  vs.  0,8  (see  Cn  pattern) 

- 

1001-1200 

CLL 

Q 

Ci  vs.  0,8  (see  Cn  pattern) 
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A.  EXAMPLE  PROBLEMS 


This  appendix  presents  two  example  missile  configurations  that  have 
been  synthesized  as  Missile  Datcom  input  files.  These  examples  may  be  used 
as  a  model  for  the  setup  of  inputs  of  similar  configurations. 

The  first  example  is  a  simple  tangent  ogive  nose-cylinder  circular  body. 
It  has  a  planar  wing  (two  panels)  and  a  cruciform  set  of  tails  orientated  in  the 
"plus"  configuration.  There  are  two  input  cases  for  this  problem.  The  first 
generates  output  consistent  with  the  inclusion  of  the  PART  control  card.  The 
second  creates  trimmed  output  (note  that  trim  partial  output  has  been 
requested  with  the  PRINT  AERO  TRIM  control  card). 

The  second  problem  is  a  body-tail-inlet  configuration.  Three  Mach 
numbers  are  requested,  one  subsonic,  one  transonic,  and  one  supersonic. 
Although  all  three  Mach  numbers  could  be  run  in  a  single  case,  they  have 
been  divided  into  three  separate  cases  to  illustrate  the  "SAVE"  feature  as  well 
as  selecting  particular  output  for  illustration. 
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A.l  EXAMPLE  PROBLEM  1 


The  first  example  problem  is  shown  in  Figure  A-l.  It  is  comprised  of  a 
3-caliber  tangent  ogive  nose  attached  to  a  cylindrical  body;  a  triangular 
monoplane  set  of  wings;  and  a  crudform  set  of  tails  orientated  in  the  "plus" 
position.  The  first  case  is  a  simple  angle  of  attack  sweep;  component  buildup 
data  and  partial  output  are  requested.  The  second  case  is  a  trim  of  the 
configuration  using  the  two  horizontal  tail  surfaces.  The  inputs  are  shown  in 
Figure  A-2. 
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Figure  A-2  Example  Problem  1  Input 
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A.2  EXAMPLE  PROBLEM  2 


The  configuration  for  this  example  is  sketched  in  Figures  A-3  and  A-4. 
The  figure  is  a  modified  copy  of  the  wind  tunnel  model  drawings  from 
NASA  Technical  Memorandum  84557.  The  model  definition  in  these  figures 
is  representative  of  the  detail  normally  found  on  design  drawings. 

This  example  has  been  divided  into  subsonic,  transonic,  and 
supersonic  cases.  Each  case  is  run  for  one  Mach  number.  Although  all  three 
Mach  numbers  could  have  been  run  in  one  case  they  were  run  separately  to 
demonstrate  the  SAVE  capability. 

This  example  provides  a  check  case  for  the  inlet  option.  It  can  be  used 
by  the  user  to  make  sure  that  he  understands  the  inputs.  Figure  A-5  shows 
the  inputs  required  to  run  this  example. 
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Figure  A-3  Example  Problem  2  Configuration 


220 


Figure  A-5  Example  Problem  2  Input 
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Figure  A-5  Example  Problem  2  Input  (Continued) 
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B.  PLOT  FILE  FORMAT 


When  the  PLOT  control  card  is  used,  a  formatted  data  file  is  written  to 
unit  3.  The  first  line  in  the  file  will  be  the  word  MISDAT.  Lines  with  the 
following  format  will  be  repeated  for  each  run: 


LINE 

COLUMN 

1  CONTENT 

i 

1-4 

Wo-d  RUN 

5-7 

For  each  case,  a  sequential  run  number 
beginning  with  one  (1)  and  incrementing 
by  one  (1) 

S-10 

Number  of  angles  of  attack 

11-14 

word  CASE 

15-17 

A  sequential  run  number  beginning  with 
one  (1)  and  incrementing  by  one  (1) 

2 

1-10 

Mach  Number 

11-22 

Reynolds  Number 

23-32 

Aerodynamic  roll  angle,  phi 

33-34 

Units  system:  1-inches,  F-feet,  C- 
centimeters,  M-meters 

35-41 

I  D.  Code  (see  below) 

3 

1-10 

Reference  Area 

11-20 

Longitudinal  Reference  Length 

21-30 

Longitudinal  C.G.  Location  (X  CR) 

31-40 

Lateral  Reference  Length 

41-50 

Vertical  C.G.  Location 

4 

1-10 

Angle  of  Attack 

(Repeated  for 

11-20 

CN  or  “TRIM 

each  angle  of 

21-30 

Cm  or  CNTRIM 

attack) 

31-40 

Ca  orCATRIM 

41-50 

Cy  or  CyrRIM 

1-60 

Cn  or  CnTRlM 

1-70 

Ci  or  CjTRlM 

LAST 

1  |  Character  R  (denotes  end  of  data  for  set)  | 

The  I.D.  field  identifies  the  configuration  for  which  the  data  is 
applicable. 


I.D.  Code 

Configuration  /  Condition 

B 

Body  Alone 

FI 

Most  forward  fin  set  or  fin  set  #1 

F2 

Second  most  forward  fin  set  or  fin  set  #2 

F3 

Third  most  forward  fin  set  or  fin  set  #3 

F4 

Fourth  most  forward  fin  set  or  fin  set  #4 

BF1 

Body  plus  most  forward  fin  set 

BF12 

Body  plus  two  most  forward  fin  sets 

BF13 

Body  plus  three  most  forward  fin  sets 

BF14 

Body  plus  four  most  forward  fin  sets 

TRIMMED 

Trimmed  results 
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Component  buildup  data  is  written  to  the  file  if  the  BUILD  control  card 
is  used.  If  the  TRIM  control  card  is  used,  both  the  trimmed  and  untrimmed 
results  are  output.  Otherwise,  only  the  final  configuration  data  is  output.  It  is 
important  to  note  that  output  via  WRITE  control  cards  is  also  written  to 
unit  3. 
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